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PARALLACTIC MEASUREMENTS ON AURORAE OVER 
CHURCHILL, CANADA! 


D. J. McCEWEN? AND R. MONTALBETTI? 


ABSTRACT 


Height measurements on the lower edges of homogeneous and rayed arcs and 
bands are reported from 404 photographs taken of aurorae over Churchill (58.8° 
N., 94.2° W.) during the winter of 1957-58 as part of the I.G.Y. programme. The 
frequency distribution in height of 1975 points measured shows a double peak 
at 101 and 107 km, and a mean height of 104.8 km. The chief factor affecting 
height of aurora appears to be that of intensity; the mean height of weak aurorae 
was found to be 107.7 km, while that of strong aurorae was 100.1 km. No differ- 
ence in mean height of aurorae before and after midnight was observed. Tidal 
effects in the upper atmosphere did not appear to affect auroral height greatly. 
The mean heights of aurorae measured during high and low tide conditions 
showed no appreciable difference. A mean orientation of 105.5° magnetic azimuth 
was measured for auroral arcs. Auroral motions with speeds up to 2500 meters 
per second were observed. 


1, INTRODUCTION 


The Defence Research Northern Laboratory and the University of Saskat- 
chewan, Saskatoon, undertook to provide height measurements and spectral 
observations in support of the U.S.A. rocket firings into aurorae at Churchill 
during the International Geophysical Year. Two stations were established 
along the railroad to Churchill for this purpose, one at Belcher (57.9° N., 
94.2° W.) and the other at O’ Day (57.6° N., 94.2° W.). The reasons for choosing 
these particular sites have been discussed by Montalbetti and Lutz (1956). 
Advantage was taken of the available facilities to make a more detailed study 
of auroral heights and related characteristics than was required for the rocket 
firings. The results of the study are reported in this paper. 

This is the only substantial set of measurements on auroral heights that have 
been made in North America along the zone of maximum auroral occurrence. 
Previous measurements in comparable numbers have been made by Currie 
(1955) at Chesterfield, N.W.T., about 350 km inside the auroral zone, and 
by Fuller (1933) at College, Alaska, about the same distance outside the 
auroral zone. As far as measurements in Canada are concerned, the Churchill 
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measurements are the only ones that have been made close to a maximum of 
solar activity. 


2. EXPERIMENTAL DETAILS 


Cameras for the parallactic photography used 35-mm film, and were of a 
type which permitted automatic adjustment for duration of exposure and 
number of exposures per minute. They were equipped with Dallmeyer F/1.5 
objectives. They were automatically controlled by timing units run from 
frequency standards. Radio communication between the two stations ensured 
that the exposures were always synchronized. Exposures were normally of 
5-second duration, and were made at the rate of 1 or 10 per minute as desired. 

Since information on aurorae over Churchill was the primary purpose for 
operating the cameras, they were mounted on concrete piers, and were operated 
usually with a fixed orientation that included the anticipated levels of aurorae 
over Churchill. The distance between the cameras was 39.47 km (somewhat 
greater than customarily used for measurements of this type); the direction of 
the line joining the cameras was approximately parallel to the geographic 
meridians (see Fig. 5). The long base line and the high viewing angle combined 
to give parallactic displacements of 10° or more for aurorae above Churchill. 
This, together with the favorable direction of the base line relative to the 
orientation.of auroral arcs, suggests that the measured heights are of greater 
accuracy than is normally achieved. One disadvantage of the fixed orientation 
for the cameras was the limited portion of an auroral display on which height 
measurements could be made. 

The angular field covered by each camera was about 30X40°. To ensure 
that the two cameras viewed the same aurora, the closer camera was directed 
10° higher than the farther one giving an effective field of 40 40°. 

The method used in measuring height and position of aurorae from the 
photozraphs was essentially the one developed by Stérmer and discussed by 
Harang (1951). Graphical aids called ‘‘nets’’ are used. They are projected 
portions on a plane of the circles of constant declination and hour angle on 
the celestial sphere. Such ‘‘nets’’ can be used to determine altitude and 
azimuth or base distance and parallax equally well, and it is customary to use 
the nets to determine the necessary co-ordinates for auroral points in the 
photographs. This is done by laying a tracing of the projections of the auroral 
form on the nets in turn and reading off the co-ordinates. 

The fixed orientation for the cameras permitted some modifications in this 
method. A large net was prepared and a positive print of this original was made 
on sheet film. This transparent print was placed on the original and oriented 
correctly by making a spherical co-ordinate transformation so that both 
horizontal and base pole co-ordinates were represented. Thus when the two 
auroral photographs were projected on this double net, which also served as a 
projection screen, all the necessary co-ordinates and the parallax for any 
particular auroral point could be read off simultaneously. This technique is 
illustrated in Fig. 1. 

A glass reticle, on which cross hairs had been etched, in each camera provided 
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Fic. 1. Simultaneous photographs of an auroral arc taken at the two stations. They are 
projected on a double net, with the stars superimposed, and show the parallactic displacement. 





fixed azimuth and elevation points on each exposure. These horizontal co- 
ordinates were calculated using known stars appearing in the photographs, 
and they served to orient the pictures correctly on the double net for as long 
as the cameras remained in that position. 


3. RESULTS 

The stations were in operation a total of 81 nights in the period from August 
1957 to April 1958. On 47 of these nights aurora was observed and photo- 
graphed; overcast skies accounted for practically all of the other nights. 
Cloud conditions persisted during the fall and spring months in particular 
and thus most of the successful pictures were taken during the middle of 
December and February, both periods of pronounced auroral activity. All 
of the photographs of aurorae occurring in the Churchill area which were 
suitable for measurement, a total of 404 pairs, have been reduced yielding 
approximately 2000 heights of the lower edges of homogeneous and rayed 
arcs and bands. Since the exposure time was fixed, results are confined to 
aurorae of intensity 1 or greater (based on the scale of brightness given by 
Hunten, Roach, and Chamberlain 1956). 

The frequency distribution of the measured heights is shown in Fig. 2 (0). 
The same characteristic double peak as has been observed by Stérmer (1955) 
and others is evident, and the mean height of 104.8 km is very similar to that 
obtained by comparable studies elsewhere. An analysis of the relation between 
auroral intensity and height was made by dividing the pictures into three 
groups of weak, medium, and strong aurorae. Figure 2 (a) shows the frequency 
distribution of height of weak and strong aurorae, and suggests a probable 
reason for the double peak in Fig. 2 (6). 
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Fic. 2. Height distribution of lower edges of homogeneous and rayed arcs and bands: 
(a) showing weak and strong aurorae only and (6) showing all measured points. 


Egedal (1929) has suggested that the double maxima, usually observed 
about 6 km apart, can be explained by tidal effects. While an analysis of the 
frequency curves for low and high tide conditions shown in Fig. 3 (5) indicates 
a difference in mean height of 0.6 km, when the curves are redrawn using 
only medium aurorae to eliminate any intensity effect this difference is re- 
moved entirely as seen in Fig. 3 (a). This indicates that the tidal effect may 
not be as pronounced in this region. 

To investigate a possible change in height of aurorae during the course of 
the night as has been discussed by Stérmer (1955), the distributions in height 
of aurorae occurring before and after midnight were plotted, and are shown 
in Fig. 4 (6). A study of the hourly variation was not feasible because of 
insufficient measurements distributed throughout the night. Since stronger 
aurorae tended to predominate before midnight this effect was removed by 
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Fic. 3. Distribution of heights during low and high tide conditions for (a) aurorae of 
medium intensity and (0) all aurorae. 


plotting only heights of weak aurorae in Fig. 4 (a). No significant change in 
height is observable. This is in agreement with measurements reported by 
Currie (1955). 

During many auroral displays much rapid motion is observed. Since ex- 
posures were made at the rate of 10 per minute during such displays, an 
accurate measure of velocity could be made by plotting the auroral positions 
for a sequence of pictures on a map, as illustrated in Fig. 5. In this case a 
northward motion with horizontal speed up to 1000 meters/second was 
measured. A number of sequences of this type, where motion appeared to be 
normal to the direction of the arc, were plotted. Results are given in Table I. 
It should be noted that these are short term velocities rather than average 
drift motion such as reported by Kim and Currie (1958). No marked correlation 
between time of night and direction of motion was observed. 

A plot of position of aurorae permitted also a measure of orientation of arcs. 
An average orientation was determined for each display from the results of all 
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Fic. 4. Distribution of heights measured before and after midnight of (a) weak aurorae 
and (6) all aurorae. 


TABLE I 
SPEEDS OF NORTH-SOUTH MOTIONS OF AURORAL ARCS AT CHURCHILL 


No. Max. speed, 





Direction measurements meters/second meters/second 
Northward 34 2100 780 
Southward 47 2500 930 





individual measurements made during that display. The orientations varied 
from a magnetic azimuth, measured clockwise from the north, of 70° to a 
maximum of :43°; this wide spread must be attributed to the fact that the 
majority of the displays were active rather than the quiet type to which 
orientation measurements have usually been confined. It is interesting to note, 
however, that the average magnetic azimuth of the 41 displays measured was 
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Fic. 5. A map of the Churchill area indicating the location of the observation stations and 
showing in dashed lines the geomagnetic polar distances and meridians. The change in position 
of an auroral arc, photographed on Feb. 13, 1958, over an interval of 30 seconds is illustrated. 


105.5°. This agrees very well with values obtained at Fort Rae and Chester- 
field during the Second International Polar Year, as given by Currie and 
Jones (1941). 

It is considered that the foregoing measurements are of sufficient number 
to indicate generally the distribution of heights and the factors affecting 
height, as well as the orientation and movements, of auroral arcs and bands 
in the Churchill area. The analysis of the photographs which have been taken 
in other directions from the two auroral stations and the further data anti- 
cipated from the continuing operations should provide more detailed in- 
formation, particularly of diurnal and seasonal changes of auroral heights, 
orientations, and movements. 
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INVESTIGATION OF THE DRIFTS OF THE EFFECTIVE POINT 
OF RADIO REFLECTION ALONG A METEOR TRAIN! 


M. SRIRAMA RAO? AND R. L. ARMSTRONG? 


ABSTRACT 


The concept of the effective point of radio reflection from a meteor train has 
been given on the basis of atmospheric turbulence, on a vertical scale of the order 
of 100 m to6 km, in the M region. Some experimental evidence has been provided 
to support the postulate of its drifting along the train, using meteoric body 
doppler radar records, taken on 30.02 Mc/s at South Gloucester, during the 
Geminid shower periods of 1948-50. The velocities (V;) of the above drifts in the 
case of 90 observations have been calculated, using the ranges of the observed 
meteors as obtained from the range-time records taken simultaneously on pulsed 
radar on 32.7 Mc/s at the Metcalfe Road field station (7.5 km distant). It is 
found that these velocities tend to have higher values in the case of shorter 
echo durations and vice versa. Theoretical interpretation of the observed results 
has been attempted. Reasonable assumptions of the ionization distribution along 
a meteor train and of diffusion rates at different levels in the M region have been 
made to derive the echo durations from different portions of a meteor train. 
Variation of these echo durations with position on the train has been taken 
into account to calculate the theoretical curves of vertical components of V, vs. 
total echo durations, in the case of a Geminid shower. The effect of turbulence 
on echo durations has been taken into account to explain the observed results 
successfully. 


INTRODUCTION 


The concept of the effective or average point of radio reflection from a 
meteor train has been given by one of the authors in a previous paper (Rao 
1958). It has been shown in that paper (let us refer to this paper as paper 1, 
hereafter) that atmospheric turbulence on a vertical scale of anywhere between 
100 m and 6 km would be responsible in the breaking up of the effective length 
of a meteor train into relatively smaller segments. Each of these segments 
would then reflect radio waves towards the observer. The result of a combina- 
tion of all these echoes has been visualized as equivalent to an echo from an 
effective point of reflection on the meteor train. This may in general coincide 
approximately with the mid-point of the effective length of the train under 
consideration. It has also been suggested that, owing to rapid diffusion of 
ionization, the disappearance of ionization along the train progresses rapidly 
from one end inwards. Considering the duration of echoes from the various 
segments, this is a minimum at the ends, and gradually increases inwards 
from either side, reaching a maximum somewhere on the train. This would 
give rise to the phenomenon of drifting of the end point of reflection, and 
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hence also of the effective point of reflection, towards the maximum echo 
duration level. The average velocities of drift of the effective points of re- 
flection have been found, in paper 1, to be high compared to the wind velocities 
in the M region. It has also been found that these velocities tend to be smaller 
in the case of meteor trains of higher duration. It is the object of the present 
investigation to verify this result with more observations and also to discuss 
the problem theoretically. Theoretical velocities of drift of the end point of 
radio reflection along the train may be calculated by assuming a certain 
theoretical distribution of ionization along the train, and certain plausible 
diffusion rates at different heights above the ground. 

Ideally, a meteor train may be considered as divided into two main sections, 
one above the maximum echo duration level and the other below it. Each 
of these two sections may be considered to have a separate effective point of 
reflection. We may designate these as the upper and lower effective points 
of reflection respectively of a meteor train. Assuming that the effective point 
of reflection, in each case, is approximately at the mid-point of the effective 
length of the train under consideration, it can easily be seen that the effective 
point drifts along the train with nearly half the speed of drift of the corre- 
sponding end point towards the maximum echo duration level. 

Let us first proceed to give some experimental evidences for the movement 
of the effective point of reflection, obtained from the analysis of a number of 
body doppler radar records taken at South Gloucester, on 30.02 Mc/s, during 
four nights in the December months of 1948-50. These nights have been chosen 
since the Geminid shower is known to be active at this time. 


EXPERIMENTAL OBSERVATIONS 


The experimental method of investigation of the average velocities of drift, 
V,, of the effective point of reflection of CW radio waves along a meteor train 
has already been given in paper 1. This method essentially consists in analyzing 
the body doppler radar records to determine the variation of the average 
doppler beat frequency, fy, with time. Thus, all the fast doppler records taken 
at the South Gloucester station, using CW transmission at 30.02 Mc/s from 
the Metcalfe Road field station in Ottawa (7.5 km distant) during the period 
of the Geminid shower on the nights of Dec. 10/11, 1948; Dec. 13/14, 1949; 
Dec. 12/13, 1950; and Dec. 13/14, 1950, have been analyzed. Figure 1 shows 
three such typical records, taken during the night of Dec. 10/11, 1948. All 
three of them show the doppler whistle due to the motion of the meteroid, as 
can be noticed at the left end of each record, followed by the body doppler 
beat wave. A graph is drawn, in each case, to show the variation of fa with 
time. These graphs, in general, show periodic fluctuations, as shown in paper 
1. However, a mean line drawn through the observed points shows on the 
average a linear variation of fy with time. The gradient of this variation may 
be either positive or negative, and sometimes a single observation may show 
both positive and negative gradients, one followed by the other. Four different 
typical cases of these observed df,/dt gradients are illustrated in Figs. 2 to 5. 
The interpretation of these different cases will be given shortly. It has been 
shown in paper 1 that 
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Fic. 1. Three typical meteoric body doppler records taken during a Geminid shower 
period. Each record is little over a second in length, the distance between two bright spots on a 
single line being 1 second. 
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Fic. 2. Example of an fy-t curve showing a positive gradient dfa/dt. 
Fic. 3. Example of an fy-t curve showing a negative gradient dfz/dt. 
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Fic. 4. Example of an fa-t curve showing an initial negative gradient followed by a positive 


gradient dfa/dt. : ; Wie bs ; 5 
Fic. 5. Example of an fa-t curve showing an initial positive gradient followed by a negative 


gradient df,/dt over a much longer period. 
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(1) dV,/dt = (A/2) (dfa/dt) 
and 
(2) V,? = Ro(dV,/dt) 


where V, is the radial component of the velocity V, and Rp is the range of the 
point of normal incidence on the meteor path (79 say) as observed at South 
Gloucester. As an approximation, the value of Ro has been taken to be the 
minimum range read from the pulsed radar range-time records, recorded at 
the Metcalfe Road field station. The wave-length \ in this case is 9.994 m. 
The values of V, have been calculated for a total number of 90 observations and 
presented in Fig. 6, plotted against the corresponding values of total doppler 
echo durations. It is seen from this figure that, although the points are 
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Fic. 6. Plot of the experimental values of Vs vs. total doppler echo durations. 


randomly distributed, there is a tendency for higher velocities to occur in the 
case of smaller echo durations, and vice versa. If h, is the elevation of the 
meteor radiant, then the vertical component of V, is given by 


(3) Vo = Vssimks. 


This equation has been used to calculate the velocities, V,, of the above 90 
observations, assuming the elevation of the Geminid radiant in each case. 
This assumption may not be quite correct as some of the observed meteors 
may be non-shower. However, it is expected that at least half the number of 
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cases, if not more, are Geminid meteors. These values of V, have then been 
plotted against the corresponding total doppler echo durations. This graph 
is plotted in Fig. 11 along with the theoretical V,—7,, curves to be discussed 
later. 

The sign of the experimental values of V, may either be positive or negative, 
as can easily be seen from equation (2). As the actual sign of V, or V, cannot 
be determined for individual cases, the absolute values only have been plotted, 
neglecting the sign. It has already been mentioned that the experimental 
gradients df,/dt take both positive and negative values. The sign has been 
neglected when substituting in equation (1), as negative values would give 
imaginary values for V,, which are not possible. The significance of the sign of 
df,/dt can be explained as follows. The positive gradients (or increasing fre- 
quencies) indicate movement of the effective point of reflection away from 7, 
while the negative gradients indicate the opposite case. Considering the 
meteor train on one side of the maximum echo duration level, let m be the 
point of maximum echo duration, e the end point, and a the average or effective 
point of reflection approximately midway between e and m. The sign of the 
gradient df,/dt would then depend upon the relative positions of a and 7p. 
The positions of m and ro are fixed for a particular meteor, and depend to a 
large extent on the elevation and azimuth of the radiant of the meteor, and 
also on the azimuth of the meteor position in the M region with respect to the 
observing station. Let us now try to examine the gradients observed in Figs. 
2 to 5. 

Figure 2 is a very good example showing a linear increase of the average f, 
with time. This indicates that the effective point of reflection a is steadily 
receding from ro. This happens when 7p is either at a; (the initial position of 
a) or to one side of a; towards e. Figure 3 shows quite the opposite phenomenon, 
which indicates that the effective point a is steadily approaching ro. This 
happens when 7 is either at m or to one side of m away from e. In both the 
above cases, it is suggested that the effective point of reflection is drifting 
downwards, and that the effect of the drifting of the lower effective point of 
reflection is negligible, as can easily be seen later in this paper. In Fig. 4 is 
shown a case of an initial decrease of the average f, for about a second, followed 
by a steady increase for another one and a half seconds. This indicates that 
the effective point first moves towards ro, crosses it, and then steadily moves: 
away from it. This happens when rp lies between m and a;, closer to a; in this 
case. Here also, the drift is suggested as downwards. Figure 5 shows an interest- 
ing case, where there is an initial increase of average fz for about a second, 
followed by a steady decrease over a much longer period. This is a specially 
long duration meteor persisting for about 14 seconds. It is quite likely that ro 
in this case lies at or below the initial position of the lower effective point of 
reflection on the meteor train. Then the initial rise in f, would be due to the 
drift of the lower effective point of reflection away from ro and towards m. 
The effect of diffusion of the lower portion of the train has predominated 
initially as that is nearer to ro. Later on, the drift of the upper effective point 
of reflection gradually shows up. Its drift is now towards ro and hence the 
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average f, decreases with time. It should be pointed out here that the negative 
gradient, df,/dt, becomes less steep with time. This may be explained from 
the fact that the diffusion coefficient, and hence the velocity of drift of a, 
steadily decreases with decreasing height, as will be shown later in this paper. 
This non-linearity in the variation of fy with time cannot be detected in the 
other cases, as the total durations are much smaller. However, it is interesting 
to note that the phenomenon of a steady drift of the effective point of re- 
flection is clearly shown, even in the case of durations less than 2 seconds (Fig. 
3). Thus the smoothed variations of f, with time, evident in these four typical 
illustrations, show indirectly that, to produce this observed effect of drift, 
several points of reflection must exist along the meteor train even at 1 second 
after its formation. 


IONIC DENSITIES AND THEIR RATE OF CHANGE WITH HEIGHT ALONG THE 
METEOR TRAIN 
Before one actually attempts to find the echo durations from the various 
segments of the meteor train, it is important to know how the ionic densities 
vary along its length. Herlofson (1948) has derived an expression for the re- 
lation between linear ionic density V per cm and height // in the upper at- 
mosphere, which is given by 


ss git Bs Ga = 
) ee att ee 
where p = V/Ninax, -Vmax iS the maximum ionization density present in the 


train, Hy is the height of .V,4,, and H is the scale height. Taking the values 
of H for the VW region as determined by the Rocket Panel (1952) the variation 
is found to be from 6.16 km at 80 km height to 7.26 km at 100 km height. 
Thus the average value of H for this region may be taken as 6.5 km. Using 
this value, a graph is drawn to show the variation of p with (47—H) and is 
shown in Fig. 7. A similar graph has already been drawn by Clegg (1952), but 
for H = 10 km. 

In order to calculate the instantaneous values of the rate of vertical shift 
of the end point of reflection, say dH,/dT, the values of 
i dN , dp 
(2) = Nuax 

dil dil 

are needed. By differentiating equation (4) with respect to Hf (/7) and H 
being constants), the following expression can be derived 


a eS) bo Py Je Se) 
) dH iH yexp H ( yi 3 exp rT ( yexp i l¢ : 


\ graph is now drawn to show the variation of dp/dH with (J/— 9) for the 
same average value of H = 6.5 km, and is reproduced in Fig. 8. Let us now 


proceed to see how the echo durations from the various positions on the 


meteor train may be calculated. 
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RADAR ECHO DURATIONS FROM VARIOUS SPECIFIC POINTS ON THE 
METEOR TRAIN 

Considering a specific point along the meteor train, defined by a particular 
value of (H—H)), its value of p can be read off directly from the graph of 
Fig. 7. If D is the value of diffusion coefficient at that level, as the train 
expands radially, the volume density n, at a radial distance r from the axial 
concentration, at instant ¢ after the trail formation is given by Carslaw and 
Jaegar (1947) as 
(7) n = (N/4rDt) exp(—r’/4Dt). 
Huxley (1952) has assumed that, in the case of long duration meteors, scatter- 
ing occurs, not throughout the expanding cloud of ions, but at that surface 


where the concentration at that instant has attained a critical concentration 
n,. per cc given by 


(8) n,. = 1.24X 104f? 
where f is the operating frequency in Mc/s. Using the same plausible assump- 


tion in the present case, where the frequency of transmission of the doppler 
radar at Ottawa is 30.02 Mc/s, we get 
n, = 1.1175X10' electrons per cc. 

Replacing , for n in equation (7), we have a relation between r and ¢ giving 
the time taken by the train to expand to a cylinder, of radius r, from the 
surface of which reflection takes place. As the train expands, r reaches a 
maximum value (say 7max) in time 7 and then decreases to zero in time 7 
from the start. Huxley (1952) has shown that 


(9) T, = 2.77, 

and assumed that the echo persists for a total duration of 
(10) T = 27;. 

Using the condition dr/dt = 0 for rmax, we get 

(11) In 7, = In V—In 4xDn.—-1, 

and using the condition r = 0, we get 

(12) Ts = N/(4aDn.). 

Using equations (9) and (10) we have 

(13) ar V/(5.49Dn.). 


The total echo duration 7 from a point can also be calculated trom equations 
(10) and (11), and still get almost the same values. But as equation (13) is a 
simpler form it will be used in the present paper for calculation of echo dura- 
tions at different points in the meteor train. Thus if D is constant throughout 
the M region, the echo durations 7° are directly proportional to the value 
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pb = N/Nwmax of the point in the train under consideration. But, from measure- 
ments of durations based on low density trains (V < 10” electrons/cm), 
Greenhow and Neufeld (1955), obtained a linear variation of log D with height 
above ground in the 80- to 110-km region. Using their experimental graph, 
the following empirical relation between D (in cm*/sec) and H (in km) may 
be adopted. 

(14) log D = 0.06H—8. 


Substituting this in equation (10) and using the proper value of ”,., we now 
obtain 
(15) T = $.2759x i+ 
where 6 = log Nimax- This is a very simple equation which gives directly the 
echo durations 7 for different values of H, for a particular meteor of known 
Nmax and Ho, the values of p against H—H being obtained from Fig. 7. 

It is of interest to determine the height, H,,, at which maximum duration 
occurs for any meteor, that is where JT is maximum (say, 7;,). For the same 
meteor, } is constant. Differentiating equation (15) with respect to H, we have 


(16) | an as 0.138 p). 
Now, 
dT/dH = 0 
when either 
T=0 (T aminimum) 
or 
dp/dH = 0.138p (T amaximum). 


Thus the condition for T,, is given by 


(17) SO ns ttl 
p 
Now, from equations (4) and (6), we can write this as 
(18) A/H = 0.138 
where 


1 {thaoltg)H alt) 4 


Using H,, for H in this equation, we finally get the condition for T,,, as 


(19) Hyo—Hm = 2.46 km. 


The same result may also be derived, starting from equation (11) instead of 
from (13). It is interesting to note that maximum echo durations in a meteor 
train always occur at a height H,, which is 2.46 km less than Hy, that of 
maximum ionization. This is true as long as the linear relation between log D 
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and H as given by equation (14) holds good. It remains to be seen if this 
result can be verified from observational or other evidence. 

The effect of recombination of electrons with positive ions is neglected in 
this investigation, compared with the effect of diffusion of ions, as the value 
of recombination coefficient for the meteor trains is pointed out by Greenhow 
(1952) to be as low as 10~-” cm?/sec. 

Equation (15), which is useful for the calculation of echo durations 7 from 
different segments of a meteor train, may be expressed in a simpler form as 


(20) log T = log p—0.06H+b—7.2778. 


Expressions for 7,, occurring at H,, and To occurring at Ho, may also be 
written as 


(21) log 7) = —0.06 H)+b—7.2778 
and 
(22) log Tm = log pn—0.06 H,, +6 —7.2778. 


However, as it has already been shown that H,, = (41)—2.46) km, the value 
of pm corresponding to (77,,—Ho) = —2.46 km may be read off from Fig. 7. 
Thus we have p,, = 0.865 and hence the expression for 7,, may now be written 
as 
(23) log Tn = —0.06 H,+6—7.3408, 

= —0.06 Ho+b—7.1932. 


From equations (20) and (23) we obtain 
(24) log (T/Tm) = log p—0.06(H — Ho) — 0.0846. 


This is an interesting relation showing that the ratio of (7/Tm) = 4, say, 
is independent of b, but is dependent only on (H7—F1p), as p is also dependent 
on (H—H,). These ratios of 7/7, are calculated for different values of 
(H—Hp) and presented in Table I. It can be seen from this table that the 
effective vertical thickness of the M region is nearly 22 km, which is in reason- 
able agreement with previous investigations. It can also be seen that, if e; 
and e, are the upper and lower ends of reflection from a meteor train, then the 








TABLE I 
VaLuEs OF T/Tm, dp/dH, AND V-Tm FOR DIFFERENT VALUES OF (H—Hpo) 

H-H), VeTm, 

km p q=T/Tn dp/dH km 
—7.25 0.010 0.0224 0.01 51.81 
—7.0 0.030 0.0649 0.06 8.28 
—6.0 0.150 0.2828 0.21 2.80 
—5.0 0.380 0.6240 0.24 3.25 
—3.2 0.745 0.9541 0.175 10.81 
—2.46 0.865 1.0 0.119 +o 

0 1.0 0.8230 0 —8.80 

5.0 0.745 0.3073 —0.075 —13.63 

7.66 0.560 0.1600 —0.065 — 24.60 

10.7 0.380 0.0713 —0.05 —52.03 


15.0 0.210 0.0218 —0.03 — 163.30 


i 
¥ 
| 
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Fic. 7. Theoretical curve of p vs. (H—Hpo) for an average M region in which H = 6.5 km. 


Fic. 8. Theoretical curve of dp/dH vs. (H—Ho) for an average M region in which H = 6.5 
km. 


height difference between e; and m is nearly 17.5 km and that between m and 
€2 is only 4.8 km. Thus the portion of the train below the maximum echo dura- 
tion level is considerably smaller in length compared to that above it. The 
determinations of the actual values of 7,, for different combinations of Ho 
and 6 can be made using equation (23). Some specific cases will be considered 
later in this paper. 


VERTICAL COMPONENT VELOCITIES OF THE DRIFT OF THE END POINTS 
TOWARDS THE LEVEL OF H,, 


It has been shown in Table I that the durations of radio response from 
various parts of the meteor train are low at the two extremes and increase 
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to a maximum at a height H/,, given by equation (19). Thus as time advances 
after the initiation of the ionization, the end point of radio reflection gradually 
moves towards the maximum duration level. If d7 is the difference in dura- 
tions at two levels differing in height by a small value dH at any height H, 
then the instantaneous rate of vertical drift of this point, when the end of the 
train coincides with that height, is given by dH/d7. These values may be 
calculated, in general, using the relations 


a dH _ dp / ap 
a ar” af dH’ 
(26) dp _ 54nDn, 


an Fes 


and obtaining the values of dp/dH from Fig. 8. This procedure has to be 
adopted when an average value of D is taken for the whole M region. This 
involves the assumption that 7 = 27), but the actual limits between which 
T can vary are given by 7; and 7». The instantaneous values of dH/dT in 
each of these three cases, 7, 71, and 7, may be derived and denoted as V, 
Vi, and V2, respectively. They have been calculated for an average value of 
D = 2.5X10' cm/sec, which is the value at a height of nearly 90 km. These 
values are presented graphically in Fig. 9 as continuous curves. 





dH/dT (KM/SEC)—> 


os 











H—-H, (KM) —> 


Fic. 9. Theoretical curves showing the instantaneous velocities V, Vi, V2, and V. of the 
drift of the end point of radio reflection at different positions along a meteor train, defined by 
(H—FH)). 
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As d7/dH is zero when T is maximum, the velocities are infinite at the 
level where 7 is maximum. From Fig. 9, it can be seen that for Nmax = 10", 
the magnitudes of these velocities are very high, ranging from 1 km/sec 
upwards. For lower densities, the velocities are still higher, but for higher 
densities the reverse is true. For example, when Ninax = 10", the velocities 
are reduced by a factor of 10 from those for Vmax = 10®. 

Now, it is of interest to see how the values of V are altered if the variation 
of D with height is considered. The corrected values of V, = dH/dT can now 
be calculated using the reciprocal of d7/dH in equation (16). The corre- 
sponding values of p and dp/dH for different values of (H7—H) are given in 
Table I. The value of 7 in this equation may be written as g7, so that equa- 
tion (16) may now be written as 


(27) V.Tm = ——_——__-.. 


dp _ ) 
o( 0.138 p 


Using this equation, the values of V.7;, have been calculated and shown in 
Table I. The values of g have also been taken from Table I. From equation 
(23), it can be seen that 7;,, is dependent on the value of Ho. Hence if the 
variation of-D with H is assumed, the value of V,. becomes dependent on 
H). As the meteors generally occur in the 80-110 km region, we may calculate 
the values of V, for a special case, .06 H,, = 5, so that 
T,, = 5.275), 10 =, 

= 4.563 X 10°-%, 
where p, = .865 as has already been seen. Using this value of 7;,, the values 
of V,, in this special case, have been calculated and shown in Fig. 9 as dashed 
curves. It is interesting to see that these curves are asymptotic to the Y-co- 
ordinate at (H7—H,)) = —2.46 km. That is, the velocities become infinite at 
the position of 7;,. In the previous case when D is assumed to be constant in 
the height, 7,, appeared at (7—H)) = 0 itself. That is why the continuous 
curves are asymptotic with the Y-axis itself. 

It remains to be seen how these high instantaneous velocities can be verified 
practically. At this stage, these curves are of only theoretical importance. 
Knowing the values of 7,, at height //,,, if the end point of the meteor train 
giving radar response is at height H,, then the average velocity with which 
the end point e moves towards the point m of maximum T is given by 


H,,—H 
97 = ett. 
(27 ) V. ; =e . 


where 7°, is the duration for which the echo persists at height /7,. This velocity 
is negative if the point moves downwards and vice versa. Now let us assume 
from Table I that the points of radio reflection occur at approximately the 
positions given by (1—H)) = —7.25 and 15.0 km respectively, where the 
values of 7/7, are nearly 0.02. Using the suffixes 1 and 2 to represent the 
upper and lower end points, we obtain 


H,—H., = —17.46, H,—H., = 4.79, 
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and using the values of 7/7, corresponding to the corresponding values of 
(H—H)) from Table I, we obtain 


Tm—T 2, = 0.9782 Tm, Tm—Te, = 0.9776 Tm. 


Hence we may write down the expressions for the velocities V, of the two end 
points of radio reflection, using equation (27) as 





: —17.85 4.90 
(29) Ve, = —“- ; Ve. = T,. 


As the upper end point moves down towards the maximum echo duration level 
at H,,, it is clear that V,, is negative. Similarly, as the lower end point moves 
upwards, V ,, is positive. It can be easily inferred from these equations that the 
upper end point of reflection moves with faster average vertical component 
velocities than the lower one. These velocities are calculated for different 
values of 7,, and are represented graphically in Fig. 10. 





Ve (KM/ SEC) —» 





Fic. 10. Theoretical curves of the average velocities, V., of the drift of the end point 
of radio reflection, vs. Tn, the total echo durations. Suffixes 1 and 2 stand for the upper and 
lower end points of radio reflection respectively. 


Let us now proceed to see how the theory, outlined above, can be used to 
explain the observed values of V, the vertical components of the velocities 
of drift of the effective point of radio reflection along a meteor train. 


THEORETICAL EXPLANATION OF THE OBSERVED VELOCITIES OF DRIFT 
OF THE EFFECTIVE POINT OF RADIO REFLECTION 

It has already been mentioned that the effective point of reflection may be 

assumed to coincide approximately with the mid-point of the effective length 
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of the train to one side of the maximum echo duration level. Hence the velocity 
of drift of the effective point of radio reflection may be taken to be half that of 
the end point of radio reflection from a meteor train. In that case, half the 
velocities represented in Fig. 10 would correspond to the velocities V,, referred 
to in equation (3), in the ideal case. However, in practice, underdense trains 
may be shorter in length than the overdense trains and hence may occupy 
a smaller thickness of the M/ region. Hence, velocities V, = V,/2, have to be 
recalculated for specific practical cases using different smaller values of 
(H,—H,) than those used, in equation (29), for the ideal case. 

For this purpose, the average height of occurrence of the meteors observed, 
in this investigation, should be known at least roughly. It has already been 
suggested, following paper 1, that at least half the total number of meteors 
observed may belong to the Geminid shower. Millman (1954) has shown that 
the experimental values of mean heights of meteoric ionization show a linear 
variation with the velocities of the meteors, within a particular range. The 
average velocity of Geminid meteors has been shown to be 35.3 km/sec (Mill- 
man 1954; McKinley 1951). Using this value of velocity, the corresponding 
mean height of a Geminid meteor may be read off the straight line graph 
(Millman 1954). This height is found to be 92.6 km. 

In fact, values of (/7,,— 7.) for the actual meteors are expected to be much 
smaller than those derived from Table I. In such a case, these values for the 
meteor train below the maximum echo duration level, H,, may be negligibly 
small compared with those for the upper portion of the train. Thus the effect 
of the drift of the lower effective point of reflection may be completely 
neglected, as has already been pointed out in the experimental observations. 
Now, assuming the initial position of the upper effective point of reflection 
to be situated at a height given by 


(30) H, = 92.6 km 


for most of the observed meteors, the height of the end point of reflection 
may also be calculated 


(31) H, = 2Hy—Hm. 


The position of the maximum ionization level, Ho, may occur at different 
heights for different meteors, and hence (H,—H,) may vary for different 
values of Zo. Three different cases of Hyp = 85, 90, 92.6 km may be considered 
which correspond to the effective train lengths of (H,—H,) = approximately 
20, 10, 5 km respectively. The upper limit of 20 km is higher than that given 
in Table I, which is only 17.5 km, and hence would include even the very 
bright meteors of longer duration. The lower limit of 5 km can be taken as 
corresponding to the underdense meteors of very small durations, of the 
order of 1 second or so, and the value of 10 km toa statistically average meteor. 
We may also consider one more special case where 1) = 96 km and H, = 92.6 
km corresponds to the initial position of the lower effective point of reflection. 
This case probably does not occur frequently. However, its consideration 
may give us a chance of verifying the velocities of drift of the lower effective 
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point of reflection, at least in those cases when they are prominent. These 
four cases are represented in Table II, and denoted by numbers 1 to 4. The 


TABLE II 
THE VARIOUS PARAMETERS FOR FOUR SPECIAL CASES OF Ho 








pe 











No. Ho H,—H, Hm — Ha H.—Ho 

1 85 —20.12 —10.06 17.66 145 
2 90 —10.12 —5.06 7.66 .560 
3 92.6 —4.92 —2.46 2.46 .920 
4 96 1.98 0.94 —4.34 .540 





average vertical component velocities of the drift of the effective points of 
reflection in the above cases may now be calculated from the formula 
— ane ie ta 

(32) V, = a = ee = 

& m e m ée m é 
The values of (H,,—H,) may be taken from Table II. The values of 7, have 
been calculated using equation (23) in the above four cases for different values 
of 6 ranging between 12 and 14.5. The corresponding values of 7, may be 
calculated using equation (24) substituting 7°, for 7 and H, for H. The various 
calculated values of 7,,, 7, and V, for the different cases are then shown in 
Table III. The values of V, in this table are the averages for the total echo 
durations 7, in each case. 

The values of V,, thus obtained, are represented graphically in the form of 
curves of V, vs. 7, in Fig. 11. The sign of V, is neglected in drawing these 
curves. The experimental values of V, have also been plotted as dots in the 
same figure, for purposes of comparison, against their corresponding total 
doppler echo durations. The curves numbered 1 to 4 refer to the four different 
cases of Table II. It has already been mentioned that case I is an extreme 
case for very long duration meteors and hence curve 1 should be better fitting 
by points at the longer durations. From a comparison of these curves with 
the experimental values, it can be seen that though there are a few experi- 
mental points lying in the region of the theoretical curves, in general the 
experimental values are much smaller in magnitude than those predicted by 
theory. 

It will be interesting to consider the various possibilities for this discrepancy 
between theory and experiment. One important fact is that the theoretical 
values represent drifts of either the upper or the lower effective point of 
reflection, whereas the experimental values represent the average drifts of 
the effective point of reflection for the whole train. In comparing the two, the 
length of the train below the maximum echo duration level has been neglected 
compared to that above. In other words, the drift of the lower effective point 
has been neglected compared to that of the upper one. In fact, as they both 
move in opposite directions, the total effective point moves with smaller 
velocity than that of either of the two effective points and is actually the 
difference between the two. This partly explains the smaller experimental 
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Fic. 11. Theoretical curves of the average velocities, V,, of the drift of the effective points 
of radio reflection, vs. Tn, superposed on a plot of the experimental points of Vy vs. Tn. 
Numbers | to 4 represent the four cases given in Table IT. 


values. This is especially true in the case of faint meteors, for which the two 
portions of the train on either side of the maximum echo duration level may 
be nearly equal in length, as Greenhow and Neufeld have pointed out (1957). 
Hence this may partly explain the wide discrepancy between the theoretical 
curves and experimental points in Fig. 11, in the case of durations less than 
about 3 seconds. 

There is yet another possibility which needs consideration. In the case of 
overdense trains, the assumption that the length of the train below the maxi- 
mum echo duration level is negligible compared to that above may be reason- 
able. In such cases the observed drift is mainly that of the upper effective 
point of reflection. From a close examination of the values of g (= T/T») vs. 
(41—Hp) in Table I, it will be seen that the 7,, vs. H,, curve is somewhat 
exponential in shape. But, the average velocities V,, are calculated from 
equation (28), which involves the assumption that the above shape is linear. 
The effect of the exponential shape is to make the true average velocity less 
than that calculated from equation (28). This may partly explain the lower 
experimental values, compared with the theoretical curves in Fig. 11, in the 
case of echo durations more than about 3 seconds. 

Having thus explained the discrepancy between the theory developed and the 
experiment performed to a certain extent, it appears that the present theory 
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represents, within a reasonable approach, the diffusion of ionization in a meteor 
train. It will also be interesting to try to explain the present experimental 
results from a theory developed purely from previous experimental obser- 
vations of other workers. Some of the discrepancy discussed above may be 
attributed to the atmospheric turbulence (which need not be small scale) 
in the M region. In order to estimate quantitatively the effect of turbulence 
the following procedure is followed. 

Greenhow (1952) has pointed out that the influence of atmospheric tur- 
bulence on the meteor train makes the actual durations of the radar echo 
much less than those predicted by theory. Physically, this has been explained 
as due to the breaking up of the train by turbulence. He has also pointed out 
that this effect is larger for longer wave-lengths. Since the wave-length used 
in this investigation, which is 9.994 m, is greater than either 8.4 m or 4.2 m 
used in his experiment, the turbulence effect is expected to be more pronounced 
in the present investigation. Thus the values of 7,, have to be corrected for 
turbulence. These corrected values of 7, must be estimated from a knowledge 
of the previous experimental data. 

Experimental relations between the visual magnitudes of meteor trails, 
the maximum ionizations produced in them, and the total echo duration by 
radar, have been investigated previously by various workers. The visual 
magnitude scale is defined so that a change of one magnitude represents a 
factor of 2.5 in luminosity and is referred to the point of maximum luminosity, 
which will coincide with the maximum ionization. Hence, we have 


(33) M = —2.5 log Nmaxt+C 


where M is the visual (zenithal) magnitude, and C is a constant. Kaiser 
(1953), from a comparison of visual and radio echo rates, obtained a value of 
C = 35 for Perseids with low density trails. He has also given a theoretical 
expression for the above relation as 


(34) M = 25.35—2.5 log (Kuv® Nmax/8). 


Using v = 60 km/sec for Perseid meteors and taking p = 50 a.m.u., he 
obtained K/8 = 7.5X10-*. Using the same constants, but v = 35.3 km/sec 
for Geminid meteors (Millman 1954; McKinley 1951), we get C = 36.75. 
Thus for Geminid meteors in the present case, we can rewrite equation (32) as 


(35) M = —2.5 6+36.75. 


From a study of 256 Geminid meteors, Millman and McKinley (1956) have 
obtained an experimental relation between the visual magnitude and maximum 
echo durations as 


(36) log Tm = 1.15 —0.349 M. 


Using equations (34) and (35), for each value of Nimax, the corrected values 
of Tm (= Tm’ say) can be calculated. These values of 7,,’ are listed in Table 
IV. As the values 7,’ in Table IV are much smaller than those of 7,, in Table 
III, it is reasonable to assume that the corresponding values of 7,’ = 0 in all 
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TABLE IV 
VALUES OF JT,’ AND V,’ CORRECTED FOR TURBULENCE 


V,’, km/sec 








Tn’, af 
b sec 1 2 3 4 

13 .4613 —21.81 —10.97 —5.33 +2.04 
13.25 .7638 —13.17 —6.63 —3.22 +1.23 
13.5 1.262 —7.97 —4.01 —1.95 +0.75 
13.75 2.084 —4.83 —2.43 —1.18 +0.45 
14 3.443 —2.92 —1.47 —0.71 +0.27 
14.5 9.397 —1.07 —0.54 —0.26 +0.10 
14.75 12.83 —0.78 —0.39 —0.19 +0.07 


15 25.70 —0.39 —0.20 —0.10 +0.04 


the four cases. The corrected values of V,(= V,’, say) may now be calculated 
using the new formula 


(37) Vi.’ = (Ao—95.06)/Tn’. 

These values are also listed in Table IV. Now, using these, the corrected 
curves of V,’ vs. 7,’ have been plotted in Fig. 12, together with the experi- 
mental values shown previously in Fig. 11. 


From a comparison of the curves of Figs. 11 and 12, it can be seen that the 
effect of turbulence is to shift the theoretical curves towards the origin. 
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Fic. 12. Theoretical curves of V,’ vs. Tm’, superposed on the experimental points plotted 
in Fig. 11. 
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Two things are taking place here. The practical values of 7, are lower for 
each value of .,, and the corresponding V, values increase because of the 
reduced time duration of drift of the effective point of reflection. However, 
it is seen that the shift due to the correction for turbulence is least for curve 
1 and increases progressively in the case of curves 2, 3, and 4. It may be 
noted here that the theory developed at the beginning of this paper for 
getting echo durations from different segments of the train assumes that 
scattering occurs at the surface of critical volume density ”,. This is true 
only in the case of long duration meteors. Hence it is to be expected that the 
theory will be sound for longest echo durations and that errors may increase 
as the echo durations decrease. In other words, this may explain why the 
turbulence corrections for actual cases are a minimum for curve 1 (which is 
supposed to apply to long duration meteors) and increase for curves 2 and 
3 (which are supposed to be cases of decreasing order of echo durations). 

It can now be seen from Fig. 12 that a great majority of the experimental 
points lie in the region between the curves 1 and 4. In particular, curves 2 and 
3 can be said to explain very well the observations of meteors with durations 
above 3 seconds. It is interesting to see that the majority of the observed 
meteors lie between curve 3 and the X and Y axes. This is explained for the 
shorter duration echoes by the fact that the rectangular co-ordinate axes 
represent, theoretically, the limit to which the curves tend as the echo dura- 
tions approach zero. While the experimental points below curve 3 for the 
longer duration echoes can be explained by curve 4, representing velocities 
of drift of the lower effective point of reflection, the good agreement between 
curve 4 and most of the experimental points at lower durations may be 
accidental. 

The theory outlined above, which is based on practical results of previous 
workers, seems to explain reasonably well the observed velocities of drift 
of the effective points of radio reflection. This shows that the true echo dura- 
tions from a meteor train are much smaller than those predicted by the 
original theory based on fundamental assumptions, where no turbulence 
effect is taken into account. The two possibilities discussed before to explain 
the discrepancy between the original theory and experiment will still have a 
significant effect on the new theoretical curves. The effect is again to reduce 
the theoretical velocities slightly, which only gives a better agreement between 
the theoretical curves and the experimental points in Fig. 12. It must be 
remembered that the experimental values of V, have been obtained, on the 
assumption of the elevation of the Geminid radiant in each case. As about 
half of the observed meteors were probably sporadic, this may explain the 
disagreement in ome of the experimental points. Another point is that the 
theoretical curves, corrected for turbulence, have to be modified in the case 
of sporadic meteors, by using the proper velocity for obtaining the constant 
C in equation (33). In fact, this constant C may be true only for short duration 
meteors, as Kaiser (1953) has derived the various corresponding constants 
of equation (34), using short duration meteors only. Hence the constant C 
is probably different for longer duration meteors, which may explain the fact 
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that there are very few points near curve 1 corresponding to the long duration 
meteors. The agreement between the latter theory, based on previous experi- 
mental results, and the present experiment may be considered as reasonably 
good within the expected limitations. Hence it may be concluded that the 
atmospheric turbulence in the M region has a significant effect on the duration 
of meteor echoes. 


DISCUSSION 


It will be interesting at this stage to discuss the various physical aspects of 
the present problem. In developing the theory, three basic assumptions have 
been made. The initial distribution of ionization along a meteor train has 
been assumed to follow the equation developed by Herlofson (1948). This 
ionization is then assumed to diffuse with time according to the theory pre- 
sented by Kaiser (1953). The actual values of the ionic diffusion coefficients 
at various altitudes have been adopted from the experimental results of 
Greenhow and Neufeld (1955). This assumption involves the suggestion that 
diffusion of ionization is not affected by small scale turbulence, which is 
contrary to that of Booker and Cohen (1956). Further, the mechanism of 
back scattering has been assumed, after Huxley (1952), as that taking place 
at the surface of the ionization cylinder where the volume concentration has 
attained a critical value ”, of equation (8). Furthermore, the meteor train, 
which is straight initially, is known to distort sinuously soon after its formation 
(Rao 1958) owing to large scale wind shears. This has been assumed to be the 
main reason for the formation of several points of reflection along the train. 
It has been seen, however, that atmospheric turbulence, of some scale (which 
need not be small scale of the order of 1 m), breaks up the train into smaller 
segments. It may be considered that these tend to round off, giving blob-like 
structures instead of cylindrical. This involves a certain amount of roughness 
in the train, but not to the extent assumed by Booker and Cohen. Going to 
the other extreme, Manning (1957) assumes that meteor reflection is usually 
completely specular as from a perfectly smooth mirror-like train. Manning’s 
theory involves the successive formation of ‘‘glints’’ along the meteor train 
as it is sinuously distorted by large scale turbulent eddies or wind shears. This 
means that the effective length of the train increases initially with time until 
the maximum length is established, and then it decreases owing to the diffusion 
of ionization. The assumptions of the present theory are essentially the same 
as that of Manning except for the slight roughness involved. It has not been 
possible to detect the initial increase in the train length during the present 
experiment. In fact, the present experiment is mainly concerned with the 
diffusion of ionization after the total length of the train has been established. 
It will not be out of place to mention here that several b type echoes (McKinley 
and Millman 1949) are usually observed in the back-scatter pulsed radar 
experiment during meteor shower periods, especially the Perseids. It is seen 
that most of these echoes exhibit properties of rough trains almost instan- 
taneously after the formation of the head echo. One of the authors (M.S.R.) 
has undertaken analysis of such echoes for determining how the ionization 
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diffuses with time. Preliminary results of this analysis seem to support the 
present theory. Full details of this work will be published soon. Finally it may 
be said that Booker and Cohen’s assumption, that small scale turbulence 
makes the meteor train rough almost instantaneously after its formation, 
may not be correct. On the other hand, Manning’s smooth-train theory also 
does not explain many observed properties of meteor echoes. However, it is 
believed that large scale turbulence of the order of 100 m to 6 km is mainly 
responsible for controlling the distortions in the train (Rao 1958). The higher 
scales may produce sinuous distortions and the lower ones tend to break up the 
individual segments and produce blob-like structures distributed all through 
the train. Further in the case of certain meteors, successive fragmentation 
may also produce sufficiently rough trains to give reflections from the entire 
train length, almost immediately after the formation of the train (Greenhow 
and Neufeld 1957). Such irregular distributions of ionization might be re- 
sponsible for the great scatter of experimental points in the case of meteor 
echo durations less than about 3 seconds. 


CONCLUSIONS 


Experimental evidence has been provided at the beginning of this paper to 
substantiate the postulate of the drift of the effective point of reflection along 
a meteor train towards the maximum echo duration level. From a large number 
of observations it has been found that larger velocities of this drift invariably 
occur in the case of underdense trains of total echo durations less than three 
seconds, and that in the case of overdense trains of longer echo durations 
these drift velocities are always smaller. Assuming that at least half the ob- 
served meteors belong to the Geminid radiant, the experimental values of 
vertical component velocities of this drift have been found to be in reasonably 
good agreement with the theoretical curves, if the effect of turbulence on the 
echo durations is also taken into account. On the other hand, this agreement 
may be taken to provide sufficient evidence for the presence of turbulence on a 
scale of the order of 100 m to 6 km or even less in the meteoric M region of 
the upper atmosphere. 

Theoretical curves for the instantaneous vertical component velocities of 
drift of the effective point of reflection at different positions along a meteor 
train have also been presented in this paper. From theoretical considerations, 
involving an empirical linear relation between log D and H, it is predicted 
that maximum echo durations from a meteor train always occur at about 
2.46 km below the height of maximum ionization. It remains to be seen if 
these two theoretical conclusions are verified in the future by observation. 
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THE RESOLUTION CORRECTION IN THE SCINTILLATION 
SPECTROMETRY OF CONTINUOUS X RAYS! 


W. R. Dixon AND J. H. AITKEN 


ABSTRACT 


The problem of making resolution corrections in the scintillation spectrometry 
of continuous X rays is discussed. Analytical solutions are given to the integral 
equation which describes the effect of the statistical spread in pulse height. The 
practical necessity of making some kind of numerical analysis is pointed out. 
Difficulties with numerical methods arise from the fact that the observed pulse- 
height distribution cannot be defined precisely. As a result it is possible in 
practice only to find smooth “solutions”. Additional difficulties arise if the 
numerical method is based on an invalid analytical procedure. For example 
matrix inversion is of doubtful value in making the resolution correction because 
there does not appear to be an inverse kernel for the integral equation in question. 


I, INTRODUCTION 


In the scintillation method for the measurement of continuous X-ray 
spectra, the incident X rays are absorbed in a light-producing crystal, chosen 
so that the intensity of the light flashes is proportional to the energy absorbed. 
A photomultiplier converts the light flashes into voltage pulses, which are 
then amplified and analyzed as to height in a kicksorter. There is not, however, 
a unique pulse height V for a given incident energy E. This means that there isa 
spread in the values of E which produce pulses of height V; if 7(V,E) is the 
response function describing this spread, the observed pulse-height distribution 
f(V) and the true spectrum g(£) are then related through the following in- 
tegral equation: 


(1) f(V) = JT(V,E) g(E) dE. 


The contributions to 7(V,E) are of two distinct types. The first type is due 
to the failure of the crystal to absorb all of the incident energy, once an 
interaction has taken place. Thus energy can escape in the form of electrons, 
in the form of Compton-scattered photons, in the form of annihilation radiation 
for high E, and in the form of characteristic X rays for low E. The second cause 
of spreading is the statistical nature of the emission of electrons from the 
cathode surface of the photomultiplier; emission from the subsequent dynodes 
makes a smaller contribution to the spread. The effect of the statistical spread 
can be described by the integral equation 


(2) f(V) = JdE g(E) (2ekE)—? exp{ —(V—E)*/2kE}, 


if it is assumed that there is total absorption in the crystal.* 


1Manuscript received July 31, 1958. 

Contribution from Division of Applied Physics, National Research Council, Ottawa, 
Ontario. 

Issued as N.R.C. No. 4954. 

*The present analysis of the resolution correction was made in connection with some experi- 
mental measurements of the spectra from conventional X-ray machines. In this energy region 
(up to about 200 kev) it is found that the width of the photopeak at half-height is proportional 
to E}, in agreement with equation (2). In the neighborhood of 1 Mev, however, the width 
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The problem with which we shall be chiefly concerned is the solution of 
equation (2), firstly by analytical methods, and then by numerical methods. 
A partial solution to equation (2) has been given by Palmer and Laslett 
(1951). It is 


; df kEd: 
(3) g(E) = f(z) —k Fe CS 


A different partial solution given by Morton (1952) appears to be incorrect. 
The general solution will be given in Part II below. 

Equation (2) bears some similarity to the Gauss or Weierstrass trans- 
formation 


(4) f(V) = J dE g(E) (240")~"” exp { —(V—E)*/20°}. 

In equation (2) the variance is proportional to the energy E£, while in equation 
(4) it is constant. The problem of solving the integral equation (4) is ap- 
parently well known to astronomers, the first solution having been given by 
Eddington (1913): 


a“) a) We (=1)"0" df(E) 
(5) g(E) = Se Ge 


provided the series converges. Equation (5) can also be written in the symbolic 
form 


9 9 
o a 


(6) g(E) = exp (-5 2) . f(E). 


Applications of equation (4) in astronomy have been discussed for example 
by Malmquist (1922, 1941, 1943) and by van de Hulst (1946). The latter 
author has shown that the solution of equation (4) can also be given as a 
series of Hermite polynomials, when f(V) is given as a power series. In the 
mathematical literature, Hille (1926) has discussed at length the reciprocal 
relation between powers and Hermite polynomials under this transformation, 
and other authors have considered the interpretation of the inversion operator 
in equation (6), for example Pollard (1946) and Widder (1951). Occasionally 
the first two terms of equation (5), namely 
o df 
(7) g(E) = f(E)-5 oa, 


have been rediscovered independently by physicists (e.g. Owen and Primakoff 
1948, 1950), who on the whole seem to be unaware of the contributions of the 
astronomers and mathematicians to this problem. 

Recently a number of authors concerned with scintillation spectrometry 


appears to vary roughly as E}, according to Kelley et al. (1956), so that equation (2) is no 
longer strictly applicable, although it is still a good approximation. The reason for this type 
of variation is that the photomultiplier no longer makes the predominant contribution to the 
spreading, there being also the effect of a statistical spreading in the crystal itself, the nature 
of which does not seem to be clearly understood. 
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have used a method of numerical matrix inversion for the solution of equation 
(1) (Ehrlich 1955; Theus, Beach, and Faust 1955; Motz 1955; Starfelt and 
Koch 1956; Hubbell 1958), but they have not indicated any conditions of 
validity for that method. One replaces equation (1) by a finite set of algebraic 
linear equations: 


N 
(8) i= > Tae, “Bw 8 cM. 


m=1 


If the determinant of the matrix (Jim) is not zero, then there exists a unique 
inverse matrix (7ym~!) such that 


N 
(9) gm > Take G«1,2.....®. 


m=1 


The popularity of this method can be attributed to the availability of large 
digital computers to compute the inverse matrix. It will be shown in Part III 
below, however, that the method should not be used to make the resolution 
correction, although it can be used to correct for the escape of energy from the 
crystal. 


II. ANALYTICAL SOLUTIONS OF EQUATION (2) 


A number of pairs of functions which satisfy equation (2) are given in 
Table I. We shall first verify these results, and then indicate how general 
solutions can be obtained. 


TABLE I 
SOLUTIONS OF EQUATION (2) FoR V > 0 











f(V) g(E) 
1. exp sV (1—ks) exp (sE—s*kE/2) 
(ja]| + b < 1/k, where s = b + ia) 
2. cosaV (cos aE + ak sin aE) exp akE/2 
(lja| < 1/k) 
3. sinaV (sin aE — ak cos aE) exp a®kE/2 
(la| < 1/k) 
4. (V—Vo)™ (RE)™!*[hm {(E— Vo) (RE)~*} 
—m (k/E)* hmi{(E— Vo) (RE)~4}] 
ym (RE)™!2 (k/E)* Ins: {(E/k) 4} 


lo 





Consider first the following definite integral: 
(10) I= J (24kE)~” exp { —(V—E)?/2kE} dE, 
0 
where E, V, and & are real (k > 0). By changing the variable of integration to 
(V—E)E-}, it may be shown that 


I=] for V > 0, 


(11) 
I = exp (—2|V|/k) for V < 0. 


In view of equation (11) we restrict further discussion to positive V. 
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To obtain the first result in Table I, we use 
(12) g(E) = (1—ks) exp (sE—s*kE/2) 
in equation (2), with V > 0 and k > 0. Equation (2) assumes the form 


(13) f(V) =e” J “(1—ks)(24kE)-™? exp {—(V—E+skE)?/2kE} dE. 


The integral converges for any real value of s. For complex values s = 6+7a, 
the integral converges for 


(14) a? < (b—1/k)*. 
We shall further restrict ourselves by the condition 
(15) b < I/k. 

Then it can be shown that 

(16) S(V) = et”. 


Equation (16) can be called the transform of (12), and (12) the inverse trans- 


form of (16). 
The inverse transforms of cosaV and sinaV can be obtained by setting 


$s = 7a in equations (12) and (16). The results given in Table I are valid 
for the condition |a| < 1/k. 
In order to obtain the inverse transform of (V— Vo)”, we compare the ex- 


pansion 


(17) VV) — SS s"™(V—Vo)™/m! 


m=(0 
with the expansion 
(1—ks) exp {s(E— Vo) —s°kE/2} 


(18) « m 
= (IRs) oF (RE) nl (E— Vo) (REY *}, 
m=0 . 
where the h,, are Hermite polynomials. The desired inverse transform therefore 
is the result given in Table I. The explicit form of the Hermite polynomials is 
(-— 1)" m! 2 * 


(19) hin(®) = De Sat Gm — Int 


n=0 





where the sum is taken up to the largest integer, m/2 or (m—1)/2. This 
definition differs from the more usual definition: 

Hi, (2) = 2" haens2). 
If Vo = 0, some further simplification is achieved by applying a recurrence 
property of the Hermite polynomials, and the inverse transform of V™ is the 
result given in Table I. 


If we write f(V) as a Taylor series, 


(20) §(V) = f(Vo)(V— Vo)"/m!, 
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where f(Vo) is the mth derivative of f evaluated at V = Vo, then g(£) has 
the form 
(m 


9 > = 7 (Vo) ym /2 = : —1/2 
(21) g(B) = De (bE) im (E— Vo) (RE) *"} 
—m(k/E)"” hms {(E— Vo) (RE) ""}], 


provided this series converges suitably. A special case of this result is obtained 
by setting V>o = FE. Then 


2 co hE m oil m sits : “i 
(22) g(E) = > SE pe) Pr), 


where use has been made of equation (19) to evaluate h,,(0). Equation (22) 
suggests the following symbolic representation of the inversion of equation (2): 


kE d° é) ee 
9° 3(E) = exp | —— 5]. (1-—k-]. f(B). 
(23) g(E) exp ( 5 = (1 kak f(E) 


Equations (22) and (23) give the analogue of Eddington’s solution as expressed 
in equations (5) and (6). 


Il. NUMERICAL METHODS 


It has been seen that exact solutions of the integral equation (2) can be 
obtained if f(V) is given in analytical form. What is actually observed, how- 
ever, is a set of numbers 

1 VitAvV/2 
Pa ae f(V) dV dS cc Be, 
AV Jyi-avee 
where N is the number of channels. The number f; is the actual value of 
f(V) at some point within the channel usually taken to be the mid-point of the 
channel, V;. One is now forced to make some kind of numerical analysis in 
order to obtain g(£). 

Perhaps the most simple and direct method for making the resolution 
correction is to use numerical approximations for the derivatives in the 
analytical solution shown in equation (22). In order to evaluate the deriva- 
tives of f; one can set up a system of (V—1) equations obtained by writing 
each f,; (j ¥ 7) as a truncated Taylor series in the first (V—1) derivatives of 
f;. The numerical solution of these equations then gives the desired derivatives. 

The fundamental difficulty with which one is confronted in this (or any 
other) procedure is the lack of definition of the function f(V). This lack of 
definition manifests itself in large uncertainties in the high-order derivatives, 
so that the series in equation (22) must be terminated after only a few terms. 
In practice attempts to define f(V) more precisely by increasing indefinitely 
the number of channels will be defeated by the existence of experimental errors. 
However, even in the absence of experimental errors, an increased number of 
channels would lead to difficulties owing to the accumulation of round-off 
errors in the increased number of arithmetical operations used to calculate 
the derivatives. 








ee ee 


eT PIER TRF 


SURE 
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For example it may be found that only the first and second derivatives can 
be evaluated with sufficient precision to be useful. In this approximation 
one can fit the points f;-1, f;, and fi;1 with the parabola 


(24) AV) = fit (V— Vo f£O +1(V— V)?/20 fe, 
and the derivatives are given by simple difference formulae: 


fio = (fiti—fi-r)/2AV. 
fi? = Siar — 2h i +fi-r)/(AV)?. 


The truncated series of Palmer and Laslett (1951) given in equation (3) is 
then used to make the resolution correction. 

The general effect of these considerations on the resolution correction is 
the expected result that one can only sharpen peaks and valleys already 
apparent on f(V), but that one cannot reconstruct fine structure which has 
been obliterated in the measurement process. 

In the remainder of this section a discussion is given of some other numerical 
methods which have been suggested or used for making the resolution correc- 
tion. These methods are subject to the same fundamental difficulties, arising 
from the lack of definition of f(V), as those already discussed. In addition, 
however, there may be convergence difficulties which can arise for example 
because the numerical method is based on an invalid analytical procedure. 

Let us consider first the fitting of f(V) with a partial sum f, of orthonormal 
functions ¢,(V): 


(26) fi(V) = 2d ci¢i(V), 
where 

fed(V)e(V) dV = by, 
and 


ty = SF V)e(V) dV. 


If we suppose further that the system ¢; is complete, then for any square 
integrable function f(V), 

(27) lim [ (7(V) fav = 0. 

Suppose now that the inverse transform for ¢;(V) is ®;(£). Then the transform 
of 


n 


(28) Bu(E) = 2) cvbs(E) 
i= 
is f,(E). It follows therefore that if 1 is sufficiently large, the g,(E) of equation 
(28) satisfies the integral equation arbitrarily closely, in the sense of equation 
(27). This is true even if the sequence g, diverges. In this case we refer to 
members of the sequence g, as pseudosolutions to the integral equation. 
As an example of how pseudosolutions can arise, suppose that f(V) has 
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been observed with a 20-channel kicksorter, and has been fitted with a Fourier 
sine series: 


(29) f(V) = : b; sin . (0 = V $20), 


where 


1 


(30) b; = 10 "4(V) sin ae av. 


Referring to Table I, we consider for solution the sum 


(31) gn(E) = = bj)sin = —k 7 cos niE\ exp nike ° 

, ‘\ 2 ao 95 7? one 
Since f,(V) gives a better fit to f(V) than does f,_1(V), we can say that g,(£) 
satisfies the integral equation more closely than does g,-1(£). However, the 
exponential factor increases in a disastrous way with , and it can be seen 
that the sequence g,(#) might misbehave badly in spite of the good behavior 
of f,(V). 

In practice it is difficult to establish the convergence or divergence of a 
sequence g,(/) because the coefficients c; are known only numerically. It is 
therefore difficult to distinguish between good approximations to the true 
solution, and pseudosolutions. 

The trial-and-error method suggested by Lidén and Starfelt (1954) is 
subject to these same difficulties with pseudosolutions. One selects a trial 
solution g,(£), and then performs a numerical integration to see if the in- 
tegral equation is ‘‘satisfied’’. But even if one satisfies the criterion 


SUV) —fi(V) 2 dV <e 


where ¢ is an arbitrarily chosen positive number, and f,(V) is the transform 
of g,(E), there is still the possibility that g,(£) may bear very little resem- 
blance to the true solution. 

It is clear that for practical applications of the function-fitting method or the 
trial-and-error method, it is desirable to restrict the solutions by a condition 
of “‘smoothness’’. In fact Eckart (1937) has shown that there is only one 
smooth ‘‘solution’’ possible; since this ‘‘solution’’ is unique, it is the best 
smooth approximation to the true solution which can be obtained.* Thus 
these methods of making the resolution correction are of approximately the 
same usefulness, and are subject to essentially the same limitations, as the 
derivative series in equation (22) 

Let us now consider the method of numerical matrix inversion as outlined 
in equations (8) and (9). If the determinant of the matrix is not zero, one 
can find a unique set of numbers g; which reproduce the given f;,, no matter 


*Eckart (1937) considered kernels for which there is not a mathematically unique solution, 
and he therefore did not need to distinguish between true solutions and pseudosolutions. The 
distinction is made in the present case because the analytical solutions to equations (2) and 
(4) are unique. For numerical purposes, however, the existence of pseudosolutions destroys 
that uniqueness. 
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how many channels are chosen. But it does not follow in the limit of V going 
to infinity, that there exists an inverse kernel 7-'(V,E) such that 
T~"(V,E) = lim Tw 


Na 
and such that 
g(E) = [T—(V,E) f(V) av. 


Unless an inverse kernel does exist, the g(Z) obtained by numerical matrix 
inversion must be regarded as a pseudosolution only. 

It appears unlikely that a real, inverse kernel exists for the integral equation 
(2), but we have not been able to demonstrate this with mathematical rigor. 
Our main evidence for this contention comes from two numerical examples of 
20 X20 matrices which we have set up to represent 


(24kE)—2 exp{ —(V—E)?/2kE}. 


For 1 Ss E s 20,1 S V S 20 and & = 0.26 channel, the determinant of the 
matrix was found to be 10-**; for k = 0.52, the determinant was found to be 
10-4, Since these dnietndeunie seem to vanish, it appears unlikely that there 
is an inverse matrix, and unlikely that there is an inverse kernel. We conclude 
that the method of numerical matrix inversion should not be used to make 
the resolution correction, unless of course the existence of a real inverse kernel 
can be demonstrated. 

The Fourier sine series which was introduced above to illustrate the existence 
of pseudosolutions can also be used to illustrate the kind of difficulties one 
encounters with matrix methods for the inversion of equation (2). If we 
combine equations (30) and (31), we obtain the inversion formula as an 
integral of the following form: 


(32) gn(E) = JS™(V,E) f(V) dV, 
where 


2-27 mp 
OYE aiVS , niE_keaj  tiE\ wp hE 
(33) S”(V,E) = 7 > sin 30 \sin 20 ~ 99 ©°S 90 § &XP e090 


The sequence S” clearly diverges as n increases. 

Although the method of numerical matrix inversion does not appear to be 
suitable for making the resolution correction, it can be used to correct for the 
escape of energy from the crystal. In actual practice one must choose the 
matrix channel width to include the whole photopeak. The resulting matrix 
then is triangular in shape, and the determinant is equal to the product of the 
elements along the main diagonal. Since none of these elements is zero, at 
least in the general case, the determinant does not vanish, and there is no 
difficulty with this method of inversion. 

Other authors have shown that while numerical methods can be used for 
the solution of Volterra integral equations of the first kind, of which the 
energy-escape correction is an example, and for the solution of integral 
equations of the second kind, there are difficulties with the numerical solution 
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of Fredholm integral equations of the first kind. Fox and Goodwin (1953) 
have given some simple examples of integral equations of this latter type for 
which numerical matrix inversion is demonstrably invalid. There apparently 
still remains the general mathematical problem of finding useful conditions 








which are both necessary and sufficient for the existence of an inverse kernel 
for integral equations of this type.* 


IV. CONCLUSIONS 


Analytical solutions have been given for the integral equation (2) describing 
the effect of the statistical spread in pulse height. In order to make use of 
equation (22), however, it is necessary to evaluate derivatives by a numerical 
method. While the first few derivatives can be evaluated satisfactorily, the 
lack of definition in f(V) leads to large uncertainties in the high-order deriva- 
tives, and only a few terms in the series can be used. The method of fitting 
f(V) by a function for which the inverse transform is known (e.g. a Fourier 
series) encounters difficulties with pseudosolutions. The trial-and-error method 
also encounters difficulties with pseudosolutions, so that in practice one 
must restrict trial solutions by a condition of ‘“‘smoothness.”’ The trial-and- 
error method is somewhat more tedious to use than the derivative series, 
unless special means are used to effect the numerical integration. We have 
describéd elsewhere an analogue method for performing this operation (Aitken 
and Dixon). 

In practice there is no method of reconstructing by a resolution correction 
the fine structure of the true X-ray spectrum which the spectrometer failed 
to resolve. One can only sharpen the peaks and valleys already apparent 
on f(V). This, however, is often very useful. 

The method of numerical matrix inversion should not be used to make 
the resolution correction unless it can be demonstrated that there is an inverse 
kernel for the integral equation (2). This seems unlikely. The method can be 
used, however, to correct for the escape of energy from the crystal, if the 
matrix channel width is made wide enough to include the whole photopeak. 
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A PHENOMENOLOGICAL THEORY OF 
INDUCED ELECTRICAL POLARIZATION! 


JaMEs R. Wait? 


ABSTRACT 


A brief theoretical derivation is presented for the effective conductivity and 
dielectric constant of a homogeneous medium loaded with a uniform distribution 
of spherical conducting particles. To account for the effect of induced polarization, 
the particles are taken to have a concentric membrane or film which has a 
blocking action to the current flow into the particle. The characteristics of this 
phenomenological model are very similar to the experimentally observed 
features of induced polarization in a block of compacted andesite particles which 
contains a dissemination of small metal particles and is partially saturated with 
a weak electrolyte. 

The theory is then extended to a two-layer medium where the lower region is 
polarizable. The results explain, at least ina qualitative way, the observed features 
of induced electrical polarization in rocks, soils, and clay. 


INTRODUCTION 


The conduction of alternating electrical current through rocks and soils 
is a complicated phenomenon. It is known, however, that linear media can 
be represented at least in a macroscopic sense by an effective conductivity 
o, and an effective dielectric constant ¢,. That is, for a unit volume (e.g. a 
cubic meter) much larger than the pore and grain structure, the average 


> > 
current density J and the average electric field E are related by Ohm’s law 
in the manner 


> 


J = (o.+we,)E. 


This relation is appropriate for fields varying sinusoidally with time [i.e. the 
time factor is exp (iw)]. 

It has been experimentally observed that o,+iwe,, which can be replaced 
by &,., the complex conductivity, varies with frequency in the audio and sub- 
audio range in rather pronounced ways (Smith-Rose 1934; Evjen 1948; 
Keller and Lecastro 1955). This dispersion effect is due to certain electro- 
chemical effects which are not yet completely understood. 

It is the purpose of this paper to present a rather simple phenomenological 
theory for alternating current flow in soils and rocks. It is assumed that the 
microscopic structure can be represented by a homogeneous medium loaded 
with a distribution of conducting spherical particles, each coated with a di- 
electric film. Transient conditions are also considered. The experimental 
results of L. S. Collett are then interpreted in terms of the theoretical model. 
Finally, the theory is extended to a two-layer model of the ground where 
either or both the layers can be polarizable. 

1Manuscript received July 11, 1958. 

Contribution from Newmont Exploration Ltd., Danbury, Connecticut. 


*Present address: Central Radio Propagation Laboratory, National Bureau of Standards, 
Boulder, Colorado. 
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WAIT: ELECTRICAL POLARIZATION 


BASIC CONSIDERATIONS 

A single spherical particle of conductivity ¢,, radius a, has an insulating 
film of thickness ¢,, and dielectric constant €,. The suspending medium has a 
conductivity o. A spherical co-ordinate system (r, 0, ®) is chosen with the 
origin at the center of a single particle. The applied electric field Eye“ is in 
the polar (@ = 0) direction. The resulting potentials ¢, and ¢ inside and out- 
side the particles are solutions of Laplace’s equation since the significant 
dimensions in the problem are small compared to the wavelength. The boun- 
dary conditions at the film are that the normal current density is continuous 
and that the potentials are discontinuous by an amount equal to the voltage 
across the film. These two conditions lead to the following equations: 


a¢ | _ 5, ob | 
(1 ) ' or t=a+ tm a" or r=a 
and 
to ae | | 
9 sia ooelems — = 
(2) - 1€,,0 or r=a+tm % r=a 


where it is assumed that ¢,, « a. The exterior potential function is then found 
to be 


2, Eoa*| 1-6 
(3) o@ = —Eor cos "> 1125 cos 6, 
where 
go X54 ee. 
Sp 1€,,wd 


Now if there are V such particles per unit volume arranged in a regular 
lattice, the dipole moment 1 per cubic meter per unit applied field is 


(4) A = 3v(1—6)/(1+26) 


where 
v = 4ra°N/3 


is the volume of particles per cubic meter. The Clausius-Mossotti relation 
of dielectric theory (Fréhlich 1949) can now be employed to obtain an ex- 
pression for the effective conductivity o, and dielectric constant e, of the array 
of spherical particles, so that 


(5) oetiwe, = (124) , 


To illustrate this result it is assumed that the volume loading of the particles 
is small (i.e. v <1) and that the conductivity of the particles is large (i.e. 
a, > It then follows that 


5, Tetlwee * se) 
(6) o = 1+30( 18 ° 
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Consequently, 


es a 1—2x’ 
(7) g = i t8v ay? 
and 

EeW nN _ 9X 
(8) ‘¢ 1440" 
where 


X = tno/(wema). 
If X > 1 the particle behaves essentially as an insulator and 
(9) o,/o X= 1—3v/2 and ew/o = 9v/4X. 
If the admittance of the film is very large such that X < 1, 
(10) o,/o = 1+3v and ew/o = 9vX. 


In both these cases the dielectric constant expressed as a parameter e€.w/o 
is very small compared to unity. However, if the film has an admittance of 
an intermediate value, the effective dielectric constant can be relatively 
much larger. In fact, its maximum value can be seen to occur when X = 2 
and is given by 
(11) << | 180/17. 

oF max 
In other words, macroscopically speaking, the phase angle between the 
electric field vector and current density vector can be of the order of v radians 
for a suspension of film-covered conducting spherical particles if 


X =2. 


The elementary theoretical model presented here is not expected to have 
anything more than a semblance of actual conditions existing in mineralized 
media. The particles would not be spherical in shape, nor would the inter- 
face impedance be characterized by a simple lossy condenser. However, as 
will be indicated below, many of the features of induced polarization curves 
for rocks with disseminated sulphide particles can be explained by the above 
model, 

The foregoing discussion was confined to the case of a dielectric coated 
spherical particle. To generalize the results, it is convenient to introduce the 
interface impedance (w). It is the complex ratio of the voltage drop across 
the film to the normal current flow for unit area. In the above example 


l 2 
(12) n(w) = ~™- ohin meter 
LE», W 


where ¢, may also be complex to account for dissipation and can also be a 


function of frequency. 
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Assuming that o, >, it then follows that the effective (complex) con- 
ductivity ¢,(w) of the medium as a function of frequency is given by 
(13) ea soot _(l—v)+(2+0)¢ n(@)/a_ 
G.(w) wets (1+2v)+2(1—v)o n(w)/a’ 
TRANSIENT RESPONSE 
If a step function of (macroscopic) density Jo is now applied at ¢ = 0 to 


the medium, the build-up of the (macroscopic) electric field E(t), can be 
obtained from the Fourier integral (Churchill 1940). 


+oo 9 f@ 
(14) EW = 5. J Bl) ede = = [Re.E(w)] cos wt du 
where Re denotes the real part of 
; Jo 
E(w) = Ge(w) . tw 


Therefore, if the functional dependence of n(w) with frequency is known, the 
transient response of the medium can be predicted. To illustrate this, it is 


assumed that 


(15) ante) = ace where g = —m. 

This would mean that the film is behaving as a lossy condenser. is the low 
frequency, or static value of the interface impedance, and a is the relaxation 
constant for the interface. Furthermore, if the volume loading v is small, the 
effective conductivity equation can be ascot to 


; Se —gqa/(atiw) 
(16) a) a ga / (a+ iw) * 
Consequently 

" as Pee! ae )]| 
(17) E() = lwo 41 ot E (58. \ 


and when this is inserted in the integrand of equation (14) it readily follows 


that 


. 1— er 4 ( 3q ) a(l+2e)t 
e ¥ = — 9 A @ 
(18) E(t) Ja 1—3v 1429 — 32 1+29 e | 


for t > 0 and of course 
E(t) = Ofort < 0. 


Now it is interesting to note that the initial response is 
(19) E(0) = (Jo/o)[1—30] 


and the final response is 


» 1( co = - _ 
(20) K() Jel y 3v 142 29 


The time constant for the rate of build-up is [a(1+2¢)]}"! second. 
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Actually, the above formulas for E(0) and E() are valid for any functional 
form of n(w). This statement can be readily proved from the initial and final 
value theorems as used in the operational calculus (Churchill 1940). 

In most transient measurements of induced polarization in the laboratory, 
the quantity measured is the voltage e(t) following the interruption of steady 
current to the sample. In this case 


e(t) = V—E(t) 


where V = E(@) is often called the primary voltage. Therefore, 


(21) e(t) = 9v i+ exp [—a(1+2q)é]. 
The variation of the build-up curve E(t) with time is of course identical with 
the decay curve e(¢) apart from a change in polarity. This is characteristic of 
any linear system regardless of the functional form of n(w), the interface 
parameter. 

It is interesting to note that e(t) exhibits its maximum value as a function 
of g, when 


29 = one ad 
<) d= dm = 4, 1+ are P 


When the measurement is made near the bottom of the decay curve (at > 1), 
the maximum value of e(t) occurs when 


1 
at 
On the other hand, at very short times (at < 1) it follows that 
24 ns. 
To illustrate the features of the decay curve, the normalized response P(t) 
is plotted in Fig. 1 as a function of time, and in Fig. 2 as a function the particle 
size. The quantity P(t) is defined by 
e(t) 2. 
V 9v 


P(t) = 


and is thus normalized with respect to the primary field V and the volume 
loading v. 

It can be concluded from the above that a medium loaded with small 
particles will decay more rapidly than a medium with large particles (Vacquier 
1957). Furthermore, for a given time on the response curve, there will be an 
optimum particle size in the sense that the response is maximized. 


THE EXPERIMENTAL FACTS 
The experimental tests were carried out by L. S. Collett, using a large 


wooden box containing a tightly compacted mixture of andesite particles 
(98% by volume) and small pyrite particles (2%). To simulate conditions 


WAIT: ELECTRICAL POLARIZATION 1639 


100 . ——— 


VALUES OF a/yn,o 
SHOWN ON CURVES 





o 


003} 


Decay Response 


Oot 05 0 § 20 2 


Time Parameter, at 


Fic. 1. Normalized decay response P(¢) as a function of time. 
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Fic. 2. The normalized decay response P(t), as a function of the particle size parameter 
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that are typical of mineralized rocks, a weak electrolyte* (5% of the container 
volume) was allowed to be diffused in the wooden container. The andesite 
particles had diameters in the range 0.84 to 2.0 mm. Six different samples 
were prepared with pyrite particle diameters from 0.25 to 12 mm. 

The induced polarization decay curves shown in Fig. 3 are measured at the 
potential electrodes following the interruption of a 1-second current pulse 
at the current electrodes. (The curves are traced from a photograph of the 
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Fic. 3. Observed decay voltage e(¢) as a function of time. For V = 15 volts, v = 0.02, 


and o = 5X10- mhos/m. (To avoid confusion, three of the curve lines are broken.) 


oscilloscope screen.) The quantity plotted is actually e(t)/V. The diameters 
of the pyrite particles for each individual test are indicated on each curve. 
It is evident that the rate of decay for the smaller particles is more rapid 
than that for the larger particles in accordance with the predictions of the 
theoretical model. 

In Fig. 4 the instantaneous response at specified times (.01, .0316 second, 
etc.) is plotted as a function of particle size. It is quite evident that there 
is an optimum particle size for a given time. This optimum size has a slight 
tendency to be smaller for measuring times that are short in accordance with 
the theory. 

It is interesting to note that the optimum particle size is approaching 3 mm 
for the longer times. The interface impedance constant mo can then be solved 
from the equation 


om/a = 0.5 


where ¢o Y 5X10-* mhos/m to give m = 0.15 ohm meter’, which is of the 
order of magnitude measured by Adey (1950) and others, for interfaces 
between metal and weak electrolytes. 


*(.01 norma! NaCl solution and distilled water. 
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Fic. 4. Observed decay voltage e(t) as a function of particle diameter. For V = 15 volts, 
v = 0.02, and o = 5X107 mhos/m. 


TWO-LAYER MODEL 
It is worthwhile to extend the analysis to the case of a two-layer flat earth. 
The current J is to be applied to two point electrodes, C,; and C2, on the sur- 
face of the upper layer, and the resulting voltage V is to be measured between 
two additional electrodes, P; and P2 on the same surface. Using results from 
resistivity analysis (Sunde 1949), it readily follows that the mutual impedance 
between the electrode pair Ci, C2 and the electrode pair Pi, P2 can be written 


™ in ] , . 
(25) Z (iw) a esa aes a [G(C\P1) +G(C2P2) —G(C2P1) —G(CiP2)] 
27d;(w) 
where the G’s are functions of the individual linear distance between the 
electrodes. If the upper layer of thickness h has a (complex) conductivity 4:(w) 
and the lower semi-infinite layer has a (complex) conductivity ¢2(w) then 
1 = &1(w) —G2(w) 7 1 
26 G(x) = -+2 E enn I mgr 
(26) a(x) .* » G1(w) + G2(w) [x°+ (2nh)7}} 
The mutual impedance as a function of frequency is thus specified by the 
complex conductivities of the layers in the same manner that the mutual 
resistance is to the d-c. conductivities of the layers in the corresponding 
resistivity problem. It then follows that for the time harmonic problem the 
apparent conductivity ¢.(w) can be written 


1f 1 I 1 1 | 1 
9 = = — as ousapn _— ——— es 
(27) a(w) clapton, CP, CP.) Zia)’ 
In the case of a Wenner array, the electrodes are arranged in a straight 


line in the order C,P;P2C2 and spaced at equal intervals /. Then 


C;P, as C2P» = / and CoP, = C,P2 = 21 








1642 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


and therefore 
Ga(w) = 1/[2r1Z(iw)]. 


The voltage response B(t) at the potential electrodes P; and P: due to unit 
step of current (1 ampere applied at t = 0) is given (Churchill 1940) by 


1 +00 ew 1 +to st 
9 = i cae ee s. — 
(28) Bit) sc. iwZ (iw)? ~ Bri Stn 5215) 





Therefore, if the transfer impedance Z(iw) is known as a function of frequency, 
the step function response B(t) can, at least in principle, be determined. In 
practice it would be necessary to employ numerical methods of integration. 
To illustrate the transient computation, the two-layer model of the ground 
is considered where the conductivity of the upper layer o; is assumed to have 
a negligible induced polarization. That is o; is assumed to be substantially 
constant with respect to frequency (or time). It then easily follows that 


(29) A (t) = rae [g(CiP1, t)h+g(C2P2, t) —g(CiP2, t) —g(C2P,, t)] 


where* 
ee as ee 

(30) | g(x, t) = an oie as-=1, 1G(x::s) 
and 

i a ee eee 8 om 
(31) G(x, s) = 2122 [F(s)] ie (nh) | 
with 
a P(e) = 1 Fels) 
(32) F(s) = sitaa(s) 


When the lower layer consists of a dissemination of conducting spherical 
particles with interface characterized by an impedance n(iw) [or n(s)], then 


G2 oe 1—n(s)o2/a 


(33 a : 
d2(s) 1+29(s)o2/a 


a, is the conductivity of the lower layer in the absence of the particles of 
radius a, and it is assumed to be essentially constant with respect to frequency 
and time. As before, v is the fractional volume loading of particles in the 
lower layer. Again for purposes of illustration, the interface is taken to be a 
lossy condenser, and thence 

n(S)oe 


a Noe2 
= q- where g = 
a at+s a 


With these idealizations, the function F(s)/s can be conveniently written 
F(s) G 1 
(34 = —+(H—G = 
, 5 5 s+B 
*L is the Laplace transform operator (Churchill 1940) and L~' is the corresponding inverse 
operator (Churchill 1940 
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where 
F ae RN o1(1—3v) a2 
(35) H = F(o) oi day bes 
and 
q 
a(1- 30 tat) - o2 
(36) G = FO) = —— -- 
i—q 
a(1- wo tat) 46, 2 
and 
J E (2+ iter) 
oi Bm ee elles 
Furthermore 
[F(s)]? _ G* , 2G(H-G) _f)\?— 
(38) ; = > * s+6 +(H—-G)’ ora 
[F(s)]? _ G* , 3G°T-G) , 3G(H- G)" (H-G)’ 
( solid SN De 4 See 
) s sv stp * (s+~)? (s+s)° 


and so on. It is thus apparent that [F(s)]"/s can be expanded in reciprocal 
powers of (s+). Starting with the result 


(40) — = J ee "dt = Le *' 


and by differentiating successively both sides with respect to 8 and ¢ it follows 
that 


(41) Gap = L(1—pt)e*' 
(42) iss = L(1—2pt+8't"/2)e*' 


and so on. 
Using the preceding results we finally obtain 


(43) g(x, t) = it >» fall) 5 aes yh?) 
where 
(44) fult) = L een 


or more explicitly 
(45) fi(t) = G+(T—G)e*! 
(46) fo(t) = G?+2G(H—G)e'+ (H—G)*(1— Bie! 


(47) f(t) = G°-+3G2(H—G)e*'+3G(H —G)"(1— Bt)e** 
+ (H—G)[1 —28t+8°/2le 9 








1644 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


In every case, the temporal behavior is dominated by the exponential 
exp (—£$t) where 1/8 is the time constant. As seen from equation (37), if either 
v << lor if o> 01, then 


8 > a(1+2q) 


which is identical with the value derived in the previous section for an infinite 
medium with a uniform loading of spherical particles. However, if o2 is nearly 
equal to oj, it is seen that 


1+30/ 
B= of 1 -2¢( +20/4) | 


> a{l+2¢(1+90/4)]. 


It is thus apparent that the time constant of decay 1/8 can be shortened 
relative to the “‘infinite-medium value’’ [a(1+2gq)]-! in the case where az is 


near gj. 


CONCLUSION 
The model, while highly idealized, predicts, at least in a qualitative fashion, 
the observed features of induced polarization in an electrolytic medium con- 
taining metallic particles. (The formulas developed here will provide a basis 
for the subsequent discussion of extensive experimental data to be reported 
in the near future.) 
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ENERGIES OF GRAIN BOUNDARIES IN SILVER 
CHLORIDE CRYSTALS! 


KATHRYN A. McCartTuy? AND BRUCE CHALMERS 


ABSTRACT 


The relative energies of grain boundaries in silver chloride crystals have been 
measured. The measurements were made on tricrystals grown by simultaneously 
drawing three single crystals vertically from a melt. Tricrystals of three different 
isoaxial orientations, (100), (110), and (111), were studied. The relative grain 
boundary energies were measured as a function of the orientation differences 
between two grains forming a boundary in crystals of each of the orientations 
given above. From these measurements the following conclusions are drawn: 
the grain boundary energies of silver chloride follow the equation of Read and 
Shockley; for small orientation differences, that is, differences less than 10°, 
the boundary energy is a function of the axial orientation of the grains forming 
the boundary; for large orientation differences, that is, differences greater than 
10°, the boundary energy is not a function of the axial orientation, except when 
cusps occur; the maximum energy of a boundary made of crystals each with a 
(100) axial orientation equals the maximum energy of a boundary between 
crystals of (110) axial orientations. 


Most crystalline materials are found in polycrystalline form; specimens 
of such materials usually consist of many grains separated by grain boundaries. 
The energies of grain boundaries in several metals have been reported; the 
results are reproduced in Fig. 1. Aust and Chalmers (1950¢, 6) have studied 
the variation of grain boundary energies with the orientation difference of 
the grains forming the boundary, using tricrystals of tin and lead. Dunn, 
Daniels, and Bolton (1950) have made similar measurements on silicon iron. 
Greenough and King (1951) reported similar results on silver from measure- 
ments made on bicrystals. Shockley and Read (M49, 1952) have derived 
from first principles grain boundary energy as a function of orientation 
difference for a symmetrical boundary formed by two crystals of simple 
cubic structure. There is generally good agreement between the theoretical 
predictions of Shockley and Read and the experimental results presented 
above. 

Shockley and Read considered a dislocation model of a symmetrical boun- 
dary (in a simple cubic crystal) consisting of many edge dislocations infinite 
in length and parallel to the cube axis, about which each of the grains forming 
the boundary is rotated. They showed that the energy per unit area of the 
boundary, /, is a function of the orientation difference, 0, according to the 


equation: 


(1) E = E0(A-In@), 
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where Eo = ub/{4r(1—v)}, an energy term resulting from the elastic region 
associated with an edge dislocation, and A = a constant, that is a function 
of the inelastic region near the dislocation. 
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Fic. 1. Previous measurements of relative grain boundary energies as a function of orien- 
tation difference. 


The boundary energies reported below are relative energies determined 
by comparing the energy of a boundary that is a function of 6 with the energy 
of a boundary that is independent of 0. In a tricrystal the boundary energies 
are related to the tricrystal angles by the equation: 

2) You _ Yor _ Yoo. 

sina sing sin 6 
where the boundaries and angles are indicated in Fig. 2. yo4, yor, and yop 
are the boundary energies of the boundaries indicated by OA, OB, and OD, 
respectively. The relative boundary energies, 7,4, reported below were ob- 
tained by letting yo, equal unity. 

This paper reports measurements of the grain boundary energies of silver 
chloride crystals; the results are from measurements made on tricrystals. 
The technique for growing these tricrystals is essentially that used by Kyro- 
poulos (1926, 1930) to grow single crystals of the alkali halides. The crystals 


“yen rerser 
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were pulled vertically from the melt and were rotated throughout the growth 
process; all crystals used for the measurements reported in this paper were 
grown in air at a rate of 3.2 cm per hour. The silver chloride used in the 


D 





A 


Fic. 2. Diagram of a tricrystal surface. The three boundaries formed by the three possible 
combinations of crystals meet in a line perpendicular to the surface at O. 


melt was initially commercially c.p. silver chloride powder. Many single 
crystals were grown from such a melt; the melt from which the seed crystals 
and the tricrystals were grown consisted of remelted single crystals. The 
impurity content of the tricrystals was not determined; however, essentially 
no darkening of the silver chloride tricrystals occurred when the crystals 
were exposed to daylight. 

Silver chloride is very corrosive to metals. A container for molten silver 
chloride must resist corrosive attack up to temperatures slightly above 455° C, 
the melting point of silver chloride. It was found that a glazed porcelain dish 
was a suitable container; such a dish was quite durable provided that the 
melt was not cooled to room temperature and reheated, and that the melt was 
not heated excessively; in either of these cases, the glaze on the porcelain 
would crack. 

The pulling method is suitable for the growths of single crystals, bicrystals, 
and tricrystals. The crystal seed holders were made of stainless steel both to 
resist corrosion and to provide a means of cooling the seed. The holder for 
the growth of tricrystals is shown in Fig. 3. Each seed crystal has a separate 
support; as can be seen from the photograph, the orientations of the seed 
crystals about a fixed vertical direction can be changed by rotating the 
sleeve support in the main holder. The seed crystals were mounted so that 
the separation between the seeds was at least equal to the diameter of the 
seeds. It was found that tricrystals could not be grown when the seed crystals 
were in contact or nearly in contact with each other. The boundaries in the 
tricrystals grown with this method were in equilibrium positions; when the 
crystals were annealed, there was no boundary migration. 
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Fic. 3. The seed holder for the preparation of tricrystals of silver chloride. 


All tricrystals on which measurements are reported here were annealed for 
4 hours at 440° C, and were cooled in the furnace for approximately twenty 
hours. However, annealing had little effect on the tricrystals; in only two of 
the 49 tricrystals did annealing cause a change in the tricrystal angles. The 
tricrystals were annealed for two reasons: first, to reduce the number of sub- 
boundaries, and secondly, to remove the slight formation of color centers due 
to irradiation by X rays used for determining the orientations of the crystals. 

Silver chloride does not cleave, is difficult to cut, and deforms markedly 
even when cut with a string saw. For these reasons, when the tricrystal had 
grown to a sufficient length for measurement, the temperature of the melt 
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Fic. 4 


A thermally etched tricrystal surface (800 magnitication 


of the crystals is (110). 


Fic. 5. 
of the crystals is (100 


A chemically etched tricrystal surface (800 X% maguitication 
The groove at the boundary is clearly shown 
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was increased so that the crystal would lose contact with the melt; the crystal 
was then held just above the melt for about five minutes to permit thermal 
etching. A thermally etched tricrystal face is shown in Fig. 4. A chemically 
etched tricrystal face is shown in Fig. 5; the etch used was Kodak acid fixer. 

The data necessary to determine relative boundary energies are the relative 
orientations of the crystals forming the grain boundary, and the tricrystal 
angles. The orientation difference at a boundary is determined from back- 
reflection Laue X-ray photographs taken at points on the boundary; the 
precision of these measurements is estimated to be +1°. The tricrystal angles 
are measured from photomicrographs of 300 magnification, examples of 
which are shown in Figs. 4 and 5; the precision of these measurements is 
estimated to be +1°. 

Three groups of tricrystals were grown with isoaxial orientations (100), 
(110), and (111), respectively. For the crystals with (100) orientations, the 
orientation difference across the boundary was varied from 4° to 45°; the 
reference boundary, OB(100), was made of crystals with an orientation 
difference of either 33° or 37°. For the crystals with (110) and (111) orienta- 
tions, the orientation difference was varied from 5° to 76°; the reference 
boundary, OB(110), consisted of crystals with a difference of either 35° or 
47°, and the reference boundary, OB(111), consisted of crystals with a differ- 
ence of either 35° or 53°. The results are given in Fig. 6. These curves for 
silver chloride are similar in shape to those for the metals, summarized in 
Fig. 1. 

A cusp at an orientation difference of 45° is shown in Fig. 6 for a tricrystal 
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Fic. 6. Relative boundary energy as a function of orientation difference for tricrystals 
of isoaxial orientations (100), (110), and (111), respectively. 
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of (100) axial orientation; the energy versus © curve for this orientation is 
symmetrical about the 45° position. For the other two axial orientations 
studied, the energy versus 0 curves are symmetrical about 90°. 

A coherent twin boundary should be found for crystals of a (110) axial 
orientation so oriented that their common boundary plane is a {111} plane; 
four attempts to grow such a boundary were unsuccessful, leading to the 
conclusion that the cusp at the twin orientation is narrower than one degree. 

All the data reported above are of relative boundary energies. Before 
considering the conclusions to be drawn from the yrq versus 8 curves, some 
knowledge of the relative energies of the reference boundaries of each group 
of tricrystals is necessary. Thus an attempt has been made to compare experi- 
mentally the absolute maximum energies of a boundary of (110) axial orienta- 
tion and one of (100) axial orientation. A bicrystal formed by crystals of 
(100) axial orientation and of a 35° orientation difference was grown. This 
bicrystal was then bent, and used to seed a bicrystal of (110) axial orientation. 
The (110) bicrystal was then used as a seed together with a (110) single 
crystal to form a tricrystal with one boundary containing the same boundary 
as that of the (100) bicrystal, and another boundary with a 40° orientation 
difference. The result of the measurements of these tricrystal angles is that 
the maximum energy of the crystals of (100) axial orientation is equal to the 
maximum energy of the crystals of (110) axial orientation. Furthermore, the 
energy of the reference boundaries for the (100) crystals is equal to the energy 
of the reference boundaries for the (110) crystals. From the y;e versus 0 
curves it can be concluded that for small-angle boundaries, the boundary 
energy for the (100) crystals is less than that for the (110) crystals. Although 
attempts to do a similar experiment to compare the maximum energy of a 
(111) boundary with either a (100) or a (110) boundary were unsuccessful, 
it is concluded that for boundaries in silver chloride having a small-angle 
orientation difference the boundary energy is a function of the axial orientation 
of the crystals forming the boundary. 

Dunn’s resu'ts suggest a similar effect. His measurements (plotted in 
Fig. 1) on silicon iron are for two different axial orientations, (100) and (110), 
respectively. When his data are normalized to a relative boundary value of 
unity, they show that the boundary energy versus © curve for the (110) 
crystals has a slightly greater slope for small angles than does the curve for 
the (100) crystals. 

From equation 1, it follows that y;./@ should be a linear function of In 6. 
In Fig. 7, log 0 is plotted as a function of y;«/0 for the crystals of each orien- 
tation. The dats ‘or the crystals with (100) orientation are best fitted by two 
straight lines, one for angles less than 13°, and the other for angles greater 
than 13°; the departure from linearity above 37° results from the existence 
of the cusp at 45°. The relationship for the same functions for the crystals with 
(110) and (111) orientations is essentially linear. Thus it is concluded that 
for small-angle orientation differences, an equation of the same form as the 
equation of Shockley and Read applies to the grain boundary energies of 
silver chloride. The differences in the slopes of these lines for small values of 
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® show that when the dislocation model is applicable, the boundary energy 
is a function of the axial orientation of the crystals forming the boundary. 












































FIG. 7. yrei/@ as a function of log © for tricrystals of (100), (110), and (111) axial orienta- 
tions. 
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A MICROWAVE ANTENNA WITH RAPID SAW-TOOTH SCAN! 


J. S. Foster? 





ABSTRACT 


In one of many forms of this device, a parallel beam of microwave radiation, 
originating within the rotor, passes between coaxial rotor and stator conical con- 
ducting surfaces. The successive wave fronts consequently show a uniform 
rotation about a fixed stator barrier. This barrier removes the radiation which 
is subsequently directed to the linear focus of a parabolic cylinder serving as 
the antenna. The resulting rapid saw-tooth scan covers a range nearly equal to 
x times the total conical angle. The chief advantages of this system are (1) high 
speed of scan, (2) time between successive signals from any object is constant, 
(3) insensitivity to changes in magnetron frequency, (4) beam of constant form 
and intensity, (5) occurrence of the scanning at the exit barrier only, hence there 
is no side slap inside the conical structure, (6) relatively high efficiency. 

The beam must be properly directed through the conical structure. Scanners 
cannot be placed in tandem to get multiple angular ranges. A simple scanner is 
limited to a range of about 45°, but a more complex system is described which 
will cover 80°. 


INTRODUCTION 


a a eT TT I ET LTTE 


From an earlier description, given in some detail by Fry and Goward (1951), 
it appears that the physical ideas which form the basis of a rapid radar scanner 
I devised many years ago (Foster 1945; Millett et al. 1946) have not been 
understood by these authors and possibly some of their readers. The scanner 
represents, rather closely, a direct translation from a plane optical mirror 
rotating in a parallel beam to a highly versatile radar antenna with less loss 
of time between scans. The ideas are readily illustrated in the following de- 
scriptions of a limited number of early forms. At the same time, it may be 
important to recall some known limitations in the method. 


ORIGIN OF NEED FOR SCANNER 

The need occasionally arises for an antenna which will scan a sharp beam 
so rapidly as to permit the display on a fast screen of a motion picture of the 
distant scene of special interest. This need initially arose from the difficulty 
of manually keeping a 3/4° beam on a small plane 20 miles away. The reflector 
in use at that time was half of a parabolic cylinder and the line source at its 
focus was the mouth of a parabolic cheese. This in turn was suitably fed at its 
focus by a simple K-band waveguide. A vertical section through the center of 
the system is shown in Fig. 1(A). Since the vertical height of the reflector could 
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be fixed to give a beam of desired angular range in the vertical plane, the 
main problem was to provide a rapid scan of the sharp beam in the horizontal 
plane. 


PARABOLIC 


CYLINDER 
—— 


LINE FOCUS’ 





- 
WAVEGUIDE FEED 






(A) 


i PARABOLIC RIB 


Fic. 1(A). This is a central vertical section through a parabolic cheese the mouth of which 
provides a line source at the focus of the cylindrical reflector. This reflector, like all others 
shown in this paper, is half of a parabolic cylinder. 

Fic. 1(B). This shows a similar cross section with a conical structure inserted to produce 
a saw-tooth scan of the beam, twice in each revolution. The (black) barriers which guide the 
radiation through the structure are rows of teeth, each of which serves as a choke. They are 
off-set so that those attached to the rotor may pass through the stationary barriers. 


Rapid vibration of the whole antenna about a vertical axis is prohibitive 
on mechanical grounds and the sinusoidal scan so produced is not the most 
attractive. 

It was immediately clear that the above mechanical problem would dis- 
appear if the reflector remained fixed while the beam striking it was somehow 
scanned, preferably in saw-tooth fashion. 


INITIAL IDEAS FOR RAPID SCANNER 


The first thought for this scanner was thus directed toward the possibility 
of pulling the ribbon of radiation from the cheese round through a desired 
angular range before it reached the line focus of the reflector. For this purpose 
it seemed attractive to wrap the radiation round a structure like a conical 
clothespin; as the pin turned, the path would be longer at the larger end and 
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the wave fronts would swing round while remaining straight. In practice, 
one must have an outer stationary conical shell open along one straight 
line element to receive radiation from the parabolic cheese and open along 
a diametrically opposite element to let the radiation out to the reflector. 
Coaxially spaced within this, the conical clothespin rotor would then periodi- 
cally scan the beam in saw-tooth fashion, twice at each revolution, provided 
means could be found to steer the radiation through the structure in the 
proper manner (Fig. 1(B)). 

Any rotor barriers which steer the radiation must pass through the corre- 
sponding stator barriers. This requirement, however, is readily met by making 
each barrier like a comb with deep teeth, each of which is a choke for the 
radiation employed. Such teeth must be separated by a gap ample to pass 
teeth of another barrier, but too small to pass the radiation. In the figures, 
representative teeth of all barriers are shown for explanatory purposes; 
actually only rotor or stator teeth may appear in a single cross section. 


FIRST LABORATORY MODEL 


To avoid expensive patterns and castings, the greater part of a preliminary 
model (designed to scan a beam 2 ft wide through 10°) was made by electro- 
plating copper on graphite-coated wax. Barriers were made separately and 
insertéd. This model scanned 60 times per second and obviously had certain 
“built in’’ features, viz. it was insensitive to small changes in the frequency 
of the radiation and the signals from a distant object were equally spaced in 
time. The rather surprising efficiency of this model ensured the success of the 
method. 

It should be noted that no operational troubles have ever been encountered 
with the teeth, which may be made with ample clearance for passage of other 
teeth and require not too critical clearance at the top. Suitably contoured 
barriers may be assembled in parallel and the teeth sawed all at one time to 
avoid cumulative errors of spacing. Efforts to do away with them by simply 
setting up grooves to offer high impedance have met with the greatest diffi- 
culties owing to energy running along the grooves. However, a successful 
application of this latter method (first suggested by J. Schwinger) has been 
reported by Honey and Jones (1956) for a particular setting in which the 
grooves meet the radiation at a constant angle. 


ENGINEERED MODEL 


A second model for a beam 4 ft wide to be scanned through 10° was therefore 
engineered under the supervision of Dr. E. M. Purcell. The two halves of the 
outer shell and the rotor were machined from aluminum castings and separate 
barriers were installed. Some flaws in the casting for the rotor, however, were 
removed by smoothly rounding the entrance to the rotor slit on one side. 
The scanner was consequently abandoned. To revive this model, I built up 
the faulty corner with black laboratory wax and had the rotor turned to 
give proper corners on the slit. The exposed wax was painted with DuPont 
conducting silver paint and the barriers inserted. This revived model sur- 
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vived many tests; perhaps the most impressive was the motion picture of a 
Bren gun bullet visible for 1.25 miles at Petawawa. 


MODEL FOR TWO BEAMS TO BE SCANNED OVER THE SAME ANGULAR 
RANGE 

Meantime a demand arose for a single rapid scanner which would handle 

two beams, and this required some new elements of design. If only the para- 

bolic cheese could be rolled up inside the rotor (and fed through a rotating 

joint) then the two exits from the stator (Fig. 2(A)) might be used to feed 

two reflectors tilted about their focal lines to fixed angular positions as 







CONICAL 
SCANNER 


i 
DEAD SPACE “PARABOLIC” RIB 


CYLINDRICAL ~ 
FEEDS 





Fic. 2(A). This shows in section a means of scanning two beams over the same angular 
range. The parabolic cheese is rolled up inside the rotor (see text) and barriers direct the beam 
under scan alternately toward one or the other of two reflectors. The reflectors are normally 
tilted round the focal lines to positions suitable for automatic following or other purposes. 

Fic. 2(B). This is a similar representation of a device to scan over an angular range of 80°. 
Linear waveguide arrays imbedded in the rotor are designed to send two beams in opposite 
directions and these are scanned through the same 40° range as they pass from the conical 
structure. However, the line foci and reflectors may be so arranged that their projections on 
a horizontal plane have a relative tilt of 20° and the two beams may then range through 0° to 
—40° and 0° to +40°, respectively, relative to the normal to the fixed exit barrier. 


required. It is immediately evident that if each beam is to be scanned through 
10°, then the full vertical angle of the cones must be ~20/r°, i.e. the cones 
must be made to scan a single beam through 20°. If such an arrangement can 
be made, then with only three barriers the scanner will operate properly. 
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Fic. 2(C). This is a practical form of 2(A). The upper photograph shows the scanner in 
operating position and the lower shows the same scanner with the reflectors folded for pro- 
tection during shipment. 


The central problem is: how to roll up a flat parabolic cheese of finite 
thickness on a conical surface with its mouth at any desired twist from a 
straight line element of the surface and have the thing function properly. 
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This problem seems complicated by the fact that the outer face will be 
stretched with reference to the inner and the relative stretch increases as one 
approaches the smaller end of the cone. There are, however, two facts which 
make the problem indeed very simple. The first is a limitation on what is to 
be a parabola. It is perfectly obvious that both top and bottom front edges 
of the rib cannot be parabolic and indeed, in close approximation, only the 
line where the median coaxial conical surface meets the rib needs to meet this 
requirement. The second fact is revealed when this imaginary median surface 
is moved parallel to the common axis of the cones an easily determined dis- 
tance until it coincides with the base cone (upon which the cheese is being built). 
Then we shall have the true parabola imprinted on the base cone. If now one 
comes back along the base cone, from any point on this parabola, a distance 
dcot a/2 toward the apex of the cone, one reaches a point on the bottom front 
edge of the rib and a plot of such points enables one to set up the rib (of height 
2d). The first cheese set up in this way gave smaller side lobes than commonly 
experienced with flat models. Figure 3(B) shows how the template was pre- 
pared on roofing copper sheet. 

Under extensive practical tests, the cast aluminum antenna built on these 
principles detected objects at a range appreciably greater than that realized 
with other devices using the same magnetron power. 


PASSAGE OF WAVES BETWEEN CONICAL SURFACES 

Meantime the proper course of the wave fronts between the conical surfaces 
had to be considered and it is just this aspect which has not been understood 
by Fry and Goward. Consider that the space between the rotor and stator is 
thin and let it be developed on a plane (Fig. 3(A)). Within a small factor 
concerned with the finite size of the teeth in the barriers, it will become clear 
that @ is the angle of scan. For maximum width of beam (and to reduce reflec- 
tions at the sides) the parallel beam should be directed as shown, i.e. per- 
pendicular to the bisector of the angle a. To ensure this, the front of the 
parabolic cheese built under the surface of the rotor (dotted) must be tilted 
at an angle a/2 from the generatrix, which is indeed the middle of the exit 
AA’ leading to the scanning space between rotor and stator. This is so, because 
the barrier which brings this about obviously acts as the analogue of a plane 
mirror in optics. 

Consider now that BB’ is a barrier in the fixed outer conical shell and let the 
figure otherwise be uniformly rotated about 0 in imitation of the turning of 
the rotor. The wave fronts initially meeting the fixed barrier at the angle 
—a/2 will uniformly shift in direction until the angle is +a/2 and conse- 
quently the beam emerging from BB’ will shift from a/2 to —a/2 with respect 
to the normal. That is to say, since BB’ acts as a second plane mirror, the 
beam will be scanned through the angle a and the angular range will be sym- 
metrical with respect to the normal to BB’. The device has a fortunate feature 
in that most of the radiation which is reflected from the teeth fails to find its 
way back to the focus of the parabola and so to the generator. 

If now one returns to the first model, in which the scan is exactly that of 
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Fic. 3(A). This shows a thin conical scanning space developed on a plane (ABA’B’). The 
thin lines indicate plane waves properly directed and the dotted lines represent the developed 
parabolic cheese (under) which emits parallel radiation reflected at the barrier AA’. 

Fic. 3(B). This indicates the method for finding the template for the front bottom edge 
of the parabolic rib (see text). 














a ribbon wound round a conical clothespin, and considers that the scanning 
space is developed on a plane, it becomes evident that the direction of wave 
fronts changes at barriers in a manner similar to the above and that this early 
model for double frequency of scan is less suitable for large angles of scan. 


ae 


LEADING CHARACTERISTICS OF THE SCANNER 


The following additional characteristics of the scanner are now evident: 
(i) there is only one suitable way to send a parallel beam between the conical 
surfaces and that is the one just described; (ii) apart from diffraction, the 
wave fronts have a constant length throughout the range of scan, thus tending 
to retain a beam uniform in shape and intensity; (iii) the scanning occurs at 
the exit barrier only and consequently there is no ‘‘side slap’”’ inside the conical 
structure. To these may be added the characteristics already noted, viz. (iv) 
very rapid repetition rate of scan; (v) insensitivity to changes in magnetron 
frequency; (vi) time between successive signals from any object is constant; 
(vii) relatively high efficiency. 

A single beam may thus be satisfactorily scanned through an angle as 
large as 45°, but beyond this the method gradually becomes inefficient and 
troublesome. 
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FOSTER: MICROWAVE ANTENNA 


MODFL FOR EXCEPTIONALLY WIDE ANGLE OF SCAN 


If, nevertheless, a special need arises for scanning through a much larger 
angle, e.g. 80°, it is still possible to employ the method with two beams, each 
making use of the full scanner. Adopting the general method outlined above, 
W. Gordy first tried two beams, but failed to put each beam properly through 
the conical structure. 

Let us consider, then, a suitable method of scanning through 80° by sending 
two beams through the scanner in opposite directions. Each beam will be 
scanned through 40°; but there are to be separate outlets and reflectors 
(Fig. 2(B)). Since the passage through the scanner is properly identical in the 
two cases (except for direction), each beam will be scanned through an angular 
range of 40° symmetrical with respect to the two essentially parallel exit 
barriers. Since now the reflectors, fed from the two exits, are themselves parallel, 
the two beams will make a double scan over the same 40° in scissor-like fashion. 

While it momentarily appears that no increase in range of scan has been 
gained by this complication, the fact remains that the two beams must still 
be handled in this manner inside the conical structure and consequently can 
produce only duplicate scans at the exit barriers. 

The responsibility for really wide angle scan therefore rests on the relative 
tilt of the two reflectors. The side slap following the exit barriers may be 
allowed for in the relatively inexpensive structures of the feeds and reflectors. 
Parallel plate feeds to the line foci may be slashed off at any line of focus to 
suit the reflectors, which, incidentally, have a relative tilt of 20° when the 80° 
range is covered. One reflector provides a beam which scans 0°-40°; the 
other, 40°-80°. These must be joined in the range-azimuth display. 

It is not convenient to roll up two parabolic cheeses inside the rotor of a 
wide angle scanner. If on the same level, they interfere with each other and 
it seems needlessly complicated to put them at different levels inside the 
rotor unless the complete stability of beam directions against magnetron 
frequency variation is a matter of prime importance. 

Linear waveguide arrays of resonant slots as perfected by Watson and 
Guptill (Watson 1947) were therefore employed to give the initial parallel 
beams. Suitable arrays, designed by E. W. Guptill, were made from silver 
K-band waveguide with the upper broad side reduced to 0.015 in. thickness. 
The spacing between machined resonant slots, located alternatively on 
opposite sides of the middle line of the waveguide, was determined by the 
phase velocity in the guide and the angle (20°) which each emergent wave 
front is to make with the waveguide. A roughly triangular distribution of 
energy emerging from the array was realized by careful control of the offset 
of slots from the middle line and the residual 10-15% of the magnetron 
power was dumped into an absorber at the end of the array. Many of these 
details were worked out by Dr. J. W. Dodds. 


LIMITATIONS 


It should be unnecessary to point out that it was understood from the 
earliest days that these scanners cannot be put in tandem because of side 
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slap and that for the same reason a tandem effect cannot be realized in a 
single scanner by a barrier-mirror in a rotor working between two fixed shells. 
These limitations are indeed brought to mind at this time only because Fry 
and Goward have neglected them in proposed modifications of the scanner. 
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HORIZONTAL MOTIONS IN RADAR ECHOES FROM AURORA' 


G. F. Lyon AND A. KAVADAS 


ABSTRACT 


A systematic motion of 48.2 Mc/sec echoes associated with aurora is found 
at Saskatoon. The motion is towards the west before midnight and towards the 
east after midnight, the mean velocity in either direction showing a statistical 
relation to variations in the earth’s magnetic field. No correlation is found 
between individual echo velocity and magnetic disturbance, and no period of 
zero velocity is observed. There is also evidence of an ordered relation between 
ene in the north-south direction and disturbances in the earth’s magnetic 
field. 


INTRODUCTION 


A number of attempts have been made to detect systematic motions in the 
upper atmosphere. Windlike motions at heights of about 100 km were found 
by Mitra (1949), Briggs and Spencer (1954), and others who studied the 
motion of inhomogeneities in the regular ‘‘E”’ layer. Similar motions were 
found by Greenhow (1952), who studied the motion of meteor trails. Photo- 
graphic observations of aurora by Meinel and Schulte (1953) at Yerkes 
indicated a westward motion before midnight and an eastward motion after 
midnight, while Meek (1954) at Saskatoon found westward motions occurring 
up to 2 a.m. local time and eastward motions at any time during the night. 
Recent more extensive studies by this method by Kim and Currie (1958) 
using cameras at Saskatoon, Flin Flon, and Uranium City have established 
that westward motions occur in the early evening and eastward motions in 
the late morning while motions in either direction occur in the 4 hours centered 
on midnight. 

A radar investigation of aurora by Aspinall and Hawkins (1950) at Jodrell 
Bank, Manchester, revealed motions in individual echoes. In a later more 
detailed study at the same location Bullough and Kaiser (1955) found a well- 
defined minimum in the frequency of occurrence of echoes about 2130 hours 
U.T. with a westward motion before and an eastward motion after this 
minimum. 

This paper concerns an attempt to detect systematic motions in radar 
echoes from aurora. The observations were made at Saskatoon, Canada 
(latitude 52.1° N., longitude 106.6° W.), where a similar study using photo- 
graphic methods has been in progress (Kim and Currie 1958). 


EQUIPMENT 
The basic equipment used was a receiver with its associated circuits and a 
switching unit. Two separate receiving antennae were built, each consisting 
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of a vertical dipole in a 90° corner reflector. The principal maximum of sensi- 
tivity of one antenna is directed 45° to west of geomagnetic North and that 
of the other 45° to the east of geomagnetic North. No transmitting equipment 
was installed, instead advantage was taken of the presence of an I.G.Y. 
auroral radar already in continuous operation. The transmitter pulse was 
also used to trigger all operations at the receiving end. 

The I.G.Y. radar, which has been described by McNamara (1958), transmits 
300 microsecond pulses at a frequency of 48.2 Mc/sec with a pulse repetition 
rate of 70 per second. The transmitting antenna is a vertical dipole which 
radiates equally in all azimuthal directions. 

Signals from the two receiving antennas are alternately displayed, by means 
of the switching unit, on a conventional range-time, intensity-modulated, 
double-beam oscilloscope, the time base for the west antenna being vertically 
upwards from the center of the screen and that for the east antenna vertically 
downwards. Range markers are displayed at 100-km intervals and the display 
is photographed on continuously moving 35-mm film. The film moves at the 
rate of 12 inches per hour and is also exposed to a clock face every 15 minutes. 

Any radial motion in the reflecting center from which the echo is received 
will be recorded as a change in range as time advances. The slope of this 
echo on the range-time record will give a measure of its radial velocity. If the 
echo appears from one antenna only, the direction of motion may be identified 
as either receding or approaching in the NE. or the NW. geomagnetic quadrant 
respectively. If the echo appears first from one antenna and then from the 
other, the motion may be identified as westward or eastward. Finally, if the 
echo is receding or approaching on both antennas simultaneously, the motion 
may be identified as northward or southward. 

Previous experience with Plan Position Indicator type displays (Currie, 
Forsyth, and Vawter 1953) indicates that echoes at this frequency are obtained 
at Saskatoon through all azimuth angles of interest. 


OPERATIONAL DETAILS 

The equipment was operated as continuously as possible during day and 
night hours for a total of 77 days made up of the following periods: January 
18 to March 23, April 9 to April 18, and April 25 to April 28, 1958. The gaps 
in this operation were caused by severe local interference. 

The period covered was one of fairly high auroral activity. Out of the total 
operating period of 77 days echoes were observed on 59 days, the total echo 
time being 2834 minutes. Both discrete and diffuse echoes were observed, but 
owing to the fact that the received signal was amplitude limited it was not 
always possible to distinguish between the two types of echo; and no attempt 
was made to do so in the analysis of the results. 

The duration of individual echoes varied considerably. Some echoes lasted 
for only 1 minute or less, others for more than an hour. It was not always 
possible to assign a direction and a velocity to an echo. About one-fifth of 
the echoes were either of such short duration or exhibited such deep fading 
that the records were difficult to analyze. For some echoes it was possible to 












































er 


1 
SS es re 


Fic. 1. Some typical radar records. 
(a) Fast moving, short duration echoes on east antenna only. 


(b) Part of very long duration echo on both antennas simultaneously. 


(c) Complex echoes between 500 and 900 km on east antenna only. 
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assign a direction but not to measure the velocity. In all, some 330 echoes were 
analyzed. Of this number, velocities and directions were assigned to about 200 
and directions only to an additional 60. 

Owing to the local environment a greater number of echoes was received 
on the eastern antenna than on the western antenna. 

Some typical records are shown in Fig. 1. 


RESULTS 


1. Diurnal Variation of Echo Occurrence 

The relative occurrence of echoes for each 1-hour interval, expressed as a 
percentage of radar operating time, is shown in Fig. 2 (a). The main peak in 
echo occurrence is between 2200 and 0100 hours M.S.T. A secondary peak 
also occurs between 0600 and 0800 hours. This diurnal variation is similar 
to that found by Currie, Forsyth, and Vawter (1953) for auroral radar echoes 
at Saskatoon. Figure 2 (b) shows the same data grouped in 2-hour intervals 
(solid curve) for comparison with the figures obtained by these workers 
(dashed curve). 
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_ Fic. 2. Diurnal variation of echo occurrence. Dashed curve is from data obtained by 
Currie, Forsyth, and Vawter at Saskatoon in 1953. 


The Saskatoon diurnal frequency variation is markedly different from that 
found at Manchester by Bullough and Kaiser, who report maxima at approxi- 
mately 1800 and 0100 hours with a well-defined minimum about 2130 hours 
local time. 


2. Correlation with Visible Aurora 
The radar results were compared with the visible aurora for the same period 
as recorded by an all-sky camera, of the type described by Meek (1955), 
situated at Flin Flon, a distance of approximately 450 km in a NE. direction 
from the radar. 
The Flin Flon camera records were available for the period of radar opera- 
tion between January 18 and March 23, 1958. During this time 237 individual 











1664 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


echoes were received on the east antenna and of these 46 were received during 
daylight hours and 65 were received during periods in which the sky over 
Flin Flon was overcast. Of the remaining 127 echoes for which a direct com- 
parison could be made, only six were received at times when no visible aurora 
was present on the camera record. 


3. East-West Motion of Echoes 

For this study, all echoes were discarded which were clearly moving in a 
N.-S. direction, i.e. approaching or receding simultaneously on both antennas. 
There remained 311 echoes to which a direction could be assigned. Unfor- 
tunately, the majority of these echoes appeared only on one of the antennas 
so that it was not possible to state definitely that any particular echo was 
moving along the E.-W. line. However, when a large number of echoes is 
studied it becomes apparent that a systematic motion does occur as detailed 
below. 

For each hour, the echoes were grouped into classes representing eastward, 
westward, no motion, and indeterminate motion. The percentage of echoes 
by hours for the first three classes is shown in Fig. 3. From these curves it may 
be seen that before local midnight the motion is almost wholly to the west; 
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Fic. 3. Percentage of echoes moving to the west, to the east, and stationary as a function 
of time of day. 


after midnight the motion is almost wholly to the east; and during the hour 
centered upon midnight motions in either direction are equally frequent. 
It is important to note that in the period near midnight the motion does not 
disappear. 

For the greater part of the night, particularly during the critical hours 
around midnight, the relative numbers of echoes in each class received on the 
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two antennas remained constant. But from 0400 to 0800 hours the echoes 
were received only on the east antenna. 

A scatter diagram of velocity against time of day for each individual echo 
whose velocity could be measured is shown in Fig. 4. This again reveals wholly 
westward motion in the evening hours and wholly eastward motion in the 
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Fic. 4. Scatter diagram; individual echo-velocity versus time. 











morning hours, while from 2400 to 0130 hours velocities in either direction 
and of any magnitude occur. The mean velocity of the echoes for each hour of 
the day is given in Fig. 5. A sufficient number of echoes was obtained in any 
1 hour between 1900 and 0800 hours to give fairly reliable mean velocities, 
the smallest number of echoes received in any 1 hour being 8 and the greatest 
24. This curve clearly demonstrates that the maximum of the mean westward 
velocity during the evening hours (6 km/min) is less than that of the mean 
eastward velocity during the morning hours (15 km/min). 










4. Correlation of East-West Motion with Magnetic Activity 

Radar echoes obtained with the eastward looking antenna were compared 
with the magnetic records obtained at Flin Flon. The X component (horizontal 
component in the N.-S. direction) of the earth’s magnetic field was used and 
the usual convention was followed, positive indicating the northward direc- 
tion of the component vector. The magnitude of this vector was referred to a 
quiet-day level, usually the undisturbed period about 1200-1400 hours local 
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Fic. 5. Mean hourly radial velocity of all echoes. 


time. Changes in this magnitude will be referred to as positive or negative 


disturbances. 
A strong correlation between auroral and magnetic activity has been 
pointed out by many workers in the past. In the present study, less than 1% 


of the echoes were received at * .1es when the magnetometer showed no dis- 
turbance. Most of these occurred at about 2300 hours when the magnetic 
vector was changing rapidly from a positive to a negative value. 

The daily pattern of magnetic activity on disturbed days at Flin Flon is 
very regular. In general the disturbance is positive in the evening hours and 
negative in the early morning hours, the transition taking place around 
local midnight (2300 hours M.S.T.). The negative disturbances are almost 
always of greater magnitude than the positive disturbances. From the X- 
component records for each of the 59 days on which echoes were received, 
the positive or negative magnetic disturbance was read off at 15-minute inter- 
vals. The mean disturbance for each hour of the day averaged over the whole 
59 days was calculated. This diurnal variation is given in Fig. 6, with the 
time scale shifted to Saskatoon local time. The shape of this curve bears a 
very marked resemblance to the curve of diurnal variation of mean echo 
velocity (Fig. 5). The second curve in Fig. 6 is the diurnal variation of mean 
echo velocity for echoes received on the east antenna only; and this curve, 
within the limits of experimental accuracy, is identical with the diurnal 
variation of mean magnetic disturbance. The slight time displacement between 
the two curves may be accounted for by the fact that the reflecting region is 
some 250 km NE. of the magnetometer and the fact that as the negative bays 
are greater than the positive bays the averaging process shifts the magnetic 
curve to the left. These diurnal variations and the scatter diagram of velocity 
versus magnetic disturbance for individual echoes (Fig. 7) show that the 





LYON AND KAVADAS: RADAR ECHOES 1667 






+ 
100 10 
» w 
° 
Z 50 5 
a VELOCITY CURVE z 
ec — 
5 Kw = 
® oO = 
a ° x 
o ~ 
e 
= 80 53 
- MAGNETIC A °o 
2 CURVE Pe 
> > 
100 10 
150 5 
E 
200°—"6 18 20 22 00 02 04 06 oo. 


MOUNTAIN STANDARD TIME 


Fic. 6. Mean hourly X-component magnetic disturbance at Flin Flon, and mean hourly 
radial velocity for east antenna echoes only. 
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Fic. 7. Scatter diagram; individual east antenna echo velocity versus magnetic disturbance 
at Flin Flon. 


westward moving echoes are associated with positive magnetic disturbance 
and eastward moving echoes are associated with negative magnetic distur- 
bance. Exception to this rule was found in only five cases out of 122. 

The similarity of the two diurnal variation curves shows a statistical relation- 
ship between the velocity and the magnetic disturbance. On the other hand 
examination of the scatter diagram (Fig. 7) reveals that there is no apparent 
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simple relationship between the magnitude of the velocity of individual echoes 
and the magnitude of the magnetic disturbance. 


5. North-South Motions 

The number of echoes exhibiting a definite north-south or south—north 
motion was not sufficient to allow any statistical analysis. Nevertheless, some 
interesting points arise from the observations. Generally such motions occur 
when the magnitude of the magnetic disturbance is high, but there does not 
appear to be any correlation between the direction and speed of motion and 
the sign and magnitude of magnetic disturbance. 

Echoes with a definite N.-S. motion were observed on 12 occasions. Of 
these five moved first to the north and then reversed direction to the south, 
one moved first to the south and then to the north, two moved to the north 
only and four moved to the south only. All occurred between 2130 and 0350 
hours M.S.T. It was noted that when an echo changed its direction of motion 
from northward to southward, or vice versa, the rate of change of the X 
component reversed its sign but there was no connection between the sign 
of the slope and the direction of motion. A typical example is illustrated in 
Fig. 8. This occurred on January 20, the echo first moving to the north with 
the magnetometer positive increasing and then moving to the south with the 
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Fic. 8. Sample magnetogram (X component at Flin Flon) with north-south echo motion 
shown on same time scale. 


magnetometer positive decreasing. On February 14, an echo moved first to 
the north and then to the south while the magnetic disturbance changed from 
negative increasing to negative decreasing while on February 18 the motion 
changed from southward to northward while the magnetometer changed from 
negative increasing to negative decreasing. 


SUMMARY AND DISCUSSION 

The echoes were received apparently from reflecting regions of high ionic 
density associated with visible aurora. A good proportion of the echoes observed 
showed no significant radial velocity. While an echo could move in a circular 
path around the radar and appear stationary such a path seems highly im- 
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probable. It is therefore to be assumed that some ionic irregularities are in 
fact stationary at least for short periods of time. 

When motion does occur it is to the west in the evening hours and to the 
east in the early morning hours with a period centered around local midnight 
during which motion in either direction is equally probable. The proportion 
of echoes remaining stationary in this period around midnight appears to be 
the same, or if anything slightly less than at any other time during the night. 
The westward motion is associated with positive magnetic disturbance and 
the eastward motion with negative magnetic disturbance. 

The results obtained as to direction of motion agree well with those obtained 
optically at Saskatoon by Kim and Currie though the period of time during 
which motion in either direction can occur is narrowed to 2 hours from the 
4 hours obtained by these workers. For this work the radar method has the 
advantage of being able to cover any given period more thoroughly than 
can the optical method as the radar continues to gather data through both 
overcast and daylight hours when the camera cannot operate. 

The directions of motion found also agree well with those found at Man- 
chester by Bullough and Kaiser with the major exception that these workers 
report a virtual absence of echoes (and consequently a virtual zero velocity) 
between 2100 hours and 2200 hours local time. The association of westward 
and eastward motions with positive and negative magnetic disturbances 
respectively is clear both at Manchester and Saskatoon. The Manchester 
period of minimum echo occurrence and the Saskatoon period of motion in 
either direction both occur at the time when the disturbance in the magnetic 
field is changing locally from positive to negative. Heppner (1954) in extensive 
studies of the relation between visible aurora and magnetic disturbance at 
College, Alaska, found that at the time of change from positive to negative 
magnetic disturbance the associated visible aurora follows one of two definite 
patterns referred here as a and 0. (a) It persists over the same location but 
undergoes a distinct change of form, or (0) it disappears or recedes northward 
from the location as the positive bay disappears and reappears in the north 
and advances southward as the negative bay forms. Heppner also states that 
at College, Alaska (situated just south of the auroral zone) pattern (a) occurs. 
more frequently than pattern (6). There would appear to be a very real 
possibility that the pattern generally found at Manchester, well south of the 
auroral zone is pattern (0), since it would give just such a minimum in echo 
occurrence as that found by Bullough and Kaiser. On the other hand, at 
Saskatoon, which is much closer to the auroral zone, the predominant patterr 
would be pattern (a) and no minimum in echo occurrence would be expected. 

The point at which the mean velocity curve passes through zero (Fig. 6) 
is 2340 hours or 1 hour before local magnetic midnight at 0040 hours. At Man- 
chester the reversal of direction took place about 2200 hours, again approxi- 
mately 1 hour before local magnetic midnight, which occurs shortly after 
2300 hours. 

The mean peak westward and eastward velocities of 6 and 15 km/min re- 
spectively are lower than those found by other workers. Kim and Currie 
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using the optical method at Saskatoon found the average velocity in either 
direction to be 29 km/min. These workers could not take any zero velocities 
into account when computing their mean owing to the limitations of the 
optical method, so that higher figures than those found in this paper would 
be expected. Bullough and Kaiser obtain a figure of 30 km/min in either 
direction at Manchester by the radar method. Approximately one-fifth of the 
echoes used by these workers in computing their average velocities appear to 
have had very small velocities while in the present instance the proportion is 
approximately one-third. 

The authors consider it significant that the westward velocities obtained 
are lower than the eastward velocities. The positive magnetic bays are in- 
variably of smaller magnitude than the negative bays so that it is to be ex- 
pected that the currents in the ionosphere associated with the magnetic bays 
shall be smaller for positive bays than for negative bays. The good statistical 
agreement between peak mean velocity and magnetic disturbance suggests 
that the echo velocity is proportional to the current (or rate of net transport 
of charge). While the echoes originate at E-layer inhomogeneities which ex- 
hibit a systematic motion, the motion of these inhomogeneities does not 
necessarily constitute an electric current. The reflecting clouds are of such a 
size that they must remain electrically neutral and the macroscopic motion 
of a neutral cloud would not constitute a current. Hence, a one-to-one corre- 
lation of velocity with magnetic field strength would not be expected. No 
correlation is found between individual echo velocity and magnetic disturb- 
ance but the statistical agreement between the two remains, which suggests 
the existence of large scale (possibly worldwide) currents. 

To fit the observed motions there would need to be a net drift of ionospheric 
irregularities to the west in the evening and to the east in the morning. It is 
interesting to note that such a motion can be caused by dynamo action 
(Wulf 1953; Vestine 1954). The regular E-layer drifts observed by Briggs and 
Spencer (1954) exhibit a semidiurnal variation which in the northern hemi- 
sphere is a wind to the south at 2100 hours and 0900 hours local time and 
to the north at 1500 hours and 0300 hours. These winds would produce a 
motion of the ionospheric irregularities of high ion density associated with 
aurora. A southward movement of a conducting region, cutting the lines of 
force of the earth’s magnetic field, would produce an electric field directed 
to the west. Similarly a northward wind would produce a field directed to the 
east. This would result in a westward motion of electrons around 2100 hours 
and an eastward motion of electrons around 0300 hours corresponding approxi- 
mately to the observed times of echo motion. It should be noted, however, 
that the drift velocities of the irregularities are not necessarily the same as 
those of the ionization as a whole. The drift velocity of an irregularity in the 
presence of winds and the geomagnetic field depends strongly on its relative 
ionic density with respect to the surrounding plasma (Weeks 1957). Therefore, 
only statistical agreement between directions of ionospheric winds and echo 
motions is to be expected. 

The morning peak in echo velocity occurs a little later than 0300 hours. 
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From 0400 to 0800 hours the echoes are obtained almost exclusively from 
the east antenna. It is possible that the subsidiary morning peak in echo occur- 
rence is due to meteor ionization, enhanced by sunlight. The theoretical 
maximum in sporadic meteors occurs at 0600 hours. This concentration should 
be observed in the south moving eastwards (Lovell 1954). 
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MEASUREMENT OF THE DIELECTRIC PROPERTIES OF LOW 
LOSS MATERIALS AT MILLIMETER WAVELENGTHS! 


A. G. MUNGALL 


ABSTRACT 


Measurement at millimeter wavelengths of the dielectric properties of low 
loss materials by a free space technique is described. The dielectric constant is 
determined from the Brewster angle, and the dielectric loss tangent from the 
attenuation coefficient measured at the Brewster angle of incidence. Results are 
given for bakelite at wavelengths between 5 and 10 mm. Details of the instru- 
ment, which was specifically designed for these measurements, are also given. 


INTRODUCTION 


Measurement of the dielectric constant and loss of low loss dielectrics at 
millimeter wavelengths can be accomplished by both waveguide and free- 
space techniques. In the former, the electromagnetic fields are confined within 
waveguides and changes in these fields due to the insertion of a dielectric 
sample of known dimensions are used to compute the dielectric properties 
of the sample. In the latter, the effect of a dielectric material on free space 
propagation is measured. Transmission of the electromagnetic waves through 
a sheet of the material and reflection from it may be used to compute the 
dielectric properties. 

There are disadvantages in both methods. Measurements of dielectric 
properties by waveguide techniques at high microwave frequencies involve 
the use of very small samples, and demand elaborate techniques of measuring 
standing wave patterns in the waveguide. On the other hand, free-space 
techniques require relatively large samples in order that plane wave propaga- 
tion can be assumed to take place, but do not require the measurement of 
standing waves in waveguides. At wavelengths shorter than 1 cm, advantages 
inherent in the use of a small sample with the waveguide technique begin to 
be outweighed by the disadvantages encountered in accurate measurement of 
standing wave ratio and length. At these wavelengths, the samples required in 
free-space methods are not too large for practical application in the laboratory. 
For example, at 1 cm wavelength, with a sample 10 cm square, the conditions 
for free-space transmission and reflection are approximated. Larger samples 
are of course preferable. 

This paper will describe a free-space technique involving the measurement 
of the Brewster angle for the determination of the dielectric constant, and a 
measurement of the transmission loss at this angle for the determination of 
the loss tangent. The theory of this is well known, and has been described 
fully by previous workers (Bachynski 1958; Rabenhorst 1955; Talpey 1953; 
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Montgomery 1947). The present work has been done at higher frequencies and 
involves che use of a millimeter wave goniometer designed especially for this 
purpose. The properties of either so. or liquid dielectrics can be measured 
with this equipment. 


METHOD AND THEORY 


Two measurements are made, one .o determine the dielectric constant and 
the other the loss tangent. 


A. Dielectric Constant 

It is known that the reflection ~» ficient of electromagnetic waves at a 
dielectric interface is dependent on the electric permittivity and magnetic 
permeability of the dielectric, on the orientation of the electric vector with 
respect to the plane of incidence and reflection, and on the angle of incidence 
itself. In particular, the reflection coefficient at the boundary between free 
space and a loss-free dielectric surface (magnetic permeability that of free 
space) is zero for polarization of the electric vector parallel to the plane of 
incidence when the angle of incidence equals the Brewster angle. 

Under these conditions, the dielectric constant k, is given by the following 
relation with the Brewster angle @, (Stratton 1941): 


k. = tan? Og. 


For materials with magnetic constants of free space and with a small 
dielectric loss tangent (tan 6), this relation is still valid within experimental 
limits of accuracy. For example, with a dielectric constant near 2.5 and a loss 
tangent as large as 0.2, the value of k, resulting from this relation differs by 
less than 2% from that calculated from the following exact relation, for 
minimum reflection coefficient 7. 

_ k* cos 0—V/(ke—sin’ 8) 


” Re* cos 0-+4/(ke—sin’ 8) * 
Where @ is the angle of incidence, and k,* the complex dielectric constant, 


k.* = k,(1—j tan 8). 








The value of k, is determined experimentally by varying simultaneously 
the angles of incidence and reflection, and plotting the reflectance as a function 
of angle. The angle for which a minimum occurs can then be taken as the 
Brewster angle and the value of k, computed directly. The condition that the 
radiation be polarized in the plane of incidence is satisfied by appropriately 
orienting the horns on the transmitter and receiver. For most dielectrics, 
the condition that the magnetic permeability approximate that of free space 
is usually satisfied. 


B. Loss Tangent 

Determination of the loss tangent requires the knowledge of both the di- 
electric constant of the test specimen and the attenuation per wavelength of 
radiation passing through it. The dielectric constant can be determined by 
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the method described and the attenuation per wavelength by measuring the 
transmission loss of the specimen. If the incident angle of the radiation im- 
pinging upon the surface of the dielectric is made equal to the Brewster angle, 
then no power is reflected. At this angle of incidence, therefore, measurements 
of the transmission loss, and sample thickness, allow the calculation of the loss 
tangent. The relation between these quantities is approximately given by 


tan 6 = — cere 


if tan 5K 1, and where Xp» is the free-space wavelength, and L is the loss ex- 
pressed in decibels per unit thickness of the sample. This expression is strictly 
valid only for small values of tan 6; however, even for values as large as 0.5, 
the error in tan é is less than 1% for values of k, near 2.5, typical of many 
plastic insulators in the microwave and millimeter wave region 

Experimentally, the dielectric constant is first measured, and is then 
followed by a measurement of the transmission loss of a sample of the dielectric 
under test. The sample is placed in the microwave beam between tive trans- 
mitting and receiving horns at the Brewster angle of incidence, and the re- 
ceived power measured. 

INSTRUMENT DESIGN 
An instrument was designed specifically for these measurements and is 


shown in Fig. 1. It consists essentially of two radial arms moving in a vertical 
plane over a horizontal sample table. The two arms and the sample table are 





Fic. 1. The millimeter wave goniometer. 
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firmly attached to a heavy base provided with four levelling screws. One arm 
carries the microwave transmitter and the other the receiver. On both arms 
slots allow movement of the various waveguide components and lenses along 
the optic axis of the system. The transmitter consists of a klystron, precision 
attenuator, and a lens-horn combination, producing an approximately 
parallel beam of radiation. The polar distribution of power from this com- 
bination is determined mainly by the diffraction pattern resulting from the 
lens aperture itself. The radiation leaving the transmitting horn assembly 
then travels effectively parallel to the radial arm and impinges upon a reflecting 
surface on the sample table, which is either the dielectric sheet under test in 
the case of dielectric constant measurements, or a metal plate for loss measure- 
ments. This surface is oriented perpendicular to a plane containing the optic 
axes of transmitter and receiver, and passes through the point of intersection 
of these two axes. For the case of liquid dielectrics the specimen surface is 
horizontal and consequently the transmitting and receiving axes must be in a 
vertical plane. For equal angles of the two arms, reflection then takes place 
along the axis of the receiver. For dielectric constant measurements, the 
reflected radiation passes directly into a receiver with a lens—horn combination 
identical with that on the transmitter. For loss measurements, the radiation 
reflected from the metal sheet is made to pass through the dielectric specimen 
interposed between the reflecting surface and the receiver at the Brewster 
angle of incidence. For dielectric constant measurements, the detector does 
not require calibration, and the angle at which minimum power is reflected 
yields directly the dielectric constant. For dielectric constants near 2.5 an 
accuracy of angle measurement of +.25° results in an accuracy of dielectric 
constant measurement of +2%. 

Since the dielectric loss is calculated from the diminution of received power, 
in this case either the transmitter or receiver must be calibrated, and the 
accuracy of the relative power measurement determines the accuracy of the 
value of the dielectric loss. 


EXPERIMENTAL RESULTS 

The instrument has been used for measurement of the dielectric constant 
and loss of a number of materials. In Fig. 2 is shown a typical curve of the 
reflectance of a bakelite sheet 9 in. X12 in. as a function of angle of incidence 
for radiation of 6.0 mm wavelength polarized parallel to the plane of incidence. 
In the figure, only that portion of the reflectance curve near the Brewster 
angle is shown, to demonstrate the sharpness of the minimum. Point scatter 
probably arises from a combination of interference effects due to non-uniformi- 
ties in the thickness of the specimen, and stray reflections from the instrument 
sample holder and base. A typical set of data for dielectric constant and loss 
as a function of frequency is shown in Fig. 3, for bakelite. Accuracy is within 
5% for both the dielectric constant and loss measurements. 

Ideally, in the measurement of the dielectric loss, with the specimen inserted 
in the received beam, there should be no standing wave at the receiving horn, 
so that the power loss due to the insertion of the dielectric can be measured 














CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


RELATIVE POWER REFLECTED 


60 6! 62 
ANGLE OF INCIDENCE 


Fic. 2. Reflectance of bakelite near the Brewster angle for 6-mm radiation polarized) 
parallel to the plane of incidence. 
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Fic. 3. The dielectric constant and loss of bakelite at wavelengths between 5 and 10 mm. 
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directly. In practice, owing to reflections at the lenses and horns, a small stand- 
ing wave is set up and is shifted axially when the specimen is inserted. This 
necessitates the measurement of the standing wave pattern with and without 
the specimen in order to ascertain the power absorption due to the dielectric. 
This can be done simply by moving the receiving horn in an axial direction 
and noting the power maxima and minima. For low loss specimens, standing 
wave ratios approaching unity are preferable: these can be obtained by 
inserting a lossy dielectric sheet between transmitter and reflecting surface. 
For short wavelength radiation, disturbing interference effects can result 
from small variations of sample thickness across the transmitted beam. For 
example, at 5 mm wavelength and a sample with a dielectric constant of 3, 
a variation in thickness of .057 in. over half the sample area illuminated by 
the microwave beam can result in a 180° phase shift between the two halves 
of the beam and a resulting null in detected power. Smaller thickness variations 
will not cause a null, but can result in appreciable errors of loss measurement 
of low loss specimens. The required degree of flatness and parallelism of the 
specimen can, however, be simply estimated for a given accuracy and the 
sample machined accordingly. 
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AIR-MASS REFRACTIVITY IN CENTRAL CANADA! 


D. R. Hay 


ABSTRACT 


A study of air temperature, humidity, and pressure profiles above Maniwaki, 
Quebec, for the year 1957 shows that the four types of air masses have distinctive 
refractivities. It is shown also that the air masses may be identified through 
their dry and total refractivities. For a typical air-to-ground communication 


w 


path, the mean effective-earth’s radius varies from 1.26 to 1.37 times the true 
radius for the four types of air mass. 


INTRODUCTION 

The importance of the troposphere in the refraction of short radio waves 
has been recognized for some twenty-five years. In 1931, Jouaust suggested 
that apparent changes in the height of the horizon along a radio path were 
due to variation in the weather near the earth’s surface. Schelleng, Burrows, 
and Ferrell pointed out in 1933 that the bending of radio waves by the average 
atmosphere near the earth’s surface could be simulated by assuming the 
atmosphere to be homogeneous and the earth to have an effective radius 
which is 1.33 times its true radius. A formal analysis of the concept of effective- 
earth’s radius was presented 4 years later by Eckersley (1937) and van der 
Pol (1938). 

Subsequently, refinements to this approximation have been made. Bean 
(1953) has found that in North America the definition of a standard atmosphere 
varies with geographic latitude. Further, over the U.S.A. the effective-earth’s 
radius factor to be used in tracing radio paths varies from 1.25 to 1.45 during 
February, and from 1.25 to 1.9 during August. Fannin and Jehn (1957) and 
Millman (1958) recently have demonstrated the importance of refractivity 
information for the troposphere in the radar tracking of distant targets. How- 
ever, little information of this sort is available for Canada’s climate. Grant 
(1954) has shown that the average air below 10,000 ft over northwestern 
Canada and Labrador during the summer may be classified as a N.A.C.A. 
standard atmosphere with a relative humidity of 60%. 

A more accurate classification of tropospheric refractivities? will be de- 
scribed in this paper. The need for this classification has arisen in a study of 
the effects of certain weather perturbations upon UHF radio transmission. 
It will be shown that each type of large-scale air mass in central Canada 
has a distinctive average refractivity profile. These average refractivities may 
be used to estimate the extent of the refractive index irregularity which 
accompanies a weather system. (The deviations about a mean refractivity 
are measured readily by sensitive refractometers (Crain 1955).) Alternatively, 

'1Manuscript received June 2, 1958. 
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the identity of an air mass in central Canada may be established from its 
refractivity profile. 
TOTAL REFRACTIVITY 

Refractivity profiles of specific air masses have been derived from soundings 
above Maniwaki, Quebec, of temperature, pressure, and humidity.’ This 
information was obtained for each of the four types of large-scale air mass in 
accordance with the three-front model used by Canadian meteorologists 
(Penner 1955; Anderson et al. 1955). (The identification of the air masses 
was part of the routine weather analysis carried out at the Central Analysis 
Office in Dorval.) In this region of Canada, the air masses are the Continental 
Arctic (cA), Maritime Arctic (mA), Maritime Polar (mP), and Maritime 
Tropical (mT). Each air mass has distinctive physical characteristics which 
are governed by its history. The temperature and dewpoint within an air 
mass are approximately constant along a horizontal plane; the transition 
in temperature between two adjacent air masses occurs at the frontal zone 
separating them. 

Fig. 1 illustrates the refractivity profile for each of the four types of air 
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Fic. 1. Total refractivity for air masses over central Canada. 


’Refractivity of troposphere = (7.76X10-5/T) {p+4810(e/T)}, where p = total air pres- 
sure in millibars, e = partial pressure of water vapor in millibars, and T = absolute air 
temperature (Smith and Weintraub 1953). This expression is valid to 0.5% in refractivity for 
radio frequencies up to 20,000 Mc/sec and for the ranges of pressure, temperature, and humidity 
normally found in the troposphere. Refractivity may be computed rapidly from measurements 
of temperature, air pressure, and dewpoint with the aid of a simple analogue technique (John- 
son 1953). The first term on the righthand side cf the above equation represents the contri- 
bution from dry air, the second term is the contribution from moisture in the air. 
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mass. Conventional radiosonde equipment was used for the temperature and 
humidity measurements. Information on the air profiles was obtained from 
the following radiosonde ascents made during 1957: 10 soundings during 
March, April, September, October, and November for the Continental 
Arctic air profile; 24 soundings during the months from March to November 
for the Maritime Arctic air profile; 40 soundings during the months from 
May to October for the Maritime Polar air profile, and 13 soundings during 
the months of June, July, and September for the Maritime Tropical air 
profile. The computed refractivity at heights between 0.5 and 6.5 km is given 
by: 


(n—1) = 10-* (310.00 —36.534+1.63h?)...Continental Arctic air 
10-* (321.20 —42.10h+2.17h?)...Maritime Arctic air 
= 10~-* (331.96 —44.33h+2.26h?)... Maritime Polar air 

= 10~* (346.27 —45.29h+2.14h’)... Maritime Tropical air 


where # = height in kilometers above mean sea level. It should be noted that 
these profiles apply to average conditions in each air mass. Although all of the 
available information from each radiosonde sounding was used in computing 
refractivity, details on local irregularities were removed in averaging of the 
profiles. 

Two features of the air masses may be noted in the curves of Fig. 1. From 
refractivity profiles alone, it would appear that the three Maritime types 
belong to a category which differs from that of the Continental type of air. 
Secondly, the average refractivity in the Maritime air masses decreases in 
the order Maritime Tropical, Maritime Polar, and Maritime Arctic. This is 
in the same order as their mean geographical positions from south to north. 
The N.A.C.A. standard atmosphere with 60% relative humidity, to which 
Grant (1954) has referred, has a refractivity profile which is identical with 
that for Maritime Arctic air at heights below 2.5 km. 


RADIO RAY BENDING 


One important application of the refractivity profile is in estimating correc- 
tions to the apparent path of a radio signal through the troposphere. If the 
radio path is through a known type of air mass, the average path may be 
computed by assuming the effective-earth’s radius which is appropriate for 
this air mass. For this purpose, the information on refractive index that was 
used to derive Fig. 1 has been reproduced in greater detail in Fig. 2. Local 
irregularities in temperature and humidity generally cause the refractive index 
to deviate from the mean profiles given in Fig. 1. The standard deviation 
about the mean is shown in Fig. 2 by the shaded areas about the solid lines. 
It will be noted that Fig. 2 is a composite diagram that shows a pair of profiles 
for each type of air mass. It is the upper curve of each pair, giving total 
refractivity, that is used to derive the effective-earth’s radius. 

Typical values of effective-earth’s radius may be derived .vr a radio path 
that is located mainly between 0.5 and 2.0 km in height. If the refractive 
index obeys the normal probability law of distribution about the mean profile, 
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Fic. 2. Total and dry refractivities for air masses over central Canada. 


then the extent of variation may be calculated from the indicated standard 
deviation. Finally, the effective-earth’s radius may be computed from the 
gradient of refractivity (Kerr 1951)* as shown in Table I. Since the refractivity 
gradient is not constant over this height interval, the average gradient is 
used for this calculation. 


IDENTIFICATION OF AIR MASSES 
The refractivity profiles will be examined further for the possibility of using 
them to identify air masses. It should be emphasized here that routine identi- 
fication of air masses is carried out efficiently through thermodynamic quanti- 
ties, such as wet-bulb potential temperature. Consequently, radiosonde in- 
formation is more useful for this specific purpose than is data indicated by a 


4 True-earth’s radius _ /a **) 41 
Effective-earth’s radius no dh , 

where a = true-earth’s radius, mo = refractive index of air at surface, and dn/dh = gradient 

in refractive index. 
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TABLE I 


Effective-earth’s radius 


Limits for 95% 
Air mass Mean value probability 


cA 1.26 a 1.24ato1.27a 
mA 1.32" i236 — tv” 
mP Lao” Loa”) Oa 
mT ion P30”. 1.40" 


refractometer. The increasing use of air-borne refractometers, however, does 
suggest the need for identifying air masses through refractivity profiles alone. 

Curves that show the component of refractivity due to dry air in each type 
of air mass are included in Fig. 2 to assist in the identification. These are the 
lower curves in each pair. Again, the shaded area represents the standard 
deviation°*about the mean. The equations of the dry refractivity profiles are: 


(n—1) = 10-* (289.23 —29.93h+1.07h?)...Continental Arctic air (dry) 
= 10-* (275.00 —27.10h+0.87h?)... Maritime Arctic air (dry) 
= 10-* (265.97 — 25.23h+0.79h?)... Maritime Polar air (dry) 
= 10-* (261.26 —23.56h+0.64h?)... Maritime Tropical air (dry) 


where again = height in kilometers above mean sea level. It is apparent 
that the magnitude of the contribution to mean refractivity from moisture in 


the air is greater for those air masses whose mean geographical positions are 
at lower latitudes than it is for those whose mean positi~1s are at higher 
latitudes. 

The curves of Fig. 2 suggest a procedure that may be used in air-mass identi- 
fication. If the four pairs of curves are superimposed, it will be seen that both 
dry and total refractivity profiles are needed to distinguish between the 
different types of air mass, because of overlapping in the standard deviations 
about the mean refractivity. However, if the information from an isolated 
radiosonde ascent is used to compute dry and total refractivities, the air 
mass may be classified with a high degree of probability under one of the 
following three headings, through the curves in Fig. 2: 


(i) Continental Arctic air, 
(ii) Maritime Arctic air, 
(iii) Maritime Polar or Maritime Tropical air. 


If the air mass falls into the third category, several sets of measurements may 
be required to separate the two possibilities. 


The author is indebted to the Director of the Meteorological Branch, Depart- 
ment of Transport, and to the Officers-in-Charge at the Maniwaki Radio- 
sonde Station and at the Uplands (R.C.A.F.) Weather Office for the weather 
information used in the preparation of this paper. He is grateful also to 
Dr. W. Godson and to Mr. J. Leaver for helpful advice on the meteorology. 
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THE QUENCHING OF TRIPLET POSITRONIUM BY IONS 
IN AQUEOUS SOLUTIONS! 


R. E. GREEN? AND R. E. BELL 


ABSTRACT 


Measurements of positron lifetimes in aqueous solutions have been performed 
using a fast time-to-amplitude converter with a resolving time of. ~10~° sec. 
The lifetime, 72, of the long-lived component of the complex decay observed in 
water has been found to decrease when paramagnetic salts are added. The experi- 
mental results are interpreted in terms of “‘pickoff” annihilation to explain the 
observed 72 in pure water, and a rapid triplet <> singlet conversion process to 
account for the quenching of the 72 component. The conversion is caused by the 
presence of unpaired electrons in the structure of the paramagnetic atoms. 
Triplet — singlet conversion cross sections have been calculated on the 
basis of this interpretation and they range from 4.0X10~* cm? for NdCl; to 
3.8X 107% cm? for FeCl. 


INTRODUCTION 


Measurements of the lifetimes of positrons in gases, liquids, and amorphous 
solids have shown that positronium exists in these substances. The work of 
Deutsch (1951) confirmed the existence of positronium in gases and the 
experiments of Bell and Graham (1953) showed that positronium is formed 
in liquids.and amorphous solids. The latter authors found that the decay of 
positrons in liquids and amorphous solids was complex in nature. Two distinct 
lifetimes were observed, the shorter, called 7;, being a few times 10~' sec, and 
the longer, r2, varying from ~5 X10-" sec to 3.5X10-* sec. The 72 radiation 
had a spectral distribution essentially the same as the 7; component, viz., a 
two-quantum annihilation spectrum. Since the triplet state of positronium 
can only decay by the emission of three photons (Yang 1950; Peaslee 1950), 
Bell and Graham suggested that the triplet positronium decayed mainly by 
conversion to the singlet state followed by a rapid annihilation from this short- 
lived state (7, = 1.25X10-' sec). The conversion was thought to be due to 
exchange collisions with the electrons of the neighboring atoms. Another 
and more reasonable explanation has since been suggested to explain the 72 
lifetime. This is the so-called “‘pickoff’”’ annihilation process, first suggested by 
Garwin (1953), in which the positron in a triplet positronium atom annihi- 
lates by two-quantum annihilation with an electron, from a neighboring atom, 
with a favorable spin orientation. The results of the present experiments 
favor this latter mechanism to explain the observed 72 lifetime in water. 
Recent experiments of Green and Bell (1957) also favor the pickoff mechanism, 
and we will henceforth take it to be correct. 

The quenching of triplet positronium has also been observed. Experiments 
by several workers (Deutsch and Dulit 1951; Pond and Dicke 1952) have 
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shown that the application of a strong magnetic field causes a mixing of the 
triplet m = 0 substate with the singlet state and so this component of the 
triplet positronium is quenched. The addition of small quantities of nitric 
oxide (NO) to a gas exhibiting the triplet positronium lifetime will completely 
quench out this lifetime. The mechanism is assumed to be triplet <> singlet 
conversion since the energy required to flip the unpaired electron of nitric 
oxide is small. Small quantities of the free radical DPH (diphenylpicryl- 
hydrazyl), which has one unpaired electron, when added to pure benzene, 
quenches the 72 lifetime (Berko and Zucheili 1956). Assuming that pure 
benzene is itself a sufficient triplet < singlet converter to account for the 
observed 72 lifetime, and that the added DPH further assists the conversion 
process, Berko and Zuchelli obtained a triplet — singlet conversion cross 
section for DPH of 1.18X10-" cm*. This result was obtained assuming 
that the positronium atoms were moving with thermal velocities. Berko and 
Zuchelli stated that the cross section was of the same order as that found by 
Deutsch (1955) for the nitric oxide molecule, about 10~-'* cm?. On this ground 
they concluded that the quenching effect of DPH is due to a spin exchange. 
Recent work by de Zafra (1958) on the angular distribution of annihilation 
radiation has cast doubt upon this conclusion, however, and he supposes 
that the effect of DPH in benzene is caused by attachment of positronium 
atoms to the molecules of the free radical. Berko and Zuchelli also observed 
a weak quenching by Mn** ions in water (quoted in Berko and Hereford 
1956), a result extended here. 

The purpose of this work has been to study the lifetimes of positrons in 
aqueous solutions and in particular to look for substances which would 
quench out the 72 component observed in water. It was hoped that para- 
magnetic salts containing unpaired electrons dissolved in water would produce 
a quenching effect comparable to that of DPH in benzene and nitric oxide in 
gases. Therefore, aqueous solutions of the iron group elements and the rare 
earth group elements have been studied. The results of these experiments show 
definite quenching effects, and cross sections have been calculated for the 
triplet — singlet conversion process assuming that the unquenched 7; lifetime 
is due to “‘pickoff”’ annihilation. Preliminary results from this study have been 
presented at the Electron Physics Conference, University of Maryland, 1956, 
at the Annual Meeting of the Royal Society of Canada, 1956, and at the 
Washington Meeting of the American Physical Society (Bell 1958). 


EXPERIMENTAL TECHNIQUE 

The source of positrons used in this work was Na”. The source (~ 10 uc) 
was sandwiched between two thin disks of mica of ~ 1.8 mg/cm? areal density 
and 13 in. diameter. (Mica was chosen since it is chemically inert and also 
shows no complex decay.) The mica disks were mounted on tantalum rings 
(tantalum being acid resistant), and the edges of the mica were cemented with 
Apiezon black wax to ensure that the source was waterproof. The source was 
then placed in a Lucite cell containing the solution under study, the size of 
the cell being such that all the positrons emerging from the source were 
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stopped in the liquid before reaching the cell walls. The cell was viewed by 
two diphenylacetylene-6342 scintillation counters connected to a time-to- 
amplitude converter and a 28-channel pulse analyzer. A single measurement 
yields a curve showing the distribution in time of the annihilation radiation 
from the sample measured with respect to the time of emission of the 1.28- 
Mev gamma ray from Na”. Further details are given in a recent paper (Green 
and Bell 1957). 





RESULTS AND DISCUSSION 





EXPERIMENTAL 
(a) Iron Group 

The first substances studied were aqueous solutions of the iron group ele- 
ments. These elements were chosen since they exhibit a strong paramagnetism 
due to the presence of unpaired electrons in their electronic structure. The 
unpaired electrons arise from the filling of the 3d subshell of these atoms. 
The 3d electrons do not pair off until they are forced to do so, i.e. until six are 
added, so that as many as five unpaired electrons (or holes) can be present. 
It was thought that these unpaired electrons would quench the 72 lifetime 
observed for positrons in pure water by causing a rapid conversion between 
the two states of positronium. 

The salts studied were chlorides. They were chosen since no quenching had 
been observed when the lifetime was measured in a strong NaCl solution 
(Bell and Graham 1953). This showed that both the sodium and chloride 
ions were non-quenchers, and so any quenching observed in chlorides could 
be attributed to the cation. Again the fairly large solubility of most chlorides 
made these salts convenient, since it was usually possible to make the solutions 
concentrated enough to observe complete quenching. 

The results of a series of lifetime measurements made in FeCl; solutions 
of various concentrations are shown in Fig. 1, in the form of experimental 
time distributions. The zero of time in this graph is arbitrary. The complex 
nature of the decay is evident and the variation of the 72 lifetime is quite 
definite. Other concentrations of FeCl; have been tried, but for the sake of 
clarity are not shown on this diagram. The 0.001 molar solution is essentially 
a non-quencher since the 72 lifetime is the same as that observed in water, 
i.e. 1.80X10~-* sec. The shortest 7. observed with FeCl; was ~4X 107" sec, 
the particular one shown here being 4.1 X 10~" sec. This is not an instrumental 
limit, but rather one set by the conversion mechanism. This is clear from the 
fact that there is a distinct difference between this curve and the aluminum 
curve also shown in Fig. 1. These results indicate that rz varies with concen- 
tration between the limits 4.0X10~ sec for concentrated solutions and 
1.80 X10-* sec for pure water. 

The measured variation of r. as a function of FeCl; concentration is shown 
in Fig. 2. On this diagram two theoretical curves are shown. The one labelled 
CONVERSION has been computed on the basis of water as a weak converter 
giving tz = 1.80X10~*% sec and the FeCl; as an additional converter. The 
other, labelled pickorr, is calculated assuming that pickoff annihilation 
directly from the triplet state is responsible for the 72 lifetime of water, and 
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Fic. 1. Experimental time distributions of positron annihilations in aqueous solutions of 
FeC]l;. The inner curve is for annihilation in aluminum. 
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that the quenching is due to triplet <> singlet conversion caused by the added 
FeCl;. The equations from which these curves have been derived have been 
written down by Dixon and Trainor (1955), and a brief outline of the method 
of calculation of these curves will now be given. 

Using the notation of Dixon and Trainor, let us assume that P, is the 
probability of finding a positronium atom, once formed, in the singlet state 
and P, the corresponding probability for the triplet state. The rate of change 
of these probabilities is given by: 


(1) dP, = — (3y7+As) PietrP 
(2) dP, = 3yP.—(y+aA,) P, 


where y is the triplet < singlet conversion rate, \, is the free space nnihila- 
tion rate for the singlet state, and \, is the free space three-quantum annihila- 
tion rate for the triplet state. The decay of positronium will then exhibit two 
lifetimes 71, 72 given by: 


(3) /t, = i = atv (a?—8?), 
(4) ‘tT. = de = a—V (a? 8), 
where 

2a = 4y+Ast+Ay, 


B? = y(As +31) +AAr 


Let us consider the 7» lifetime only. Using the free space values of A, and 
d,, the variation of tr. with y can be evaluated using (4). If it is now assumed 
that water is a sufficient converter to give the observed rz = 1.80 10~° sec, 
the partial conversion rate y; for water must be 7.0X10* sec~!. For a parti- 
cular quenching experiment with a given rz value, the total conversion rate y 
is given by (4) and the partial conversion rate, y2 for the added substance is 
then (y—vy1). The curve labelled conversion in Fig. 2 is a plot of yz = (y—71) 
vs. T2 normalized to the experimental results at rz = 1.0 10~® sec. The lower 
limit which this curve tends to is 5.0X10-" sec, four times the singlet life- 
time. The factor 4 arises from the fact that for rapid conversion the positron- 
ium atom spends 3/4 of its time in the triplet state, as originally pointed out 
by Garwin (1953). 

To obtain the PICKOFF curve of Fig. 2 it is assumed that the 7, lifetime in 
water is due to positronium decaying directly out of the triplet state by the 
pickoff mechanism, with a small contribution from the less probable three 
photon process. The \» corresponding to the observed r2 for pure water is then 
given by: 

(5) A2(H2O) = AV+A,_z 


where X, is the pickoff annihilation rate in water. The A, term in (3) and 
(4) must then be replaced by \,+A,. The A, term must also be replaced by 
A; +A, since pickoff can occur equally for the singlet state. Using these values 
in (1) to (4), the variation of 72. with y is obtained. In this case y = y2, and 
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this is plotted in Fig. 2 as the PICKOFF curve. It has also been normalized to 
the experimental points at ro = 1.0X10-* sec. The limit for this curve is 
given by: 


(6) 1/ro(min) = A2(max) = {(As+Ap) +3(A,+A,)} /4- 


This gives 72 (min) = 3.9X 107" sec, in good agreement with the experimental 
value of 4.0107" sec. The PICKOFF curve agrees better in shape with the 
experimental points than the CONVERSION curve owing to the lower final 
limit on 72, although there is still a discrepancy somewhat greater than the 
assigned error of the experimental points. Owing to the chemical complexity 
ot the solutions it would not be surprising to find minor differences between 
the experimental points and a curve calculated on such simple grounds. 
(Of course, when the FeCl; is added to the water, the value of A, will presum- 
ably change slightly. For the dilution used here, this effect will be small.) 
The agreement is good enough, however, to verify that we are witnessing a 
process of conversion by the FeCl. 

We have already assumed that the 7, lifetimes in amorphous substances 
are due to the pickoff process, and we therefore drop all further reference to 
the CONVERSION curve in Fig. 2. The close agreement between the experi- 
mental limiting lifetime, r2(min) = 4.0*10-"" sec, and its value calculated 
from (6), 3.9107" sec, shows that the use of the free space value of X, in 
(6) is justified. (The value of 72(min) calculated from (6) depends mainly on 
\,, partly on the experimental \,, and hardly at all on \,.) We conclude that 
the singlet positronium atom immersed in water has essentially its free space 
properties, at least so far as the electron density at the position of the 
positron is concerned, within an accuracy of perhaps 10%. Similar conclusions 
have been drawn about the triplet positronium atom immersed in teflon 
(Telegdi, Sens, Yovanovitch, and Warshaw 1956) and liquid helium (Paul 
and Grahara 1957). 

The results of similar experiments with other iron group elements are 
shown in Figs. 3, 4, 5. The FeCl; results are shown again in Fig. 3 as a com- 
parison with the other salts. The solid curves are drawn to fit the experi- 
menta! points and do not represent any particular theory. It is evident from 
these results that FeCl; is by far the strongest quencher of the iron group 
elements which have been studied. 


(b) Rare Earths 

The paramagnetism of the rare earths is due to the presence of unpaired 
electrons in the unfilled 4f subshell of these atoms. Only two of the rare earths 
have been studied to date, one of these being diamagnetic lanthanum, the 
first member of the group, which is not a quencher since there are no 4f elec- 
trons. The other substance studied was paramagnetic neodymium, with three 
unpaired 4f electrons. The results for NdCl; are shown in Fig. 6. It was not 
possible to make the NdCl; solution strong enough to follow the complete 
variation of rt. with concentration, but the results obtained show it to be a 


quencher like the elements of the iron group, but much weaker. 
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Fics. 3, 4, 5. Experimental values of 72 as a function of concentration for CoCle, CuCl, 
and FeCl; (Fig. 3); CrCl; and NiCl: (Fig. 4); FeCl, MnSO,, and MnCl, (Fig. 5). The curves 
are drawn through the points and do not represent any theory. 
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Fic. 6. Experimental values for 72 as a function of concentration for SbCl; and NdCl;. 
The curves are drawn through the points and do not represent any theory. 


The reason for the much weaker conversion by NdCl; is apparent when the 
electron configuration of the atom is considered. The 4f subshell lies deep 
in the structure of the atom, being surrounded by the filled 5s and 5p subshells. 
(The 5d and 6s electrons are assumed to be occupied in chemical combination 
with the chlorine.) The iron group elements, however, have their 3d electrons 
on the outside of the atom, once the valence electrons are taken away. Hence, 


the unpaired 3d electrons of the iron group are within easy reach of the posi- 
tronium atom, while the 4f electrons of the rare earths are protected by the 
outlying subshells. This explains why the conversion is more marked for the 
iron group elements than for the rare earths, even though the paramagnetic 
susceptibilities of the two groups are comparable. 


(c) Antimony 

An unexpected converter, antimony, was also discovered. The results for 
SbCl; are also shown in Fig. 6. The electronic configuration of antimony does 
not exhibit any obvious abnormality which would explain the observed 
effects. The ion Sb*** at first sight appears to have a closed 5s subshell, within 
which lie a vacant 4f subshell and filled 4d, 4p, 4s subshells. With this arrange- 
ment of electrons no quenching would be expected. 

A possible explanation might be the following. If it is assumed that the 
antimony exists as Sb+++, then when the 5p valence electrons are occupied 
(chemically) perhaps one or both of the 5s electrons fall back into the vacant 
4f shell, where they would be unpaired and capable of producing conversion. 
A similar sort of behavior has been noted for non-trivalent rare earth elements 
(Selwood 1943). For example, the normal rare earth ions Eut*+*+, Gdt+*+* have 
six and seven 4f electrons respectively and exhibit different paramagnetic 
susceptibilities. The divalent europium ion, Eut+, however, has a para- 
magnetic susceptibility almost identical with that of gadolinium, suggesting 
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that the 5d electron of Eut+* has dropped back into the 4f subshell to give a 
structure identical with that of Gdt++. This might also occur for the 5s 
electrons of antimony. 

Another possible explanation is that the antimony is not present in the 
solution as Sbt** at all, but rather as an ion in a different valence state, 
which might leave the antimony with one or more unpaired electrons. Sidgwick 
(1950) has suggested that an SbC1; solution contains much of the complex 
ion SbO*, so that the situation is perhaps rather more complicated than 
suggested here. 


(d) Cross Sections 

For all the quenchers found, cross sectiors for triplet — singlet conversion 
have been calculated. This was done assuming that all the atoms of the 
quenching agent are equally able to cause conversion. If the solutions are 
fully ionized this would be correct, and even if this were not so, little difference 
is expected for an unionized atom, since the valence electrons will be occupied 
in the formation of the compound and should not affect the conversion pro- 
duced by the unpaired electrons. 

Assuming a “‘pickoff plus conversion’? mechanism, the conversion rate y 
for triplet < singlet conversion for rz = 1.0X10~-* sec is obtained from an 
equation like (4) modified to include the A, term. This is combined with the 
molar concentration required to give rz = 1.010~-* sec. The cross section 
o is then given by 
(7) y = Naod 


where .V, is the number of atoms per cc, and @ is the Maxwell—Boltzmann 
average velocity of the positronium atom, assuming thermal equilibrium. 
This reduces to: 

(8) Ope = 1.21XK10—'9/ My cm? 


where My is the molarity required to give rz = 1.0X10~° sec. (li the ‘“‘com- 
plete conversion’”’ picture had been accepted, the result would have been 


(9) Seo = 2.84 10-"/ My cm, 


i.e. about double that calculated assuming pickoff plus conversion.) In Table 
I are shown the values for the cross sections of the quenchers found here 
using equation (8). Apart from assumptions about the mechanism, the error 
in the cross sections is about 10%. 
CONCLUSIONS 

The results in Table I show large variations in ¢,, for the ions studied. It 
can be said for the iron group that the quenching is definitely associated with 
the presence of unpaired electrons in the electronic structure of these atoms. 
This becomes evident when it is realized that a 4-molar solution of ZnCl, 
exhibited no quenching; the Zn**+ ion has a filled 3d shell with no unpaired 
electrons. 

There appears to be little relation between the number of unpaired electrons 
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VALUES OF THE CROSS SECTION dpc 


Ion Compound No. unpaired electrons Op-(cm?) X10" 


6.3 
2.5 
1.8 
38 
4.8 
5.3 
5.6 
11 
0.40 
3.7 


Crttt CrCl; 
Mntt MnCl, 
MnSO, 
Retr FeCl, 
Fett* FeCl. 
Cott CoClh 
Nitt NiCl, 
Cut* CuCh 
Ndt*++ NdCl; 
Sb**+*(?) SbCl; 


"UO mm DO Oo me ONO OO 


and the observed cross sections; e.g. Cut* with only one unpaired electron 
has a cross section second only to Fet+*+*, and is a much better quencher 
than other ions with two, three, four, and five unpaired electrons. Again, the 
apparent difference between the Mn** cross sections for the two different salts 
MnCl, and MnSO, is puzzling, since the only ion contributing to the con- 
version process in each case is Mn**. (Both NaCl and Na2SO, are inert in this 
respect.) The lack of correlation between the number of unpaired electrons 
and the observed conversion cross sections must be due at least in part to 
the complex chemical nature of the solutions. The formation of hydrates and 
complex ions would reduce the conversion effect by preventing close approach 
of the positronium atoms to the unpaired electrons of the paramagnetic ions. 
The low value for the cross section of Nd+** reported here is to be understood 
on similar grounds, for the outlying subshells of the atom surround the 4f 
subshell where the unpaired electrons lie. The result for SbCl;, which at first 
sight would not be expected to show any effect, is not fully explained, but 
two possible mechanisms were suggested earlier. 

The triplet — singlet conversion cross sections quoted here have all been 
derived on the assumption that the positronium atoms are fully thermalized. 
This step has been taken for definiteness, in spite of the fact that angular corre- 
lation measurements on water made recently by de Zafra (unpublished) do 
not in fact show the very narrow angular component that would be ex- 
pected on this assumption. Instead de Zafra’s results suggest positronium 
mean velocities a few times the thermal value, in agreement with Wallace’s 
estimate (1955) that positronium atoms may well not have time to thermalize 
completely within their lifetime. 

The mechanism of the triplet <> singlet conversion can be either that of 
electron spin exchange, or that of electron (or positron) spin flip by inter- 
action with the magnetic field of the paramagnetic ion. The latter process is 
not expected to furnish quenching cross sections so large as those observed 
here (Ore 1949). Conversion by electron exchange has been discussed recently 
by Ferrell (1958), who shows that conversion can take place whether or not 
electron spin flip is energetically possible. His argument presumably applies 
to the converters used here, because to flip the spin of a single one of the 
five unpaired Fet+* electrons, for example, requires an energy of many 
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electron volts. The cross sections found in this study for the iron group of 
paramagnetic ions range from 38X10~-' cm? for Fet*+ to 1.8X10- cm? 
for Mn++ from MnSQ,. These results are in rough average agreement with 
Deutsch’s value of 4X10~" cm? for the oxygen molecule. Deutsch’s value 
for nitric oxide is much larger, about 10~'* cm?; on the other hand Ferrell 
points out that Heinberg and Page (1957) seem to find oxygen a slightly better 
converter than nitric oxide. This point needs clarification. 

We note that Berko and Zuchelli (1956) computed their triplet — singlet 
conversion cross section for DPH in benzene on the basis that benzene is a 
converter, i.e. essentially from equation (9). On the pickoff view, we would 
use equation (8), reducing their cross section to about 610-8 cm?. In view 
of the results of de Zafra (1958) quoted earlier, this cross section may very 
well be meaningless, since the effect of DPH is probably not one of conversion 
at all. 

The disagreement in detail between the shape of the theoretical quenching 
curve and the experimental points (Fig. 2) could be explained as follows. As 
the concentration of the converter is raised, the mean lifetime of the positron- 
ium atoms becomes smaller. If these atoms are indeed not thermalized, but 
are still in the process of slowing down, then at higher concentrations the 
mean velocity of the positronium atoms is higher. The conversion rate varies 
as ov, and if o varies less rapidly than 1/v, as seems reasonable, the conversion 
rate will be relatively greater at higher concentrations. This will make our 
simple theoretical curve look insufficiently steep, as shown in Fig. 2. The other 


experimental curves (Figs. 3, 4, 5) show shapes very similar to that of FeCl; 
in Fig. 2. The above suggestion has the virtue of applying equally to all the 
converters, i.e. of depending on the conversion rate introduced into the 
solution, and not on the absolute concentration of the converting substance, 
and hence gives the observed similar shape for all. 


We wish to thank Professors J. S. Foster and P. R. Wallace for valuable 
discussions on different aspects of this work. One of us (R.E.G.) was supported 
by National Research Council Studentships. 


REFERENCES 

BELL, R. E. 1958. Bull. Am. Phys. Soc. Ser. II, 3 (3), 182. 
Bell, R. E. and Graham, R. L. 1953. Phys. Rev. 90, 644. 
BERKO, S. and HEREFORD, F. L. 1956. Rev. Modern Phys. 28, 299. 
BeErKOo, S. and ZucHELu, A. J. 1956. Phys. Rev. 102, 724. 
DeutscH, M. 1951. Phys. Rev. 82, 455. 
—— 1955. In Beta- and gamma-ray spectroscopy, edited by K. Siegbahn (North-Holland 

Pub. Co., Amsterdam). 
Deutscu, M. and Dutit, E. 1951. Phys. Rev. 84, 601. 
DE ZaFRA, R. 1958. Bull. Am. Phys. Soc. Ser. II, 3 (3), 229. 
Dixon, W. R. and Trarnor, L. E. H. 1955. Phys. Rev. 97, 733. 
FERRELL, R. A. 1958. Phys. Rev. 110, 1355. 
Garwin, R. L. 1953. Phys. Rev. 91, 1571. 
GREEN, R. E. and BELL, R. E. 1957. Can. J. Phys. 35, 398. 
HEINBERG, M. and Pace, L. A., 1957. Phys. Rev. 107, 1589. 
Ore, A. 1949. Univ. Bergen Arbok, Naturvitenskap. Rekke, No. 12. 
Paut, D. A. L. and Grauam, R. L. 1957. Phys. Rev. 106, 16. 
PeasLEE, D. 1950. Helv. Phys. Acta, 23, 845. 





GREEN AND BELL: TRIPLET POSITRONIUM 1695 


Ponp, T. A. and Dicker, R. H. 1952. Phys. Rev. 85, 489. 

SELWoop, P. W. 1943. Magnetochemistry (Interscience Publishers, Inc., New York). 

Sipcwick, N. W. 1950. Chemical elements and their compounds (Oxford University 
Press). Vol. I, p. 793. 

TELEGDI, V. L., SENs, J. C., YovANovitcH, D. D. and WarsHaw, S. D. 1956. Phys. Rev. 
104, 867. 

WaLvace, P. R. 1955. Phys. Rev. 100, 738. 

Yanc, C. N. 1950. Phys. Rev. 77, 242. 





NUMERICAL DETERMINATION OF THE FACTOR OF 
ENCIRCLED ENERGY RELATIVE TO A DIFFRACTION 
PATTERN OF REVOLUTION! 


Guy LANSRAUX AND GERMAIN BOIVIN? 


ABSTRACT 


An optical system which is symmetrical about some axis forms a diffraction 
pattern around its focus F in an observation plane II. 

The factor of encircled energy E(W) is, apart from a normalization factor, 
the energy spread over the II plane inside a circle of radius W, centered on F. 

Although it is a good test for the quality of an optical system, the factor of 
encircled energy has not been used a great deal because of the lengthy calcula- 
tions that its numerical determination requires. Two rapid computing methods 
are exposed here, with tables allowing a complete determination of E(W). 


INTRODUCTION 
Here we consider the most general optical system symmetrical about some 
axis. Besides allowing for spherical aberrations and focusing error the trans- 
parency will be allowed to have radial variations (Fig. 1). 


Fic.1. Schematic diagram of the formation of a diffraction pattern in an optical instrument. 


When the system has no spherical aberration the incident wave 2 from an 
axial object point // is transformed into a spherical wave S whose center of 
curvature is C. From the geometrical optics viewpoint the path of a light 
ray is HQOPC. 

But with the hypothesis of a more general system, the incident wave 2 
is transformed into an emergent wave & still possessing symmetry about the 
Z'Z axis but otherwise deformed in some arbitrary fashion. The chosen 
reference sphere S is then osculatory to 2 at the point O, and C becomes the 
paraxial focus or principal focus. The aberration due to the system is measured 
by the optical path difference PAJ taken along the radius vector CP. The 


corresponding phase is: 


¢(P) = (29r/A) PM 


‘Manuscript received May 5, 1958. 
Contribution from the Physics Department of Laval University, Quebec, Que 
*Holder of a National Research Council Studentship. 
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(A being the wavelength of light in the medium of propagation where the 
image is formed). 

When the optical system is perfect and exactly focused, the image of the 
object point H is the Airy pattern formed in the plane perpendicular to the 
optical axis at the point C. 

When the observation plane II is perpendicular to the axis of the system 
(aberrant or not) at some arbitrary point F, then it is convenient to indicate 
its position by 

a = Il (FC/OC), 
this being called the ‘‘focusing error’. 

Finally, a further generalization is obtained by allowing for a radially 
modulated amplitude on the emergent wave front 2. This will be described 
by the transmission coefficient 

ide amplitude at P 
TEP) « See 
amplitude at Q 

Amplitude modulation{has been*dealt with by various authors and is now a 
familiar concept. It can be achieved in ordinary instruments by incorporation 
of an amplitude filter, optically flat but having a non-uniform transmission 
(Lansraux 1946, 1947, 1948, 1952, 1953, 1956). 

The complex amplitude at a point P in the aperture may be written as: 

S(P) = T(P) e#*™, 
Since radial symmetry is assumed, this may also be expressed as a function 
simply of the radial co-ordinate x referred to the center O of the aperture and 
normalized so that x = | is the radius of the aperture; then 
f(x) = T(x) ec. 


Ina diffraction pattern whose radial co-ordinate is W, the complex amplitude 


r(w) = fre eflest+0@)] 7, (Wx) d(x’) 
where Jo is the Bessel function of the first kind and of order zero. This classical 
expression for ['(W) can be written more simply: 
r(wW) = fre) e* 7. (Wx) d(x’) 
with P(x) = ax?+ (x). 


I. EXPRESSIONS FOR INTENSITY AND FOR THE FACTOR OF ENCIRCLED 
ENERGY 
The intensity, in the diffraction pattern, is: 


el 
(3) I(W) = |r(W)|’ = T(x) e* Jo (Wx) d(x") 


el = 
(4) J(W) = | T(x) e™® Jy (Wx) d(x”) | T(x) © Jy (We) d(x’). 
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Changing variables x to s and x to ¢ respectively in the first and second integral 
of the product (4) allows J(W) to be written as a double integral 


(5) 1) = ff 1%) TO cos [V(s)— VO] Ju (Ws) Jo(WA) dls) de). 


The factor of encircled energy (Lansraux 1946, 1952, 1955) introduced as 
a test of the quality of an optical system is (apart from a normalization co- 
efficient .V) the energy spread over a circular domain of the diffraction pattern 
of radius W, and centered on the observation focus F. 


w 
(6) E(W) = Nx2n f I(W) W dW. 
0 
We choose the coefficient NV so that E() = 1. 


(7) E(o) = x2" f I(W) W dW = 1. 
0 


Applying Parseval’s theorem to Hankel’s integral (2) which connects 
T(x) e¥™ and I'(W) we have 


~ 1 1 
(8) J Ir(W)|? WdW = af T(x) e* |? x dx = af T’(x) x dx. 
0 0 0 


Hence 


(9) 


1 
co vis foe} i 1 2 
an [ I(W) W dw an [ IT(W)|’W dW sr | T(x) x dx 
0 0 0 


We also can give another physical meaning to E(W) by considering the 
image of a circular disk of radius W, iiluminated by incoherent light of uniform 
intensity and contrasting with a dark background. E(W) is then the contrast 
of the center of the image with respect to the background. 

Replacing J(W) by its expression (5) we get a triple integral in which the 
orders of integration may be interchanged according to Fubini’s theorem, 


(10) E(W) aan J If J T(s)T(t) cos [W(s) — (t)] 


x 10s) JoWwoa(sace) | W dW, 


(11) E(W) = 2xN J T(s)T(t) cos [p(s)— ¥()] 


x | J ” Jaws) J(WOW aw | d(s*)d(t’). 


But: 


Ww We es 7 
(12) 7 To ( Ws) Jo(Wt) W dw = W([sJi(Ws) Jo(Wt) — t Ji(Wt) Jo(Ws)] : 
0 


2 2 
s—t 
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Substituting this expression into (11) it follows that 


(18) EW) = 2eN Jf 76) T@ cos [¥(s)- VO] 


x WEE Las) Zo =P TS) d(s*)d(t"), 


(14) E(W) = 24N ff[rore cos [¥(s)—y(t)]D(W, s, t) ds dt, 


with 
D(W, s,t) = 4Wst [sJ:(Ws) AU) Jo(Ws)] | 
The domain of integration 0 for the double integral (14) is the square 
0<s<¢ land 0 <t <1. The diagonal s = ¢ is an axis of symmetry for the 
integrated differential. It divides the domain of integration 1 into two tri- 
angles [7 and Z] respectively characterized by s < t and s > t. If we choose 
one determination, for example s > ¢, we have: 


(15) E(W) = 2xN Sf =4nN Sf 
O A 


“aan . f= '7)7O costv(s)— VOW, sed. 


The function D(W, s, t) will be called prRECTRIX; it depends on three vari- 
ables W, s, t but involves only the products S = Ws and T = Wt. In fact, 
we can write 
(16) D(W,s,t) = D(Ws, Wt) 


_ 4 WsWt (Ws Ji(Ws) Jo(Wt)— Wt Ji(Wt) Jo(Ws)] 
7 (Ws)?— (Wt)? 


4ST [S Ji(S) Jo(T)—T I(T) Jo(S)] 


(17) D(W, s, t) = D(S, T) = SP : 


The function D(S, 7) is represented graphically by the directrix surface 
which is shown on Fig. 2. A table of D(S, T) = D(Ws, Wt) is given in Table I. 


II. NUMERICAL COMPUTATION OF E(W) 


For the calculation of E(W), the usual method consists of finding the com- 
plex amplitude of [(W) and from that the intensity for a sufficient number 
of values of W to ensure the desired accuracy for the integral 


Ww 
(18) E(W) = aan f I(W) W dW. 
0 
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LS 
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Fic. 2. View in perspective of the directrix surface D(S,,T)}= D(Ws4Wt). The thinner 
lines represent the sections S = constant, and T = constant over the range Ofto 25 steps of 0.5. 


TABLE I 
D(S, T) 





0.5 1.0 1.5 2.0 2.5 3.0 3.5 


+0.4697 +0.8541 +1.0851 +1.1257 +0.9768 +0.6762 +0.2901 
+1.5583 +1.9928 +2.0909 +1.8519 +1.3392 +0.6660 
+2.5797 +2.7634 +2.5371 +-1.9687 +1.1849 
+3.0619 +2.9693 +2.5329 +1.8535 
+3.1180 +2.9887 +2.6200 
+3. 
+4. 


IZ 


+3.2866 3809 
0026 . 
+5.3651 
+6.4488 
+6.9424 


Cobh OOM Ne HO 
|; COMOoOanonovon 





/ 


6.0 6.5 7.0 7.5 8.0 y R : 10.0 





—0.5291 —0.3002 —0.0193 +0.2463 +0.4369 51e -4643 . +0.0878 
—0.9165 —0.5544 —0.0917 +0.3578 +0.6912 8 -7781 s 
—1.0408 —0.7112 —0.2407 +0.2469 +0.6351 : - 8267 
—0.8182 —0.7122 —0.4399 —0.0920 +0.2371 46 .5596 
—0.2156 .4963 —0.6080 —0.5701 —0.4213 ; .0065 
-7426 .0106 .§6231 —1.0185 —1.1568 . .7234 
-9757 7737 .38515 —1.2258 —1.7229 .790: -4585 
3.3557 -8456 .3132 —0.9913 —1.8670 ‘ -9954 
-7230 +3.1429 -4135 —0.1808 —1.3930 d 2.1377 
.9080 5489 .9010 +1.2292 —0.2142 . 236% .7359 
+6.7552 .9021 -6302 +3.1301 +1.6140 . -7208 
+7.1459 -0173 5.3740 +5.2885 +3.8930 . . 8766 
-7162 +7.8609 +7.3844 +6.3105 : 2.9222 
+8.8268 +9.0697 +8.4956 +7. 5 . 1940 
+10.0383+ 10.0901 +9.176- -4145 
+ 10.8183 + 10.543 . 2929 ‘ 
+11. 5721 +9. 
+11.0708+ 10. 
+11. 


| 


ST OMMBNNAOUVUPRWOENNHKHO 
a 


oucnowe 


on 
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o 


on 
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TABLE I (Concluded) 





10.5 11.0 11.5 12.0 12.5 13.0 13.5 


—0.1429 —0.3287 —0.4280 —0.4211 —0.3140 —0.1360 +0.0677 
—0.2035 , —0.6939 .6958 —0.5306 —0.2447 +0.0889 ‘ “3ea2 
6834 -7174 —0.5755 -2999 +0.0373 352% 5698 
.3658 .4494 ; -2794 —0.0806 6 - 2964 
-1927 .0614 ° -1701 —0.2240 22% -1738 
.8418 -6920 A -0266 —0.3294 ‘ 5 -7076 
-3831 - 2705 P -2879 —0.3320 é - 1362 
-6227 -6163 : . —0.1913 . .3028 
-4259 -5868 “ . +0 .0883 6 .1100 
-7626 .1179 : . +0.4480 ‘ -5553 
-2705 - 2504 ’ .8139 +0.7820 . - 2573 
-4566 -8660 . . +0.9621 é 
-4942 -9923 : -0822 +0.8740 
3.0569 -8403 . , +0.4573 
.8639 —3.1283 . a —0.2612 
—1.7460 —2.6389 . . — 1.1522 
+0.3057 -2690 —2.1642 ; — 1.9927 : 
-7185 + +3.1231 -9389 -8219 ‘ —2.5011 .3945 
4360 +8.6915 +6.3778 +3.7917 . 2676 ; — 2.3920 - 1402 
+12.5588+ 11.5799 +9.6326 +6.9678 -9545 : —1, -1126 
+13.7194+ 13 .6865+ 12.4155+ 10.0654 9748 , +0. -0588 
+14.6553+ 14 .3049+ 12.6665 +9 .9819 5.6 +3. - 1069 
+15.0061+ 14.4062+ 12.5953 +9. 5 +6. +3 .2627 
+15.0344+ 14 .4589+ 12.7608+ 10. .0787 
+15.2999+ 14.9762+13. -0616 + 
+16.1439+ 16. -6346+ 12.071: 
+17. .2403+15.6046+12. 
+18 .4456+ 18 .0262+ 16. 
+ 18.9558+ 18.2: 
+19.0: 


I 





SCOCMBBNNAAMPeOWONNK HO 


CASUSMOMONOMOUOMOMOMNOUONOUOUOR 





16.0 17.0 2 20.0 21.0 22.0 23.0 24.0 25.0 


| 


+0.2971 —0.2823 575K .3381 +0.1903 +0.5219 +-0.3691 -1087 —0.4682 -0. 3911 
.2670 —0.0839 .3e .2635 +0.0431 +0.2931 +0.2627 .0057 —0.2561 .2617 
.1541 +0.4405 612 2417 —0.3176 —0.5668 —0.3044 -2131 5174 3492 
.6620 +0.6123 - 25 .7370 —0.4056 . 1207 7940 . 2282 9966 8342 
-7527 +0.0215 6 .6750 +0.0386 -5838 .6283 0955 4917 -5952 
.1265 —0.8581 , .0830 +0.6312 -7492 -2216 -4552 6936 3181 
.9269 —0.9992 . 86 .0194 +0.6487 .6193 .0899 3731 .4170 . 1390 
.5375 +0.0818 231k . 2692 —0.1602 9771 +1.1213 . 2356 7935 .0248 
.8409 +1.5447 . 898 .3476 . 1048 8513 .0658 .8145 -7996 .1167 
-1474 +1.7348 +2.69% . 2526 .0703 . 2162 .3327 - 6308 . 8830 . 3883 

3.0648 —0.2051 ; 2.1908 .3027 .5709 .8025 .3974 -2151 -5794 
.7428 —3.0086 -946 . 2499 -9185 -9145 -6230 .3731 8722 2711 

+1.3332 —3.6748 302% .4487 .9327 +3.1196 .5368 .1181 —2.5033 5978 

+8.5719 +0.3390 3624 3.9811 .4568 +2.6230 +2.8503 5760 8655 — 2.3361 

+16.0262 +8.7048 + : 3.5419 —3.6540 9546 -7049 +2.3269 9255 9435 

+19 .8458+ 17 .6739+ 10.3538 +1.6118 2254 -9504 7067 +2.1238 +3.4940 .7728 

+ 22.1892+19.3251+10.4900 +0.6161 1637 8961 7512 +3.0398 +3.5018 
+ 23 .0173+ 19 .3370+ 10.3707 +0.9794 . 3828 3416 9662 +2.1983 

+ 23 .5899 + 20.7123+4+ 11.9776 +1.8445 .7805 5569 -9156 

+25 .8911+422.5904+ 12.4006 +0.9311 9162 .7964 

+ 27 .0070+ 22 . 6858+ 12.0835 .9981 —5.2565 

+ 27 .3822+ 23 .7641+ 13.5672 +2.0272 

5602+ 25 8062+ 14.2903 

30 .9768+ 26 .0622 

+31 .2250 


| 
| 
| 


But this process is too long to be practical. A simpler approach is obtained 
by numerical integration of the expression 


easel 
(19) E(W) = inv | fir ) T(t) cos [¥(s)— ¥(t)] D(Ws, Wt) ds dt. 
s=() 

In particular we can use Cotes’ method with the values of the directrix func- 
tion D(S, 7) given in Table I. 

When the complex amplitude in the aperture can be represented (with a 
sufficient accuracy) by a complex polynomial of the fifth degree, the following 
method is much preferable. 
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We consider the six fixed abscissas 
(0) m=0, m=V/02,..., x =V(n—1)X02,..., x= 1. 


There is a unique set of polynomials II,(x) which take the value +1 for 
x = x, and vanish for any other abscissa of the chosen set. These are the 
Lagrange polynomials of the fifth degree relative to this chosen set of abscissas. 

Let F(x,) be the complex value of the function 


(21) F(x) = T(x) e#@ = T(x) etlaz*+9@1 


at the point whose radial co-ordinate is x, in the aperture; the polynomial 
n=6 


(22) P(x) = » Th, (x) F (xn) 


is then an approximate expression for F(x) which equals F(x) for any x = Xp. 
Inasmuch as the approximation made by substituting P(x) for F(x) agrees 
with the accuracy of the physical determination of “— we have 


n=6 


(23) T(W)~ f P(x) Jo(Wx)d(x*) = DO Fes) Mh, (xe) Jo(Wae)d (x) 


TABLE II 
COEFFICIENTS Kyo(W) 


| Kye. 5) Kpg(1.0)  Kpg(1.5)  Kyg(2.0) Kyg(2.5)  Kyg(3.0)  Kpg(3.5) 


+0.00027 +0.00107 +0.00234 +0.00403 +0.00605 +0.00834 +0.01081 +0.01342 
‘ +0.00420 +0.00919 +0.01571 +0.02333 ; 5 +0.03988 +0.04778 
+0.00259 +0.00511 +0.00751 +0.00898 . +0.00690 +0.00290 
+0.00274 +0.00580 +0.00948 +0.01328 : +0.01915 +0.02046 
+0.00379 +0.00725 +0.01012 +0.01124 , +0.00602 +0.00024 
+0.00096 +0.00183 +0.00257 +0.00289 ‘ 55 +0.00188 +0.00077 
+0.01654 +0.03608 +0.06126 +0.08995 ; +0.14759 +0.17152 
+0.01021 +0.02009 +0.02934 +0.03492 -03466 +0. ‘ae +0. 01592 
+0.01076 +0.02277 +0.03701 +0.05129 5 
+0.01494 +0.02847 +0.03958 +0.04381 
+0.00378 +0.00720 +0.01004 +0.01125 ; 
+0.00628 +0.01131 +0.01478 +0.01595 -016 +0.01784 +0.02405 
+0.00667 +0.01274 +0.01797 +0.02079 ‘ +0.01876 
+0.00921 +0.01611 +0.02030 +0.02126 -03 
+0.00234 +0.00408 +0.00515 +0.00541 
+0.01437 +0.02245 +0.02957 
+0.01807 +0.02438 +0.02654 
+0.00457 +0.00618 +0.00680 
; +0.02298 +0.02805 -+0.02890 x : 
+0.00101 a +0.00581 +0.00711 +0.00734 +0.00876 +0.01123 
+0. -00026 +0. 00086 +0. 00147 ~ 00180 —— 00186 +0. 00190 Bios 00215 +0. 00262 





| 
| AAMNAMAEAMNEWOAMNEONOUP WN | a 


p 
1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
3 
3 
3 
3 
4 
4 
4 
5 
5 
6 


Kyl 5) Kya. 0) Kyg(5.5 Kya(6. 0 Kyq(6.5)  Kpg(7.0) — Kpg(7. 5) Kyq(8.0) 


; 


a: 01612 +0.01885 +0. 02156 +0. 02421 ~ +0.02672 +0.02904 ~+0.03111 -+0.03288 
+0.06052 +0.06471 +0.06722 +0.06809 +0.06751 +0.06580 +0.06336 

—0.00991 —0.01767 —0.02547 —0.03270 —0.03888 —0.04370 —0.04702 

+0.01924 +0.01719 +0.01485 +0.01287 +0.01180 +0.01204 +0.01363 

.01253 —0.01702 —0.01905 —0.01842 —0.01558 —0.01145 —0.00718 

-00128 —0.00166 —0.00141 —0.00061 +0.00055 +0.00176 +0.00275 

- 20192 -20845 +0.21092 +0.21129 +0.21151 +0.21297 +0.21639 

-01379 —0.02572 —0.03314 —0.03565 —0.03414 —0.03020 —0.02566 

-07343 -07029 +0.06809 +0.06737 +0.06758 +0.06731 +0.06506 

-02162 —0.02864 —0.03070 —0.03061 —0.03142 —0.03488 —0.04080 

; 5 .00044 —0.00105 —0.00085 —0.00060 —0.00090 —0.00191 —0.00326 
+0.03642 -05471 .07684 +0.09988 +0.12066 +0.13714 +0.14837 +0.15467 
+0.01694 -02010 -02599 +0.03283 +0.03859 +0.04145 +0.04096 +0.03802 
+0 .04783 .06370 07647 +0.08247 +0.08068 +0.07316 +0.06359 +0.05561 
+0.01028 .01257 -01364 +0.01292 +0.01062 +0.00755 +0.00473 +0.00286 
+0.03796 -03851 .04003 +0.04216 +0.04370 +0.04423 +0.04433 +6.04590 
+0.01800 .02069 .02406 +0.02581 +0.02538 +0.02416 +0.02470 +0.02854 
+0.00516 .00553 -00580 +0.00560 +0.00497 +0.00445 +0.00459 +0.00547 
+0.06527 .07917 .08670 +0.08849 +0.08886 +0.09250 +0.10079 +0.11109 
+0.01407 .01609 +0.01675 +0.01662 +0.01686 +0.01831 +0.02073 +0.02310 
+0.00312 -00341 +0.00347 +0.00351 +0.00378 +0.00435 +0.00507 +0.00562 


GU nm ie i OS Go Go GS 8S NS NS ee ee oo 
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TABLE II (Concluded) 





| 
| 


Kpq(8.5) Kpq(9.0)  Kpg(9.5)  Kpg(10.0) Kpg(11.0) Kyg(12.0) Kyg(13.0) 


+0.03434 +0.03549 +0.03635 +0.03697 +0.03768 +0.03796 +0.03804 
+0.05777 +0.05521 J +0.04999 +0.04838 +0.04773 
—0.04955 —0.04918 a —0.04506 —0.04221 —0.04051 
+0.01952 +0.02270 +0.02530 +0.02777 +0.02716 +0.02561 
—0.00178 —0.00134 —0.00212 —0.00513 -—0.00720 -—0.00749 
+0.00343 +0.00316 +0.00268 +0.00174 +0.00140 +0.00149 
+0.22859 +0.23618 +0.24380 +0.25719 +0.26646 +0.27136 
—0.02059 —0.02172 —0.02547 —0.03912 -—0.05572 —0.06861 
+0.05199 +0.04256 +0.03343 +0.02287 +0.02671 +0.03864 
—0.05209 —0.05331 -—0.05051 —0.03712 —0.02524 —0.02320 
—0.00459 —0.00377 —0.00207 +0.00202 +0.00398 +0.00325 
+0.15757  +0.15782  +0.15974 +0.17427 +0.20454 +0.23875 
+0.03296 +0.03436 +0.03884 +0.04891 +0.04099 +0.00873 
+0.04998 +0.05004 +0.04854 +0.03433 +0.01319 +0.00841 
+0.00202 +0.00187 +0.00101 —0.00320 -—0.00656 —0.00444 
+0.06003 +0.07175 +0.08267 +0.09179 +0.09547 +0.12706 
+0.04023 +0.04176 +0.03845 +0.02848 +0.03314 +0.03639 
+0.00675 +0.00574 +0.00371 +0.00051 +0.00100 —0.00130 
+0.12262 +0.12300 +0.12420 +0.13815 +0.15337 +0.15542 
+0.02468 +0.02459 +0.02525 +0.02938 +0.03194 +0.03197 +0.03508 
+0.00587 +0.00597 +0.00636 +0.00761 +0.00807 +0.00829 +0.00926 


» 


| 
PER CORO | 
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Kpq(20.0) Kyq( ©) 


+0.03806 +0.03808 +0.03811 +0.03813 +0.03813 -03814 +0.0382195 
+0.04758 +0.04745 +0.04725 +0.04714 +0.04713 .04707 +0.0465499 
—0.04002 —0.03980 —0.03925 .03890 —0.03886 -03873 —0.0372963 
+0.02499 +0.02500 +0.02435 .02377 +0.02371 -02359 +0.0214509 
—0.00739 —0.00783 —0.00773 .00733 —0.00729 -00736 —0.0057307 
+0.00157 +0.00162 +0.00181 -00193 +0.00194 00203 +0.0027789 
+0.27293 +0.27363 +0.27502 -27584 +0.27593 . 27628 +0.2796442 
—0.07441 -—0.07571 —0.07944 -08206 —0.08288 -08321 — 0 .0925362 
+0.04765 +0.04772 +0.05213 -05645 +0.05709 -05783 + 0.0718682 
—0.02696 —0.02460 —0.02527 -02828 -02853 —0.0393788 
+0.00199 +0.00168 +0.00043 -00042 3 -00107 —0.0057307 
+0.26219 +0.26477 +0.27482 - 28317 ; - 28649 +0.3125251 
—0.03127 —0.03199 —0.04398 3 . .06270 —0.1024005 
+0.02943 +0.02573 +0.02776 , a -04022 +0.0718682 
+0.00114 +0.00180 +0.00518 : : -00941 +0.0214509 
+0.20109 +0.20318 +0.21779 . : - 25046 +0.3125251 
—0.00891 —0.01111 —0.01412 —0.03060 : 5 -03999 —0.0925362 
—0.01163 —0.01195 —0.01600 —0.02043 4 .02183 —0.0372963 
+0.18965 +0.20048 +0.20203 +0.21422 +0 .22847 +0.2796442 
+0.03828 +0.03756 +0.03828 +0.04136 +0.04067 +0.0465499 

+0.00974 +0.00993 +0.01106 +0.01193 +0.01204 +0.01296 +0.0382195 


| Kpg(15.0) Kyg(16.0) Kyg(17.0) Kyg(18.0) Kpg(19.0) 


-_ 


D Ot Co So Orme COND Crh Cot 





1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
3 
3 
3 
3 
4 
4 
4 
5 
5 
6 
pP @ 
1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
3 
3 
3 
3 
4 
4 
4 
5 
5 
6 


AACS Ore 





n=6 


T(W) =~ DOT (xn)e™ ff 1.0) 200%) de") 


n=1 


n=6 2 


I(W) = |TK(W)P?~ | DO Tlmne”™ J II, (x) Jo(Wx) d(x) 
p=6 g=6 


I(W)~ dX YX T(x,)T(x,) cos [¥(x,) — v(x,)] f Tl, (x)Jo( Wx) d(x’) 


p=ml gel 


x J “Tg(x) Jo(We) d(x’) 


ow 

(27) E(W) = ann f I(W)W dW 
0 
p=6 g=6 


(28) E(W)~4aN 3) DD T(xp)T (xq) cos [W(xp) — (x¢)]Kpg(W) 


p=1 g=l 
with 


K,,(W) = ; sl J (x) Jo( We) ac’) |] . T(x) Jo(Wx) ae?) |W aw, 


When the computation of the coefficients K,,(W) has been done once and 
for all, the calculation of E(W) is reduced to the quadratic form (28). The 
values of K,,(W) are given in Table II. 
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AKu(W) 


0} 


4K,.(W) 


0d 


Coefficients of Kpg(W). 
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4K(W) 


OF 


3B. Coefficients of Kpa(W) 
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4Ks(W AKs.(W) AKe(W) 
03 03 03 
2 02 02 
01 O01} 
w w w 
> O- > —_—_—$——— > 


Fic. 3C. Coefficients of Ky,(W). 


In the double sum (28) the two terms with subscripts pg and gp being 
equal, E(W) can be expressed by a reduced form containing 21 terms instead 
of 36. 

p=6 g=p—1 


(29) EW) ~ 4a | T° (xp)Knn(W) +2 ie YS T (xp) T (xq) 


n=1 =I q=1 


X cos [py (x,) — Hr) Kya) | 


In the particular case of a classical instrument of uniform transparency, 


7(x) = 1, the normalization factor (9) is N = i and we have 
Tv 
p=6 gq=6 
(30) E(W)~ >> D> cos [v(x,)—w(x_)] Kp, (w), 
p=1 g=l 
=§ g=p—l 
(31) E(W)~ . Kan(W)+2 : > cos [W(x,)—W(x,_)] Kpg(W). 
q=l 


But the Airy pattern formed by perfect systems (~(x) = 0) has a well- 
known factor of encircled energy 


(32) E(W) = 1-Je(W) —J2(W) 


where Jo and J; are respectively the Bessel functions of the first kind and of 
orders zero and one. This may be used to obtain a further simplification. 
From (31) we have 


p=6 g=p—l 


(33) E(w) = F Kan(W +220 2, Ky (W). 


Introducing 4 (W) into the expression (31) for E(W) the double sum is 
reduced to 15 terms 


p=6 gq=p—l 


(34) E(W) ~ &(W)-2 d, 2d {1 —cos [b(x,) — W(x_)]} Kpg(W). 


q=! 


This applies to an unfiltered aperture. 
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In the general case of a system with amplitude filter a similar reduction to 
15 terms can be obtained. If <°(W) is the factor of encircled energy for the 
diffraction pattern formed by a system without aberration, exactly focused 
and having a transmission 7\.) it is given by 


n=6 - p=6 g=p—1 
é(W) = tan u T?(X%n)Knn+2 2d 2 Te) Poe dKa(W) | 
oa >= = 
so that 
(35) 


p=6 g=p— 


E(W) ~ &(W)—8aN YY T(%p)T (xe) (1—cos [(xp) — ¥(%4)] Ke). 


& (W) is not known for the general case. In some particular cases which 
may be of interes., however, it should be possible to calculate it analytically. 


ERROR ON E(W) 


There are two sources of error in the method of calculation of E(1W); one is 
in the approximation of F(x) = T(x) e*@ by the polynomial P(x); the other 
is a purely numerical error from the summation (28) in which the accuracy 
of each term is limited by the error on K,,(W) of Table II. 

The substitution of P(x) for F(x) generates an error on E(W) whose order 
of magnitude can be obtained by the following calculation. We will use the 
notations ['(W), 7(W), E(W), N for the diffraction patterns corresponding 
to the pupil F(x) and y(W), 1(W), e(W), m for those corresponding to the 
pupil P(x); the subscript 1 and 2 will refer to the real and imaginary parts 
respectively of F(x), P(x), [(W), y(W). Finally & and & will be the maximum 
error on F(x) and on F(x), that is to say 


(36) | F(x) -~ P(x) | < &, | Fo(x) — P3(x) | < bo. 
Krom that 


| vl el 
(37) |T.(W)—v1(W)| = lJ F\(x) Jo(Wx) a(x’)— f P,(x)Jo(Wx) d(x’) | 
e/ 0 | 





fw- P,) Jo(Wx) d(x*) | < f \cr—Po| rows) d(x’) 
it el | 7 ; vl 
< J |Fi—P,| d(x") < af d(x*) = &. 
Similar calculations give 
(38) \P2(W) —y2o(W) | <é. 
But we have 
(39) 1(W) = Ti(W)+P3(W) and i(W) = yi(W)+72(W) 
(40) |Z(W)—i(W)| = [Pi—yal [Pittyi!+|P2—vyel|P2+y2 


< 2/Pi— yi] +2/la—ye}] < 2 
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and 


w ow 
(41) E(W) = sen f I(W)W dW, e(W) = 4x0 J 1(W)W dW 
Ww 
(42) E(W)—e(W) = arf (NI—n1) WdW 
0 
WwW ow 
E(W)—e(W) = 22 (v+n) (I-1)WdW+2r (vn) (+i) WdW 
0 0 


ow 
J (I-i) waw| -+2x\N—n 
0 





(43) |E(W)—e(W)| < 2x (N+n) 


x faa W aw| 


(44) 24(V+n) if (I- wa] < 2a(N+n) (E:+8:)W? ~ 49 N (E:4+8)W* 
| pw 
(45) 2e|N—ni| (+i) Wwaw| <29|N —n|[E(W)+e(W) ]~40|N—n|E(W). 
Hence the order of magnitude of AE(W) 
(46) |AE(W)| = |E(W)—e(W)|< 24(N +n) (1+) W?+22(N —n)[E(W) 
+e(W)] ~ 4a N (Ei +h) W?+40|N—n|E(W). 


This last expression vanishes for W = 0 and increases with W; it is no 
longer valid at least for |AE(W)| = 1 because (whatever may be F(x) and 
P(x)) we always have 0 < E(W) < land 0 < e(W) < 1 hence |AE(W)|<1. 

Generally, the error given by (46) is an upper limit which is rather too 





large. 

On the other hand, the error which is made on the double sum is related 
to the accuracy of the tabulated values K,,(W). Using Table III of |AK,,(W) |, 
the error on K,,(W), we have 


p=6 q=6 
(47) |AB(W)| <4av SY Te»)T (4) 008 [4% — Hes) II AK a) 


The reduced forms (34) and (35) consisting of 15 terms are not only simpler | 
to calculate, but the associated error of summation is appreciably reduced. 


TABLE III 
ORDER OF MAGNITUDE OF |AK,,(W)| 


a q 
Me 1 2 3 4 5 6 


1 0 .000,005 0 .000,007 0 .000,010 0.000,010 0.000,005 0.000,005 
2 0.000,078 0 .000,035 0.000,040 0 .000,025 0.000,012 
3 0.000,075 0. 000,085 0 .000,055 0.000,018 
4 0.000, 100 0 .000,065 0 .000,025 
5 0.000,040 0.000,015 

0. 000,008 
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NOTES 





OSCILLATIONS IN THE OMEGATRON! 


L. R. McNARRY 


Oscillations in the electron collector current of an omegatron were first 
reported by Redhead (1953). Pierce (1948) has discussed the effects of positive 
ions in electron beams and he shows that growing space charge waves of the 
order of the ion plasma and ion cyclotron frequencies are to be expected when 
the beam is confined by a magnetic field. At the magnetic field used in the 
omegatron (of the order of 3000 gauss) the ion cyclotron frequency for nitrogen 
would be near the ion plasma frequency. Redhead postulated that the low 
frequencies (of the order of 10 kc/sec) observed were due to beats between the 
ion plasma frequencies and the ion cyclotron frequencies. Experiments per- 
formed by varying the entrance and exit hole sizes and so varying the ion 
plasma frequencies, due to changes in electron beam densities, tended to 
confirm this hypothesis. It was found subsequentiy that when oscillations 
were present in the electron beam, similar oscillations could be observed on 
any electrode common to the interaction space. 

One of the primary limitations in the use of the omegatron is due to the 
presence of a background ion current (of the order of 10~- ampere) which 
may be present when no radio-frequency (r-f.) field is applied to the omegatron 
(McNarry 1956). This background ion current appears to be associated with 
the oscillations mentioned above. Under these conditions it may be postulated 
that ions formed in the immediate vicinity of the electron beam can arrive 
at the ion collector by two fundamentally different processes, first, in the 
classical manner where the c-bits expand continually with time, and second, 
by a mechanism which causes the ion orbits to drift across the magnetic field. 
In normal operation the time required for ions to reach the ion collector is 
several milliseconds, and since the drift mechanism might well be a relatively 
slow process, pulse measurements were carried out on a simple double r-f.- 
plate omegatron in an attempt to differentiate between the two processes. 

The electron beam was pulsed either ‘‘on”’ or ‘‘off’’ by means of a grid 
between the filament and the shield box. Observations were made with a syn- 
chronized oscilloscope at either the ion collector, the electron collector, or an 
r-f. plate. Pressure changes were made by introducing mercury vapor from 
a pellet of mercury through a greaseless stopcock. No radio-frequency energy 
was applied externally during any of the experiments. Figure 1 shows the 
effects of pressure and the build-up of oscillations when the beam was pulsed 
“on’”’ for 1 millisecond. For the bottom curve, the pressure was 610-7 mm 
Hg and the oscillation frequency 6 kc/sec. (The waveform is essentially a 
sawtooth.) The middle curve was obtained at a pressure of 7X10~-7 mm Hg 
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and the frequency increased to 16 kc/sec, while for the top curve the pressure 
was 4X10-§ mm Hg and the frequency had increased to about 700 kc/sec. 





Fic. 1. Effect of pressure on oscillations and background ion current when electron is 


pulsed ‘‘on” for 1 millisecond. 
Fic. 2. Decay in background ion current when electron beam is pulsed “‘off’’ for 1 milli- 


second. 


The decay in the steady background ion current to the r-f. plate when the 
beam current is pulsed ‘‘off’’ for 1 millisecond is shown in Fig. 2. Tne lower 
trace was obtained at a pressure of 3X10-’ mm Hg when the background 
ion current was zero. The upper trace was obtained at a pressure of 410-* 
mm Hg. The current sensitivity in these two figures was about 10-* ampere 
per large division. 

A slow change in pressure over several seconds produced the effect shown 
in Fig. 3. The oscilloscope was d-c. coupled to an r-f. plate and at the same 
time a-c. coupled to the electron collector. Thus it was possible to observe 
simultaneously the gradual build-up of oscillations and the associated increase 
in background ion current. The time scale is 1 second per large division. 
Three seconds after the pressure rise started, from 310-7 mm Hg, the oscil- 
lations abruptly ceased at a pressure of 2X10-* mm Hg, and the background 
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ion current dropped to zero. This phenomenon was quite reproducible. Figure 
3(b) is a record of the sine-wave oscillations on the electron beam at a pressure 
of 10-*§ mm Hg having a frequency of 660 kc/sec. 





: Fic. 3.(a) Effect of a slow change in pressure on electron beam oscillations and background 
ion current. 

(b) An expanded portion of the curve in (a) showing a sinusoidal oscillation at a fre- 
quency of 660 kc/sec. 


Positive ion current to the ion collector or an r-f. plate is only detectable 
on the oscilloscope when the oscillations are coherent. Random oscillations or 
noise (which may have amplitudes of several millivolts across a 1000-ohm 
resistor in the electron collector circuit) do not produce background ion current, 
but immediately the oscillations become coherent, through, say, a slight 
change in pressure, background ion current is observed. This phenomenon is 
not frequency dependent over the range from 1 kc/sec to 1.2 mc/sec. 

It was found that the time of the onset of instabilities was pressure depen- 
dent, and Fig. 4 illustrates the delay in the start of instabilities, observed 
at an r-f. plate as a function of pressure, after the electron beam had been 
pulsed ‘‘on’’. The minimum observed delay was 5 microseconds, which is in 
good agreement with Mihran’s (1956) observation of 4 microseconds in a 
travelling-wave tube structure. 
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IN| MICROSECONDS 


OELAY 





io? 2 3 5s 7 10°6 2 3 5 7 105 
PRESSURE IN mm Hg 


Fic. 4. Delay in start of instabilities in positive ion current as a function of pressure. 


In summarized form the following effects are observed: 

1. The minimum delay before the onset of instabilities is 5 microseconds. 

2. Electron current oscillations stabilize in about 50 microseconds. 

3. Electron current oscillations reach maximum amplitudes in about 80 
microseconds. 

4. Positive ion current starts at about 90-100 microseconds. 

5. Positive ion current reaches its maximum amplitude in about 1 milli- 

second. 

Probe measurements indicate that oscillations, in the range of 200 to 
700 kc/sec in the electron beam, violently disturb its profile at a low fre- 
quency of 50 to 100 c.p.s. The profile disturbances do not appear at the electron 
collector, though the high frequencies are observed there. This phenomenon 
may be similar to the radial beam explosions observed by Mihran (1956). 

While the original problem of the cause of background ion current is still 
unresolved, these and subsidiary experiments tend to disprove the ion-orbit 
drift hypothesis. 
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CORRECTIONS 
Vol. 35, 1957. 
Page 615. 3v1+2r2.+7; R(32) Read 6.595 instead of 6.495 
618. 2v, +2. +73 — v2 R(15) 4.784 4.769 
620. Tableau VI @1 1354.91 1354.94 
we 673.00 673.02 
w3 2396.49 2396.40 
X12 3.65 3.62 
621. Equation (3) WwW’ —(Wo... (Wo... 
623. Line 17 Wo 51.33 51.31 
624. 04°1 Wo 4958.64 4968.64 
629. Line 13 B'«, D' B’4, D4 
Line 14 Be. D’a B’:, D’« 
648. Reference 6 40, 339 40, 1950 


Vol. 35, 1957. Pages 1361-1379. In Table VIII the neutron-binding energies as determined 
by mass measurements of some of the hafnium isotopes were misquoted. The binding 
energies of nuclei formed by neutron capture in Hf'§, Hf"?, Hf, and Hf!” should read 
6.28, 7.55, 6.17, 7.32 Mev. respectively. The remarks in the final paragraph concerning 
the comparison between the y-ray energies and the binding energies determined from the 
mass measurements therefore do not apply. We are grateful to Dr. W. H. Johnson for 
bringing this error to our attention. 


Vol. 36, 1958. Page 710. The following correction was submitted by M. S. Sodha, Aerophysics 
Section, Propulsion and Fluids Research Department, Armour Research Foundation, 
Chicago 16, Illinois, concerning the paper Low field mobility of carriers in nondegenerage 
semiconductors. “Dr. R. Stratton of Manchester pointed out in a private communication 
that Table I was based on a consistent arithmetical error. Accordingly, the calculations 
were revised and the correct results are presented in the table below. The data notation 
and the relevant equations used are the same as those given in the paper. 





T (CK). Xo I (xo, 1) B/Bo x0? F (x0) 
800 0.4 0.782 0.0414 0.0021 
533 0.6 0.789 0.0296 0.0052 
400 0.8 0.798 0.0240 0.0099 
320 1.0 0.808 0.0210 0.0169 
267 1.2 0.819 0.0191 0.0266 
229 1.4 0.830 0.0181 0.0401 
200 1.6 0.842 0.0175 0.0578 





These results are in substantial agreement with experiment. From somewhat different 
considerations Stratton (private communication) arrives at values of 8/8» of the same 
order of magnitude as given in the above table.” 
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Vol. 36, 1958. Page 803. The following illustration should appear above (c) and (d) in Fig. 2. 
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